AD  748424 


21 


( 


AD 


FUNDAMENTAL  STUDIES 

OF 

SEMICONDUCTOR  HETEROEPITAXY 

FOURTH  SEMIANNUAL  REPORT 

R.  P.  Ruth,  A.  J.  Hughes,  J.  L.  Kenty,  H.  M.  Manasevit 
and  A.  C.  Thorsen 


ARPA  Support  Office 

Research,  Development,  Engineering  and  Missile  Systems  Laboratory 
United  States  Army  Missile  Command  AMSMI  RND 
Redstone  Arsenal 
Huntsville,  Alabama 
Contract  No.  DAAH01-70-C-1311 


Distribution  of  this  document 
is  unlimited. 


Sponsored  by: 

^  Advanced  Research  Projects  Agency 
ARPA  Order  No.  1585 


DIStRIBUrrON  STATEivtEirr  A 

Approved  for  public  release; 
Distribution  Unlimited 


Reproduced  by 


NATIONAL  TECHNICAL 
INFORMATION  SERVICE 

U  S  Deparfrnent  of  Commerce 
SprIngfiplH  VA  ??151 


Security  Cl«»»ificition 


DOCUMENT  CONTROL  DATA  •  R  &  D 

_ (Security  clmB9lllemtion  ol  body  ot  •6«rr«el  mnd  Indowtng  mnnoiotion  mutt  bo  mnttrmd  lA'iRn  tho  ovoroU  rtport  tt  OtOtmlUmd^ 


1  ORtGtNA  TING  ACTIVITY  fCofpor«f«  «utflor>  T2«,  R  C  POR  T  SE  C  U  R I  T  V  CLASSIFICATION 

North  American  Rockwell  Corn,  Electronics  Group  I  TTnrlfiRctifipH 
Research  and  Technology  Division 
Anaheim,  California 


3  REPORT  TITLE 


|26.  CROUP 


Fundamental  Studies  of  Semiconductor  Heteroepitaxy 


4  OEScnierivc  notes  <7>pa  oi  raporl  and  Ineluaira  dafaaj 

Semiannual  Report,  January  1972  through  June  1972 


s  Au  THORiSi  (FIrtt  iMSM,  middlb  Initimit  lott  nomo) 

Ralph  P.  Ruth,  A.  James  Hughes,  Joseph  L.  Kenty,  Harold  M  Manasevit, 
Arthur  C.  Thorsen 


S  REPORT  OATE 

July  1972 


ta.  coNTnacT  on  onaNT  no. 

DAAH01-70-C-1311 

6.  PnOJEC  T  NO. 


7a.  TOTAL  NO.  OF  PACES  176.  NO.  OF  PEFS 


25  .+  8  (Appendix) 


Sr.  ORIG«NATOR*t  REPORT  NUMSERIS| 


C 7 1-86. 4/501 


S6.  OTHER  REPORT  NOIS)  (Any  othof  numbmto  tfiRl  mmy  bo  momi0»od 
fills  topott) 


10  OISTRIRUTION  STATEMENT 

Distribution  of  this  document  is  unlimited. 

Details  of  illustrations  in 
this  document  may  be  better 
studied  on  microfiche 

II  Si;RPlEMENT  ARY  NOTES 

II.  ABST.R^CT  ......  .  .  ... 

12.  SPONSORING  MILITARY  ACTIVITY 

Advanced  Research  Projects 

Agency,  ARPA  Order  No.  1585 

Washington,  D.C. 

'  otO 

been  modeled  by  means  of  interatomic  Morse  potentials  and  computer  simulation,  with  only 
anion-anion  and  cation-anion  potentials  used  for  the  A£203'  (2)  The  electron-on-network 
technique  has  been  used  to  calculate  work  functions  and  surface  double-layer  potentials  of 
monovalent  metals. (3)  Theoretically  calculated  changes  in  mobility  and  resistivity  in  the  plane 
of  the  film  in  the  SiAj2203  system  caused  by  stress  effects  arising  at  the  interface  due  to 
differential  thermal  contraction  and  to  the  piezoresistance  effect  are  in  excellent  agreement 
with  experimental  results.  Theoretical  formulas  have  been  developed  for  Si  and  Ge  and  for 
A £903  and  MgA£204  substrates,  for  a  variety  of  crystallographic  orientations.  Experimental 
studies  have  been  concentrated  in  the  (221)  Si/(1122)  A^gOo  ^ad  the  (100)  Si/(0112) 
systems.  The  mobility  anisotropy  factor  was  found  to  be  m>out  39  percent  for  the  (221^  ^ane  [ 
and  about  9  percent  for  the  (100)  plane.  Results  indicate  that  (221)  Si  probably  exhibits 
hi^er  electron  mobilities  than  o&er  more  commonly  used  orientations. 

Attempts  have  continued  to  optimize  the  Si  deposition  process  for  growth  on~(1120)  and 
(0ll2)  A£203  surfaces.  The  effects  of  post-nucleation  annealing  have  been  examined,  but  no 
significant  improvement  in  Si  film  quality  has  been  demonstrated.  The  effects  on  film  growth 
of  gas-phase  etching  of  A  220^  surfaces  prior  to  deposition  have  been  evaluated  further.  No 
significant  advantage  to  the  Si  growth  process  has  been  observed  to  result  from  cooled 
reactor  chamber  walls  in  the  vertical  systems  used  in  this  work. 

Additional  reactant  analyses  have  been  carried  out  by  mass  spectrometric  techniques, 
with  significantly  impurity  concentrations  found  in  some  of  the  reactants  (especially  SiH^. 
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13.  ABSTRACT  (Coit) 

Study  of  the  chemistry  and  reaction  kinetics  of  CVD  processes  used  for  growing  Si  and  GaAs 
films  in  heteroepitaxial  systems  has  been  initiated. 

Evaluation  of  polishing  methods  for  MgAf204  surfaces  has  indicated  that  surface  fill-in 
occurs  for  this  material  just  aa  for  Af203;  gas-phase  etching  experiments  have  been 
initiated.  N-type  Si  films  with  electrical  properties  at  least  as  good  as  those  grown  on  Af203 
have  been  obtained  on  Czochralski-grown  stoichiometric  MgAf  2O4  when  He-H2  gas  mixtures 
are  used  for  the  gi'owth  environment;  autodoping  appears  to  be  operative  in  the  Si/MgA£204 
system  (in  He-H2  atmospheres)  at  approximately  toe  same  temperatures  as  for  Si/Af203. 

Ion-beam  sputtering  techniques  have  been  developed  for  preparing  ultra-thin  (-20oA) 
A2203  wafers  of  (0001),  (1014),  and  (0112)  orientations  as  substrates  in  the  to  situ  CVD 
experiments  with  Si.  Fabrication  of  the  modified-design  CVD  microchamber  is  essentially 
complete.  Numerous  in  situ  PVD  experiments  have  been  carried  out  with  All  and  Au 
deposited  onto  amorphous  carbon  substrates  to  delineate  experimental  techniques  required 
for  the  CVD  experiments.  Unusual  nucleation  effects  were  observed  in  Au  films  exposed  to 
toe  electron  beam. 

The  variation  of  the  electrical  properties  of  Si/A2203  with  temperature  has  been 
studied;  some  effects  attributed  to  high  defect  densities  or  inhomogeneous  strains  were 
observed.  Photoelectric  effects  in  the  Si,  A ^2^3  system  were  studied  further,  verifying 
that  the  observed  phenomena  result  from  photoinjection  of  electrons  from  metal  films  into 
the  Af203*  The  work  functions  of  several  metals  and  the  heights  of  the  metal-Af203  inter¬ 
face  barriers  were  determined.  High-field  transport  properties  of  heteroepitaxial  films 
have  been  continued. 

Work  on  the  determination  of  carrier  lifetimes  using  the  MOS  pulsed-capacitance 
technique  and  the  fabrication  of  Schottky-barrier  FET's  in  GaA8/Af203  has  continued.  An 
analysis  was  made  to  evaluate  the  effect  of  the  impurity  redistribution  in  the  Si  in  the 
region  near  the  oxide  interface  on  the  interpretation  of  lifetime  data  obtained  by  the  MOS 
technique. 

A  summary  of  the  work  planned  for  the  next  six  months  is  included. 
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■in.l  film  fT  f  program  is  to  carry  out  a  fundamental  study  of  nuclcation 
li  e  result  mechanisms  in  hetcroepitnxial  semiconductor  thin  films,  and  to  apply 

lubstrn cs  f  of  /mproved  films  and  thin-film  dciiccs  on  insulating 

substi.itcs.  Both  theoretical  and  experimental  investigations  are  involved,  with 

S.™on  o  '"’*7  ■'  a'’”''  ''Trl'”"  to  the  Sl-on-A^Os, 

5  y  ^  t’;'As-on-A«203  syslcms.  The  aeeomplishnients  of  the  fourth 

month  period  are  described  in  terms  of  seven  program  subtasks. 

Theoretical  studies  have  been  devoted  to  three  main  areas:  (1)  The  Si/A^oOo 

Suhur  w-fl!  "'T  interatomic  Morse  potentials  and  computer 

nhTT  ^  \  ^  anion-amon  and  cation-anion  potentials  used  for  the  AfoOi  and 

^nstraints  of  the  physical  lattice.  Investigation  of  surface 
r>)  Th?Hpl.T  basa  -plane  Ai’20;j  has  also  been  initiated,  using  these  potentials, 
snrf-m  tccliniquc  has  been  used  to  calculate  work  functions  and 

surface  double-layer  potentials  of  monovalent  metals.  (3)  Theoretically  calculated 
changes  in  mobility  and  resistivity  in  the  plane  of  the  film  in  the  Si/A^203  system 
and^i^ih^  stress  effects  arising  at  the  interface  due  to  differential  thermal  contraction 
and  to  the  piczoresistance  effect  are  in  excellent  agreement  with  experimental  results 
I  foimiulas  have  been  developed  for  Si  and  Ge  and  for  and  MgAf.d), 

substiates.  for  a  variety  of  crystallographic  orientations.  ^  ^ 

The  experimental  studies  of  anisotropy  of  electrical  properties  have  been  eon- 
centra  ed  in  the  (221)Si/  (1122)AM).j  and  the  (100)  Si/(0ll2)Af<>Oo  systems.  The 
mobility  anisotropy  factor  A  (defined  as  the  ratio  of  the  difference  between  the  maximum 
ml  minimum  values  of  carrier  mobility  m  the  plane  of  the  film  to  the  average  laUuHf 
Uk  mobility  in  that  plane)  was  found  to  be  about  39  jiereent  for  the  (221)  plane  and  about 
I  r  The  calculations  and  experimental  results  indicate  that 

oHentmion?*^^^  «xhjljits  higher  electron  mobilities  than  other  more  commonly  used 

-  ^^^‘-’•J^P^s  havc  continued  to  optimize  the  Si  deposition  jiroeess  for  growth  on 
-(1120)  and  (0112)  A^^G^  surfaces.  The  effects  of  post -nuclcation  annealing  have  been 
examined,  but  no  significant  improvement  in  Si  film  quality  has  been  demonstrated. 

The  effects  on  film  growth  of  gas-phase  eteliing  of  Ae20;.  surfaces  prior  to  deposition 
ha\e  been  evaluated  further.  No  significant  advantage  to  the  Si  growth  process  has 
been  observed  to  result  from  cooled  reactor  chamber  walls  in  the  vertical  systems 
used  in  this  woik.  Additional  reactant  analyses  have  been  carried  out  by  mass  speetro- 
metrie  techniques,  with  significant  impurity  concentrations  found  in  some  of  the 
reactants  (especially  Sillq).  Study  of  the  chemistry  and  reaction  kinetics  of  C\'1J 
pi  oeesses  used  for  growing  Si  and  GaAs  films  in  heteroepitaxial  systems  has  been 
*u*^*»fl  ^  ^  experiments  undertaken  arc  directed  tow'ard  determining  the  role  of 

the  surface  in  catalyzing  the  pyrolysis  of  Sillq. 

Kvaluation  of  polishing  methods  for  MgA?2G4  surfaces  has  indicated  that  surface 
M  material  just  as  for  Af203.  Gas-phase  etching  experiments  with 

MgAl'204  surfaces  have  been  initiated.  N-type  Si  films  with  electrical  proiierties  at 
least  as  good  as  those  grown  on  A?203  can  be  obtained  on  Czochralski-grown  stoiehio- 
metrie  spinel  (MgAK2G4)  when  IIe-Il9  gas  mixtures  are  used  for  the  growth  environment¬ 
autodoping  appears  to  be  operative  in  the  Si/MgA^204  system  (in  IIe-Il2  atmospheres)  ’ 
at  approximately  the  same  temperatures  as  for  Si/A%03. 


Ion-beam  sputtering  techniques  have  been  dev’elopecl  for  preparing  ultra -thin 
(-200  A)  A2203  wafers  for  use  as  substrates  in  the  in  situ  CV'D  experiments  with  Si. 
Wafers  of  (0001),  (1014),  and  (0112)  orientations  thinned  to  -0.002  in.  by  mechanical 
polishing  techniques  have  been  subsequently  thinned  by  ion  etching  to  the  point  of 
perforation  in  some  areas,  resulting  in  adjoining  regions  of  thicknesses  suitable  for 
transmission  electron  mieroseopy.  Fabrication  of  the  modified-design  CVI)  miero- 
chamber  is  essentially  complete,  with  testing  and  the  first  CVD  in  si^  experiments 
planned  for  early  in  the  third  year.  Numewms  in  si^  PVD  experiments  have  been 
carried  out  with  A2  and  Au  deposited  onto  amorphous  carbon  substrates  to  delineate 
experimental  techniques  and  problems  to  be  encountered  in  the  CV'D  experiments. 
Unusual  nuelcation  effects  were  observed  in  those  portions  of  the  Au  films  exposed 
to  the  electron  beam. 

The  variation  of  the  electrical  properties  of  Si/A2203  with  temperature  has  been 
studied;  some  effects  attributed  to  high  defect  densities  or  inhomogeneous  strains  were 
observed.  Photoelectric  effects  in  the  Si/A?203  system  were  studied  further,  verifying 
that  the  observed  phenomena  do  result  from  photoinjeetion  of  electrons  from  metal 
films  into  the  A22t)3-  The  work  functions  of  several  metals  and  the  heights  ofthemetal- 
A22O3  interface  barriers  were  determined.  The  high-field  transport  properties  of 
heteroepitaxial  films  have  been  further  studied. 

The  device  effort  has  centered  about  the  determination  of  carrier  lifetimes  using 
the  MOS  pulsed-eapaeitanee  technique  and  attempts  to  fabricate  Sehottky-barrier  FET's 
in  GaAs/A?203.  An  analysis  was  made  to  evaluate  the  effect  of  the  impurity  redistri¬ 
bution  in  the  Si  in  the  region  near  the  oxide  interface  on  the  interpretation  of  lifetime 
data  obtained  by  the  MOS  tcelinique. 

A  summary  of  the  work  planned  for  the  next  six  months  is  included. 
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SECTION  I 


INTRODUCTION 


This  is  the  fourth  Semiannual  Technical  Report  for  this  contract.  It  describes 
work  carried  out  during  the  period  1  January  -  30  June  1972.  Earlier  semiannual 
reports  (Refs  1-3)  described  work  done  in  the  first  18  months  of  the  program. 

1.  PROGRAM  OBJECTIVES 

The  overall  objective  of  the  program,  unchanged  from  that  originally  proposed, 
i^s  to  carry  out  a  fundamental  study  of  the  nucleation  and  film  growth  mechanisms  in 
heteroepitaxial  semiconductor  thin  films,  leading  to  new  knowledge  and  understanding 
of  these  processes,  and  then  to  apply  these  results  to  the  preparation  of  improved 
semiconductor  thin  films  and  thin-film  devices  on  insulating  substrates. 

The  specific  technical  objectives  of  the  three-year  program  are  the  following: 

1.  Investigation  of  the  many  aspects  of  the  mechanisms  of  heteroepitaxial 
film  growth,  to  establish  (through  accumulation  of  basic  knowledge)  sets 
of  techmcal  guidelines  for  the  preparation  of  better  films  which  can  then 
be  applied  to  real  situations. 

2.  Preparation  of  improved,  high-quality,  device-grade  heteroepitaxial  films 
of  Si  and  GaAs  on  insulating  substrates  by  chemical  vapor  deposition  (CVD) 
methods. 

3.  Development  of  methods  of  characterizing  heteroepitaxial  films  as  to  their 
suitability  for  subsequent  device  fabrication. 

4.  Design  and  fabrication  of  selected  thin-film  devices  which  take  advantage 
of  the  unique  properties  of  such  films. 

The  general  plan  for  accomplishing  these  objectives  involves  as  the  primary 
effort  the  study  of  the  fundamentals  of  heteroepitaxial  semiconductor  film  growth  on 
insulating  substrates.  Specialized  device  fabrication  is  used  both  as  a  means  of 
evaluating  certain  properties  of  the  films  (and  thus  as  a  measure  of  film  quality  as 
the  program  progresses)  and  as  a  means  of  exploiting  certain  unique  properties  of 
heteroepitaxial  semiconductor-insulator  systems.  The  determination  of  which 
fundamental  mechanisms,  properties,  and  processes  to  investigate  is  based  on  extensive 
background  knowledge  of  epitaxy  and  its  problems  and  from  the  thin-film  device 
difficulties  encountered  over  a  period  of  several  years  in  maiy  laboratories. 

The  problems  studied  are  in  no  way  restricted  to  those  identified  a  priori; 
experimental  (and  theoretical)  attention  can  be  shifted  as  the  program  progresses 
in  order  to  achieve  the  goal  of  a  better  understanding  of  heteroepitaxial  processes  and 
the  resultant  improvement  in  thin-film  active  semiconductor  devices. 


2.  PROGRAM  SCOPE 

The  program  involves  boih  theoretical  and  experimental  investigation  of  the 
nucleation  and  growth  mechanisms  of  heteroepitaxial  films  in  semiconductor-insulator 
systems,  the  development  of  improved  techniques  for  preparation  of  heteroepitaxial 
semiconductor  films,  and  the  fabrication  of  some  devices  utilizing  these  films, 
the  latter  primarily  for  the  purpose  of  evaluating  the  heteroepitaxial  film  materials 
but  also  to  exploit  the  special  properties  of  the  films. 

The  theoretical  studies  consist  of  two  types  of  activity.  First,  there  is  direct 
interaction  with  the  experimental  program  involving  data  analyses,  suggestion  of 
definitive  experiments,  and  postulation  of  specific  models  to  explain  experimental 
observations.  Second,  there  is  development  of  original  contributions  to  the  theory 
of  heteroepitaxial  growth,  the  goal  of  which  is  the  generation  of  significant  advances 
in  fundamental  epitaxy  theory. 

The  experimental  investigations  are  also  of  two  types.  First,  fundamental 
explorations  are  carried  out  to  delineate  mechanisms  and  general  empirical  principles 
of  the  heteroepitaxial  growth  process.  Second,  practical  studies  accompany  the 
fundamental  Investigations  so  that  useful  developments  can  be  immediately  applied  to 
the  improvement  of  semiconductor  films  and  thin-film  devices  on  insulating  substrates. 

The  work  has  emphasized  the  CVD  method  of  growing  semiconductor  thin  films 
because  of  its  importance  in  the  semiconductor  industry.  One  of  the  unique  aspects 
of  the  program  is  this  emphasis  on  the  study  of  fundamental  mechanisms  of  CVD 
growth;  most  previous  fundamental  studies  of  epitaxy  have  concentrated  upon  physical 
vapor  deposition  (PVD)  methods,  partly  because  such  studies  are  easier  with  PVD 
techniques. 

The  program  emphasis  is  on  filmc  of  Si  and  GaAs  and  substrates  of  sapphire 
(Af203)  and  spinel  (MgA?204);  not; stoichiometric  spinel  and  beryllla  (BeO)  may  also 
be  included  as  substrate  materials  later  in  the  program.  The  initial  emphasis  has 
been  on  the  Si-on-A2203  system,  with  increasing  attention  being  given  to  the  Si-on- 
MgA£204  GaAs-on-A2203  systems.  Si  and  GaAs  have  been  chosen  because  of 
the  preeminence  of  the  former  in  the  semiconductor  industry  and  the  high-frequency 
and  high-temperature  attributes  of  the  latter;  in  addition,  they  represent  the  elemental 
and  compound  semiconductors  for  which  most  comparative  information  exists. 

The  program  as  described  is  carried  on  primarily  at  facilities  of  the 
Electronics  Group  of  North  American  Rockwell  Corporation  (NR)  by  NR  personnel. 
Parts  of  several  of  the  subtasks  are  performed  by  personnel  of  the  University  of 
California  at  Los  Angeles  (UCLA),  in  the  Department  of  Electrical  Sciences  and 
Engineering  and  in  the  Chemistry  Department.  In  addition,  work  on  one  of  the 
subtasks  is  being  done  in  part  in  the  Department  of  Chemistry  of  California  State 
University,  San  Diego  (CSUSD).  Both  the  ULCA  and  the  CSUSD  programs  are 
supported  by  subcontracts  from  NR. 
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3.  PROGRAM  DESCRIPTION  BY  SUBTASK 

The  three-year  program  was  originally  divided  into  nine  subtasks  -  hvo 
theoretical  and  seven  experimental  (Refs  1,  2).  However,  at  the  start  of  the  second 
year  it  was  decided  that,  on  the  basis  of  the  way  in  which  the  work  of  the  first  twelve 
months  had  develojjed,  the  contract  work  would  be  more  accurately  described  in  terms 
of  seven  main  subtasks  -  one  theoretical  and  six  experimental. 

The  seven  subtasks  are  as  follows: 

Subtask  1:  Theory  of  Epitaxy  and  Heteroepitaxia!  Interfaces.  Theoretical 
examination  of  CVD  kinetics  and  the  processes  of  nucleation,  surface  migration,  and 
film  growth  with  emphasis  on  crystallographic  relationships  between  overgrowth  and 
substrate  to  attempt  to  identify  mechanisms  and  establish  general  principles  of 
heteroepitaxial  growth;  theoretical  modeling  of  the  heteroepitaxial  interface  using 
appropriate  potentials  to  determine  surface  configurations  and  interfacial  binding 
energies  in  real  and/or  simplified  systems. 

Subtask  2:  Deposition  Studies  and  Film  Preparation.  Investigation  of  the  effects 
of  various  experimental  parameters  upon  the  properties  of  deposited  semiconductor 
films;  investigation  of  the  delivery  kinetics  of  CVD  processes,  to  improve  the  detailed 
understanding  and  control  of  the  chemical  reactions  involved  in  the  preparation  of 
heteroepitaxial  semiconductor  films  by  CVD;  preparation  of  films  for  use  in  other 
parts  of  the  program. 

Subtask  3:  Analysis  and  Purification  of  CVD  Reactants.  Analysis  of  impurity 
content  of  reactant  materials  used  in  metalorganic-hydride  and  other  CVD  processes; 
preparation  of  research-sample  quantities  of  improved-purity  reactants  for  use  in 
film  growth  experiments. 

Subtask  4:  Preparation  and  Characterization  of  Substrates.  Preparation  of 
substrate  wafers  and  characterization  of  surfaces  and  impurity  content  of  substrates 
used  for  semiconductor  heteroepitaxy;  development  of  reproducible  new  and/or 
improved  substrate  polishing,  cleaning,  and  handling  methods. 

Subtask  5:  Studies  of  ip  aitu  Film  Growth  in  the  Electron  Microscope.  In  situ 
observation  and  study  of  the  early  stages  of  growth  of  CVD  films  in  the  electron 
microscope,  to  develop  additional  fundamental  knowledge  of  the  epitaxy  process. 
Results  of  these  experimental  observations  will  be  incorporated  into  the  theoretical 
studies  wherever  possible. 

Subtask  6:  Evaluation  of  Film  Properties.  Measurement  of  the  electrical, 
optical,  crystallographic,  and  thermal  properties  of  heteroepitaxial  semiconductor 
films  on  insulators,  by  a  variety  of  measurement  techniques.  Standard  techniques 
will  be  eniployed  and  new  methods  developed  where  required  for  measurement  of 
those  film  properties  which  appear  best  to  characterize  ultimate  device  performance. 
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Subtask  7:  Design  and  Fabrication  of  Devices.  Design  and  experimental 
fabrication  of  certain  types  of  devices,  using  heteroepitaxial  films  produced  in  the 
above  studies.  Some  devices  will  be  used  to  evaluate  material  properties  and  others 
to  exploit  semiconductor  film  characteristics  unique  to  heteroepitaxial  systems. 

These  subtasks  will  be  modified  as  needed  as  the  program  progresses. 

The  results  obtained  during  the  fourth  six-month  period  of  the  contract  are 
discussed  by  subtask  in  Section  11.  An  outline  of  the  work  planned  for  the  next 
six  months  is  contained  in  Section  111.  A  program  summary  to  date  is  given  in 


contains  a  pre-publication  version  of  a  paper  based  on  theoretical 
studies  at  UCLA  which  were  supported  in  part  by  this  contract,  as  part  of  the  effort 
on  Subtask  1. 
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SECTION  II 


RESULTS  AND  DISCUSSION 


The  work  of  the  fourth  six-month  period  of  the  contract  is  discussed  in  this 
section.  The  effort  of  each  of  the  seven  subtasks  is  described  separately.  Because 
of  the  interrelationships  among  the  various  subtasks  the  discussion  in  several  instances 
necessarily  involves  activities  and  results  which,  primarily  are  part  of  another  subtask. 

1.  SUBTASK  1.  THEORY  OF  EPITAXY  AND  HE  TE  ROE  PI  TAXI  A  L  INTERFACES 

Several  possible  approaches  to  the  theoretical  modeling  of  heteroepitaxial 
systems  have  been  investigated  to  date  under  this  contract,  The  general  criteria 
originally  adopted  for  determining  suitability  of  a  given  technique  are  that  the 
theoretical  treatment  and  associated  calculations  should  (1)  relate  explicitly  to 
heteroepitaxy;  (2)  be  as  nearly  as  possible  a  "first-principles"  approach;  (3)  relate 
as  closely  as  possible  to  an  actual  system  such  as  Si/A^203;  and  (4)  represent  an 
original  contribution  to  the  theory  of  heteroepitaxy. 

Attempting  to  meet  these  criteria  is  both  difficult  and  ambitious.  However, 
the  criteria  are  important  in  providing  a  general  framework  and  goals  for  the 
theoretical  studies.  Although  the  Si/A£202  system  is  of  considerable  commercial 
and  practical  importance,  theories  and  calculations  available  in  the  literature  fall 
drastically  short  in  that  there  is  no  real  connection  with  systems  such  as  Si/AioOo. 

The  emphasis  of  previous  theories  has  instead  been  on  simple  and  unrealistic 
one-  or  two-dimensional  cubic  lattices  which  unfortunately  cannot  be  applied  to 
Si/A£20Q«  The  intent  of  the  theoretical  studies  of  this  contract  is  thus  to  develop  an 
approach  which  relates  specifically  to  the  Si  and  the  A£203  lattices. 

During  the  first  year  of  the  contract  a  formal  theoretical  method  of  replacing 
overgrowth  atoms  on  a  substrate  with  Gaussian  mass  distributions  was  further 
developed  for  those  cases  where  the  effective  interatomic  potential  is  known.  The 
technique,  applicable  to  irregular-shaped  islands  or  films  of  finite  extent,  was 
applied  to  a  simplified  model  to  determine  preferred  orientation  relationships  from 
calculated  film-substrate  interaction  energies.  The  method  was  not  pursued  further, 
however,  because  it  was  not  sufficiently  adaptable  to  real  systems.  Several  other 
possible  approaches  to  the  theoretical  modeling  of  heteroepitaxial  systems  were 
critically  reviewed,  including  the  Frank-Van  der  Merwe  model,  a  Green's  function/ 
Wannier-function  approach,  a  contrived  potential-energy  model,  and  the  two-body 
interatomic  potential  method.  It  was  concluded  at  that  time  that  most  existing 
theories  are  inadequate  for  application  to  real  systems. 

The  feasibility  of  a  molecular  orbital  development  of  the  heteroepitaxial  interface 
was  then  investigated.  However,  early  in  the  second  year  it  was  determined  infeasible 
to  apply  this  technique  in  a  manner  directly  relevant  to  heteroepitaxy,  so  this  effort 
was  terminated.  The  interatomic  potential  approach  to  heteroepitaxy  was  then  re¬ 
instigated,  with  the  goal  being  the  computer  simulation  of  growth  of  Si  on  fsjt20^. 
Mechanical  stability  conditions  for  an  A^oOg  lattice  modeled  with  two-body  potentials 
were  investigated  and  determined  to  the  depth  required  for  these  applications. 

Computer  programming  of  the  Af20g  lattice  energy  and  elastic  constants  was  begun 
for  use  in  determining  appropriate  empirical  potentials  required  for  modeling  this 
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lattice.  Successful  modeling  of  the  AfoOo  lattice  is 
Si  growth  on  A 1203. 


a  major  requirement  for  modeling 


In  addition,  dunng  the  first  half  of  the  second  year  the  application  01  ihp  electron- 
on-network  theory  to  the  problem  of  determining  surface  configurations  and  interfacial 
billing  energies  in  heteroepitaxial  systems  where  the  surface  structure  is  allowed  to 
relax  was  investigated  and  initially  appeared  promising  for  real  systems  of  interest 
Normalized  eigenvectors  were  developed  as  a  basis  for  a  secular  equation  whose  ’ 
solution  IS  fundamental  to  the  solution  of  the  total  problem. 


During  the  past  six  months,  the  theoretical  studies  at  NR  aixl  UCLA  have  been 
devoted  to  three  main  areas:  (1)  modeling  of  the  Si/A^gOo  system  by  means  of 
interatomic  potentials  and  computer  simulation;  (2)  use  of  the  electron-on-network 
techmque  to  calculate  work  functions  and  surface  double-layer  potentials  of  mono- 
valent  metals;  and  (3)  calculations  relating  the  mobility  of  Si  films  on  AgoOo  to  stress 
effects  arising  at  the  heteroepitaxial  interface  and  caused  by  the  differential  thermal 
contraction  of  film  and  substrate.  These  investigations  will  be  discussed  in  sequence 

in  rhiQ  Q<iori/\rt  ~ 


a.  Studies  of  the  Si/Af203  System  Using  Interatomic  Potentials 


It  now  appears  that  a  general  solution  to  the  heteroepitaxial  Si/A^oOn  problem 
cannot  be  hoped  for  at  this  time  and  may  indeed  be  several  years  in  themmre  What 
can  be  done  is  to  incorporate  the  geometrical  and  structural  features  of  the  Si  and 
f  crystal  lattices  into  a  simplified  model  based  on  interatomic  potential  modeling 
ot  Si  and  of  A2203.  This  should  represent  a  reasonable  way  of  incorporating  those 
structural  features  of  Si  growth  on  A^203  which  are  caused  by  the  A^2^3  lattice. 


fi,  «•  interatomic  potentials  to  represent  the  interaction  between  atoms  in 

the  Si/Ai;20g  system,  the  cohesive  energy  of  a  given  configuration  is  obtained  by 
summing  over  all  interaction  atom  pairs.  As  background,  a  few  comments  ort  the 
theoretical  foundations  of  interatomic  potentials  will  be  included  here. 


^e  concept  of  an  interatomic  potential  in  the  present  context  is  best  regarded 
as  a  phenomenological  construct.  While  there  are  quantum  mechanical  foundations 
which  provide  a  basis  for  interatomic  potentials  at  the  philosophical  or  qualitative  level 
there  do  not  appear  to  be  ai^  quantitative  and  rigorous  calculations  of  interatomic 
potentials,  with  the  possible  exception  of  the  case  of  Na.  Even  if  rigorous  atomic 
potentials  existed  and  were  available  for  the  bulk  "perfect-crystal"  materials  of 
interest,  these  potentials  would  be  specific  to  the  bulk  crystal  and  would  have  to  be 
altered  in  the  presence  of  defects  or  at  a  heteroepitaxial  interface. 


interatomic  potentials  to  be  useful  in  practice  it  is  necessary  to  assume 
that  the  basic  interatomic  interaction  is  unaffected  by  atomic  rearrangement.  ITiis 
assumption  is  therefore  fundamental  in  the  present  studies.  A  second  consideration 
lies  in  the  particular  type  of  potential  to  be  employed.  These  investigations  assume 
a  simple  two-body  central-force  interaction  between  atoms.  More  complicated  non¬ 
central  potentials  could  be  employed,  in  principle,  but  in  practice  they  are  out  of  the 
question,  particularly  for  the  complicated  lattice.  In  this  work  the  atomic 

potential  modeling  of  Si  and  A2203  is  regarded  simply  as  a  phenomenological  treatment. 
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There  are  two  attractive  features  which  argue  for  application  of  atomic  potentials 
to  the  Si/Af203  system:  (1)  Computer  modeling  with  interatomic  potentials  has 
apparently  been  successful  in  the  past  for  structural  problems  in  simple  metals.  (2) 
There  is  at  present  no  theory  of  epitaxy  which  makes  it  possible  to  predict  the  orienta¬ 
tion  of  the  overgrowth  on  the  substrate  in  the  Sl/Af203  system.  The  orientation  aspects 
of  heteroepitaxy  can  be  treated  fairly  directly  in  terms  of  interatomic  potentials. 

The  theoretical  modeling  of  the  Sl/Afi203  system  with  interatomic  potentials 
can  be  divided  into  five  parts:  (1)  the  modeling  of  the  A^-Og  lattice  by  determination 
of  phenomenological  potentials  meeting  appropriate  crlteiia  or  constraints;  (2)  the 
investigation  of  AfgOg  surface  reconstruction  on  the  basis  of  the  determined  potentials; 
(3)  modeling  of  the  Si  lattice  with  an  interatomic  potential;  (4)  the  parameterization,  in 
lieu  of  useful  empirical  data,  of  the  Si/A^20g  adatom  interaction;  and  (5)  the  simulation 
of  Si  growth  on  AfgOg  by  computer  computation  of  energy  for  small  clusters  of  Si  atoms. 

The  modeling  of  the  Af20g  lattice  has  been  carried  out  in  terms  of  Morse 
potentials.  This  choice  was  suggested  by  the  following  factors:  (l)The  Si  is  to  be 
modeled  by  means  of  a  Morse  potential  available  in  the  literature  (Ref  4).  (2)  The 
Si/A^20g  adatom  potential,  which  is  to  be  treated  essentially  as  a  parameter,  can  be 
conveniently  expressed  with  Morse  potentials.  (3)  Analytical  expressions  for  the 
atomic  force  can  be  easily  obtained  for  Morse  potentials. 

^2s  two  types  of  atoms,  would  in  general  require  three  Morse 

potentials: 

1.  D(,(,  (exp(-2<.(,p 

2.  {exp(.2.^(,  (R-Rj^(,)) -2  exp  (- a^p(R-R^p)), 

3.  {exp  (-2  (R-R^^))  -2  exp  (-  »aa<«-«aa»» 

The  subscripts  CC,  AC,  and  AA  denote  cation-cation,  anion-cation  and  anion- 
anion  interactions,  respectively.  A  general  determination  of  these  three  interatomic 
potentials  would  require  treatment  of  a  nine- parameter  system.  A  straightforward 
determination  is  therefore  impractical.  It  thus  becomes  important  to  incorporate  as 
much  physical  and  empirical  data  as  possible  in  order  to  simplify  the  problem. 

The  appropriate  cohesive  energy  (energy  of  formation),  elastic  constants,  and 
lattice  parameters  are  known  experimentally.  Gieske  and  Barsch  (Ref  5),  in 
considering  the  elastic  constants,  have  noted  that  is  somewhat  unusual  compared 

with  other  oxides  such  as  MgO  in  that  the  Cauchy  relations  Cj2=Cj0  and  C23=C44  are 
approximately  fulfilled.  In  addition,  although  A^2®3  is  trigonal  (rhombohedral),  the 
second-order  elastic  constant  data  indicate  that  the  deviation  from  pure  hexagonal 
symmetry  (C44=0)  is  small.  (C44  is  about  an  order  of  magnitude  smaller  than  the 
other  elastic  constants.)  This  suggests  that  the  elastic  properties  are  related 
primarily  to  the  nearly-hcp  oxygen  fr.amework. 

Anderson  (Ref  6)  earlier  expressed  this  same  point  of  view  -  that  the  oxygen 
framework  is  dominant  in  determining  elastic  properties  of  tliis  oxide.  Since  the  M 
atom  sites  in  Af203  are  believed  to  play  an  important  role  in  heteroepitaxial 
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film-substrate  orientation,  8  model  based  solely  on  anion-anion  interactions  would 
not  be  reasonable.  Intuitively,  the  deviation  from  hep  symmetry  would  be  exhibited 
most  strongly  through  the  cation-cation  interaction,  which  relates  explicitly  to  a  non- 
hexagonal  lattice  with  defects.  It  is  therefore  concluded  that  the  cation-cation 
interactions  do  not  play  an  important  role  in  the  lattice  and  that  the  exclusion  of  cation- 
cation  interactions  is  a  reasonable  and  Important  simplification.  Consequently,  only 
anion-anion  and  cation-anion  potentials  have  been  employed  throughout  the  modeling. 

The  phenomenological  potentials  have  been  determined  by  requiring  that  the 
modeled  lattice  provide  (1)  the  correct  cohesive  energy  of  about  34.54  eV  per 
rhombohedral  unit  cell;  (2)  that  the  lattice  energy  be  a  local  maximum  and  the 
derivative  vanish  at  the  observed  lattice  spacing;  (3)  that  the  forces  on  each  of  the 
cations  in  the  unit  cell  vanish;  and  (4)  that  the  forces  on  each  of  the  anions  in  the  unit 
cell  also  vanish. 

Potentials  meeting  these  four  constraints  have  been  determined.  To  facilitate 
more  detailed  discussion  of  the  procedure  used  in  establishing  these  potentials  the 
first  few  interatomic  distances  in  the  ^2^2  lattice  (anion-anion,  cation-anion,  and 
cation-cation  distances)  arc  given  in  Table  I.  The  first  four  anion-anion  distances 
and  the  first  two  cation-anion  distances  relate  to  what  would  be  first- neighbor  inter¬ 
actions  in  an  idealized  Afo^i  lattice. 

2  *j 


TABLE  1. 

Interatomic  Distances  for 

Lattice  (A) 

Anion- 

Cation- 

Cation* 

Anion 

Anion 

Cation 

2.524 

1.857 

2.649 

2.619 

1.969 

2.792 

2.725 

3.222 

2.869 

3.223 

3.499 

3.430 

3.843 

3.790 

3.560 

4.762 

3.957 

3.862 

5.148 

4.334 

4.162 

5.449 

4.325 

4.616 

4.380 

4.762 

4.429 

4.813 

4.675 

4.945 

5.064 

5.072 

— 

5.244 

5.128 

5.366 

5.252 

As  a  first  attempt,  only  the  four  first-neighbor  anion-anion  and  two  first- 
neighbor  cation-anion  interactions  were  included  in  the  analytical  expression  for 
cation  and  anion  forces,  potentials  were  determined  which  gave  the  required  zero 
forces.  However,  examination  of  the  derivative  of  the  cell  energy  showed  that  a 
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maximum  could  never  be  obtained  with  only  first  neighbors.  All  subsequent  discussion 
therefore  relates  to  a  model  including  13  different  anion-anion  separations  (up  to  5.072  A) 
and  12  different  cation-anion  separations  (up  to  5.064  A).  This  cutoff  point  in  the  model 
is  arbitrary  but  is  believed  to  be  reasonable  for  present  purposes. 

Potentials  meeting  the  four  constraints  and  including  the  interactions  listed  above 
were  then  determined  for  a  range  of  values  of  the  potential  parameters  and  R^a* 

It  was  originally  planned  that  the  remaining  parameters  Rac  ^^d  Raa  would  be 
further  restricted  by  calculation  of  elastic  constants.  The  Af203  lattice,  viewed  as 
rhombohedral,  contains  10  atoms  per  unit  cell  and  involves  internal  (unequal)  displace¬ 
ments  or  strains  of  the  various  atoms  in  the  unit  cell.  The  accurate  calculation  of 
elastic  constants  is  an  extremely  involved  task  in  the  case  of  internal  displacements. 

A  calculation  could  be  made  in  the  approximation  which  neglects  internal  strain,  and 
this  was  begun.  Subsequent  considerations  indicated  that  such  a  procedure  is  regarded 
as  providing  only  an  order-of-magnitude  estimate.  Therefore,  elastic  constant 
calculations  have  not  been  pursued  further,  since  they  would  divert  the  effort  from  the 
main  task  as  related  to  heteroepitaxy. 

One  aspect  of  the  use  of  elastic  constant  data  should,  however,  be  discussed. 

In  interatomic  potential  calculations  for  simple  cubic  metals, as  found  in  the  literature, 
the  calculation  of  elastic  constants  and  the  requirement  that  the  strain  energy  be 
positive-definite  is  usually  imposed  for  lattice  stability.  In  such  cases  this  is  an 
important  constraint,  because  the  force  on  an  atom  vanishes  identically  due  to  lattice 
symmetry  and  therefore  imposes  no  constraint  at  all  in  determining  an  appropriate 
potential.  The  A^203  lattice  is,  perforce,  unstable  unless  the  potentials  are 
determined  such  that  the  required  forces  are  zero. 

As  mentioned  above,  potentials  for  Af203  have  been  determined  for  a  range 
of  values  of  the  parameters  Rac  and  Raa*  general  grounds  and  from  examination 
of  the  constraint  analytical  expressions,  it  appears  that  Raa  ^AC  would  lie 
in  the  range  of  1.969  A  -  Raa  Rac  -  2.524  A.  The  plan  is  to  carry  through 
subsequent  investigations  with  perhaps  two  different  potentials.  For  example,  two 
different  acceptable  potentials  can  be  distinguished  from  one  another  on  the  basis  of 
the  relative  amounts  of  anion-anion  and  anion-cation  energy  represented  in  the  energy  per 
unit  rhombohedral  cell.  At  present,  the  lattice  potentials  being  employed  for  surface 
reconstruction  mvestigations  to  b^  discussed  subsequently  are  for  Dac  “  0-  77352  eV, 

®AC  “  1.31819  A“l,  Rac  =  2.32  A  for  the  anion-cation  interaction  and  Daa  =  0.14484 
eV,  =  1.99225  A-l,  Raa  =  2.32  A  for  the  anion-anion  interaction. 

Early  work  in  heteroepitaxy,  rather  naturally  but  without  other  foundation,  made 
extensive  use  of  the  concept  of  minimum  mismatch  as  a  determining  factor  in  establish¬ 
ing  relative  film-substrate  orientations.  However,  low-energy  electron  diffraction 
(LEED)  studies  of  a  number  of  film  and/or  substrate  materials  have  already  established 
that  some  sort  of  surface  "reconstruction"  sufficient  to  yield  modified  LEED  patterns 
is  the  rule  rather  than  the  exception.  Thus,  the  surface  periodicity  and/or  symmetry 
can  be  altered  from  that  of  the  bulk  crystal. 
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An  investigation  of  surface  reconstruction  in  basal-plane  Ai203  has  been 
initiated,  using  the  potentials  just  discussed.  Although  basal-plane  ALO^  is  not  the 
substrate  orientation  of  greatest  practical  importance,  it  is  particularly  convenient  for 
beginning  studies  of  reconstruction.  The  probable  goal  of  these  studies  is  to  predict 
or  calculate  the  relative  energies  of  the  various  (1x1),  s/'3,  3  sTS,  and/^i  9  deg 

LEED  patterns  described  by  Chang  (Ref  7).  However,  it  is  not  clear  if  all  of  these 
patterns  can  fst  within  the  framework  of  the  present  model.  Also,  there  is  some 
evidence  that  impurities  and/or  surface  stoichiometry  may  be  determining  factors 
producing  the  various  LEED  patterns.  Thus,  additional  experimental  information 
would  be  helpful  in  establishing  which  atomic  surface  configuration  should  be  examined 
theoretically  and  in  identifying  the  necessary  theoretical  ingredients. 

Somewhat  more  basic  but  still  related  to  surface  stoichiometry  is  the  question 
of  the  composition  of  the  Ai203  substrate  at  Si  deposition  temperatures.  Early  workers, 
in  applying  mismatch  concepts  at  the  interface,  generally  considered  the  deposited  Si 
atoms  to  bond  onto  a  bare  plane  of  O  atoms  with  the  Si  atoms  occupying  the  missing 
A£  atom  sites.  On  the  other  hand,  since  each  O  plane  in  A£203  is  separated  by  two 
Ai  planes  with  long  and  short  bonds  to  the  O  planes,  Chang  (Ref  7)  suggests  that 
mathematical  cleavage  could  most  easily  take  place  between  two  A  planes,  leaving 
one  Ai  plane  at  the  surface  of  each  of  the  cleaved  half-crystals.  Intuitively,  it  appears 
that  there  would  be  considerable  rearrangement  of  Ai  surface  atoms  if  this  is  the  ease. 
Thus,  one  goal  of  the  present  work  will  be  to  examine  both  simple  alternatives  for  the 
surface  configuration  and  determine  what  conclusions  ean  be  reached. 

Studies  of  surface  reconstruction  of  basal- plane  Af203  have  been  started 
by  considering  the  relaxation  of  the  four  planes  nearest  to  the  surface.  Table  II 
lists  the  planes,  atom  type,  and  hexagonal  coordinates  of  the  15  atoms  which  can  be 
considered  to  represent  the  upper  half  of  a  30-atom  hexagonal  unit  cell.  For  this  half 
cell  the  origin  of  coordinates  is  at  an  Ai  atom  hole  (i.e. ,  the  site  of  a  missing  atom). 
The  upper  boundary  of  the  half  cell  then  also  contains  an  Ai  hole  and  defines  a 
hexagonal  network  of  holes. 

In  treating  surface  reconstruction  it  will  be  assumed  for  definiteness  that  the 
surface  is  an  O  plane,  labeled  plane  2  in  Table  II.  Plane  1  is  therefore  excluded  and 
planes  2,  3,  4,  and  3  are  included  in  the  relaxation  treatment. 

Next,  consider  a  mathematieally-cleaved  ideal  (0001)  Ai203  crystal.  Since  this 
half-crystal  still  possesses  three-fold  symmetry  normal  to  the  surface  and  along  the 
z  (e-axis)  direction  and  in  addition  still  retains  translational  periodicity  in  directions 
lying  in  the  surface,  it  is  clear  that  there  must  be  appreciable  relaxation  taking  place 
which  possesses  this  rotational  and  translational  symmetry.  There  will  also 
probably  be  relaxation  which  is  not  symmetric,  but  its  origin  is  somewhat  obscure  and 
may  relate  to  impurities  or  surface  stoichiometry. 

Consider  relaxations  of  atoms  in  planes  2-5  which  are  consistent  with  symmetry. 
Two  such  relaxations  are  possible.  The  first  type  is  associated  with  the  parameter  x 
in  Table  II,  which  defines  the  positions  of  C)  atoms  in  the  unit  cell  and  is  equal  to  0.306 
for  bulk  Ai20j,.  With  the  absence  of  translational  symmetry  in  the  a4  or  z  direction, 
an  allowable  change  would  be  for  each  O  plane  to  relax  toward  a  new  value  of  x  differ¬ 
ing  from  the  bulk  value.  Intuitively,  it  appears  that  the  (1  plane  at  the  surface  would 
show  the  largest  change.  The  magnitude  of  this  type  of  relaxation  has  not  yet  been 
numerically  determined.  Preliminary  indications  are,  however,  that  it  is  small. 
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Atom 

Type 

Hexagon&l 

Coordinates  a^,  a 9 

84 

Initial 

Coordinate 

34 

Relaxed 

Coordinate 

At  hole 

(0,0) 

0.5000 

Af 

(2/3,  1/3) 

0.481333 

- 

O 

(1/3-x,  2/3);  (1/3,  2/3-x);  (x+1/3,  x+2/3) 

0.41666 

0.43853 

Af 

(0,0) 

0.3520 

0.34913 

Af 

(1/3,  2/3) 

0.31466 

0.31629 

O 

(x,0);  (O.x);  (1-x,  1-x) 

0.250 

0.24300 

Af 

(2/3,  1/3) 

0.18533 

- 

Af 

(0,0) 

0.1480 

- 

0 

(2/3-x,  1/3);  (2/3,  l/3-x);  (2/3+x,  l/3+x) 

0.0833 

- 

Af 

(1/3,  2/3) 

0.01866 

- 

Af  hole  j 

(0,0) 

0.0000 

type  of  relaxation  is  that  in  which  each  of  the  various  planes 
considered  is  allowed  to  move  individually  in  the  54  or  z  direction.  This  type  of 
relaxation  appears  sigmficant  for  (0001)  AioOg,  an^  numerical  results  have  been  ob- 
tained.  The  rclaxcd-planc  positions  were  determined  by  an  iterative  procedure 
which  relaxed  planes  2,  3,  4  and  5  in  turn  toward  the  position  of  maximum  energy 
until  convergence  was  obtained.  The  results  obtained  are  nearly  as  expected  and  are 
shown  in  the  last  column  of  Table  II.  Plane  2,  the  surface  plane  of  O  atoms  In  this 
case  shows  the  greatest  (outward)  displacement.  Planes  3  and  4  are  displaced  only 
slightly  a^  tend  to  move  somewhat  closer  together.  Plane  5  is  seen  to  move  inward 
and  away  from  the  surface  a  moderate  amount.  This  negative  displacement  could 
have  been  expected  on  general  grounds  from  the  work  of  Tick  and  Witt  (Ref  8)  who 
found  a  similar  negative  displacement  for  Ni-like  metal  systems. 

Other  types  of  relaxation  are  to  be  investigated  in  the  near  future.  The  next 
case  already  being  examined  is  identical  to  that  reported  above  in  that  the  constraint 
of  symi^try  is  imposed.  However,  in  this  instance  the  surface  plane  will  be  an  Af 

configurations  determined  with  the  constraint  of  symmetry 
^11  then  be  employed  as  the  basis  for  more  complicated  types  of  surface  reconstruc¬ 
tion  in  which  symmetry  is  not  maintained.  In  this  way.  it  is  hoped  that  energetics  of 
various  types  of  substrate  surface  reconstruction  can  be  determined.  These 
reconstruction  may  play  a  major  role  in  achievement  of  the  goal  of  predicting  the  Si 
overgrowth  orientation  on  Af20g  substrates. 
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b.  Application  of  Electron-on- Network  Technique 

In  the  theoretical  studies  at  UCLA  the  possible  application  of  the  electron-on- 
network  technique  of  Montroll  ^  al  (Ref  9)  to  heteroepitaxial  sy:'tems  has  been  under 
investigation  (Ref  3).  The  procedure  for  adapting  the  technique  to  cases  where  lattice 
parameters  near  the  surface  are  allowed  to  vary  was  developed  and  applied  to  calcula¬ 
tions  for  the  metal  Na;*  the  plan  was  to  then  apply  the  procedure  to  the  case  of  AfnOo. 

Although  the  method  is  simple  in  principle,  it  became  clear  that  significant 
difficulties  would  be  encountered  in  applying  it  to  heteroepitaxial  systems.  For 
example,  in  the  Na  case  the  evaluation  of  electron  density  using  network  wave  functions 
was  found  extremely  difficult  in  view  of  the  transcendental  distribution  of  parameters 
indexing  the  wave  functions,  although  the  computation  was  completed  in  spite  of  this. 
However,  this  cast  considerable  doubt  on  the  feasibility  of  using  the  method  for  an 
accurate  assessment  of  interfacial  binding. 


Further  investigation  of  the  method  for  use  in  heteroepitaxial  systems  has  led 
to  the  conclusion  that  it  is  not  adequate  and  should  not  be  further  pursued  for  this 
purpose.  However,  the  use  of  the  technique  for  computing  variations  in  work  function 
and  surface  double-layer  potentials  for  single-crystal  monovalent  metals  has  led  to 
a  technical  paper  on  the  subject,  and  this  is  reproduced  in  Appendix  I. 

In  summary,  the  model  uses  an  electronic  wave  function  which  is  defined  to 
be  nonzero  only  along  the  lines  connecting  first  nearest  neighbors  in  the  metallic 
lattice.  The  electrons  are  summed  to  move  freely  along  the  lines  between  nearest 
neighbors.  No  electron-electron  or  electron-nucleus  force  is  included  in  the  model 
calculations  (except  for  forces  arising  from  the  Pauli  exclusion  principle). 

The  work  function  is  defined  as  the  amount  of  energy  required  to  move  an 
electron  from  a  point  slightly  inside  the  crystal  to  a  point  slightly  outside.  The 
contribution  of  the  electronic  double  layer  is  included  in  the  calculation  of  the  work 
function  as  well  as  the  dependence  of  the  double-layer  potential  on  the  surface  geometry. 
Surface  states,  where  the  electron  is  localized  in  the  neighborhood  of  the  face  of  the 
crystal,  are  found  to  have  energies  sufficiently  above  the  Fermi  level  to  eliminate  the 
possibility  that  they  make  any  contribution  to  the  double  layer  potential  for  the  case  of 
the  (100)  crystal  plane.  Consequently,  surface  states  were  ignored  in  all  the 
calculations. 


The  surface  double  layer  is  assumed  to  be  caused  by  the  presence  of  a  finite 
potential  barrier  at  the  surface  of  the  crystal.  Bulk  electronic  wave  functions  can 
penetrate  this  barrier  and  decay  exponentially  outside  the  crystal.  The  only  parameters 
required  by  the  model  are  the  nearest-neighbor  distance  for  the  lattice  and  the  height 
of  the  potential  bari’ier  at  the  surface.  The  former  quantity  is  fixed  by  the  lattice 
structure  (body  centered  cubic  for  the  alkali  metals)  and  by  the  density,  while  the  latter 
quantity  can  be  adjusted  to  give  the  best  agreement  between  the  model  calculations  and 
experiment. 


♦The  specific  calculations  for  Na  were  carried  out  at  UCLA  on  a  separately-funded 
study  not  related  to  this  contract. 
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For  the  alkali  metals,  Li  through  Na,  the  best  value  of  the  barrier  height  is 
about  50  percent  of  the  sum  of  the  ionization  potential  energy,  the  heat  of  vaporization, 
and  the  calculated  Fermi  level  for  the  corresponding  metal.  In  addition,  the  value  of 
the  double-layer  potential  found  for  Na  agrees  very  well  with  a  more  sophisticated 
calculation  in  the  literature  and  is  reasonably  close  to  the  experimental  measurement. 
Further  details  are  given  in  Appendix  I. 

c.  Calculations  of  Substrate-Induced  Anisotropy  in  Electrical  Properties  of  Si/Af20g 

It  has  long  been  recognized  that  differential  thermal  contraction  between  Si  film 
and  Ai203  substrate  will  produce  substantial  stress  In  the  Si  film.  Through  the  piezo¬ 
resistance  effect,  this  stress  acts  to  modify  the  electrical  properties  of  the  semi¬ 
conductor,  including  the  carrier  mobility. 

An  assumption  which  is  usually  implicit  in  most  semiconductor  thin-film  material 
evaluations  is  that  the  electrical  properties  are  isotropic  and  thus  do  not  depend  upon 
orientation  in  the  plane  of  the  film.  However,  within  this  reporting  period,  it  was 
determined  experimentally  that  most  Si/Ai203  films  exhibit  some  degree  of  mobility 
anisotropy  in  the  plane  of  the  film.  The  experimental  results  are  reported  in  detail 
under  Subtask  6  in  a  subsequent  section  of  the  report. 

In  order  to  explain  theoretically  the  anisotropy  in  carrier  mobility,  the  changes 
in  mobility  and  resistivity  due  to  the  stress  originating  in  the  differential  thermal 
contraction  of  film  and  substrate  have  been  calculated.  For  the  sake  of  completeness, 
the  difference  in  thermal  expansion  coefficients  parallel  and  perpendicular  to  the  c-axis 
in  AfgOg  has  also  been  taken  into  account.  Theoretical  formulae  have  been  developed 
which  are  valid  and  applicable  for  both  Si  and  Ge  films  and  for  AfgOg  and  spinel  sub¬ 
strates.  The  theory  has  been  formulated  generally,  and  (100),  (110),  (111)  and  (221)  Si 
film  orientations  have  been  treated  specifically  as  well  as  certain  "off-orientation"  cases. 

Some  of  the  numerical  results  are  presented  in  the  discussion  of  the  work  of 
Subtask  6.  A  comprehensive  theoretical  paper  will  be  submitted  for  publication  within 
the  next  few  months,  so  the  detailed  description  of  the  theory  will  not  be  reported  here 
at  this  time.  Some  of  the  results  obtained  from  the  theory  of  (100)-  and  (221)-oriented 
Si  films  are,  however,  contained  in  Ref  10  and  Ref  11,  respectively. 

WTiile  effects  of  stress  on  carrier  mobility  in  heteroepitaxial  films  have  been 
previously  recognized  (e.g. ,  Ref  12),  the  investigations  reported  here  are  the  first 
which  have  explicitly  recognized  and  treated  the  orientational  anisotropy  in  mobility 
which  arises  from  differential  thermal  contraction  and  the  relative  orientation  of  film 
and  substrate. 


This  study  of  mobility  anisotropy  is  potentially  veiy  significant  in  that  it  promises 
a  more  detailed  understanding  of  mobility  in  heteroepitaxial  semiconductor  films  than 
it  has  been  possible  to  obtain  from  previous  investigations  of  the  mobility. 
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2.  SUBTASK  2.  DEPOSITION  STUDIES  AND  FILM  PREPARATION 

During  the  first  year  of  the  program  the  emphasis  on  this  subtask  was  placed  on 
determining  the  effect  of  experimental  growth  parameters  on  the  quality  of  Si  epitaxial 
films  grown  by  the  CVD  method  of  pyrolysis  of  SiH4  o”  substrates  of  various  orienta¬ 
tions  of  A#20g  and  MgAf204.  It  was  established  for  the  growth  system  used  that 
autodoping  occurs  in  Si  on  A#203  at  temperatures  greater  than  about  1050  C,  so  a 
concerted  study  was  made  which  considered  the  effects  of  such  factors  as  growth 
temperature,  growth  rate,  and  nucleation  phenomena  at  or  below  this  temperature 
(Ref  2). 

Much  of  this  work  was  rendered  invalid  by  a  very  important  effect  which  had 
not  previously  been  reported  in  the  technical  literature  dealing  with  the  properties  of 
heteroepitaxial  semiconductor  films.  Specifically,  it  was  found  that  the  electrical 
properties  of  Si  films  on  insulators  with  carrier  concentrations  <10^6  cm“3  are 
dominated  by  surface-state  conduction.  Consequently,  further  studies  were  made  on 
n-type  films  doped  to  cm“3. 

Simultaneous-growth  studies  on  substrates  which  produced  (100)  Si  and  (111)  Si 
below  the  autodoping  temperature  revealed  the  electrical  properties  of  the  films  to 
be  very  dependent  on  substrate  orientation.  Although  (111)  Si  growth  on  (10l4)  AfgOg 
was  found  to  be  essentially  equivalent  electrically  to  (100)  Si  growth  on  (0112)  Af203'^ 
below  the  autodoping  temperature  (~1050  C),  changing  the  substrate  orientation  to 
near  the  (1120)  A(20^  plane  (not  previously  used  in  epitaxy  studies)  increased  the 
mobility  of  (111)  Si  films  by  about  20  percent  for  2  fim-thick  films  and  even  more 
(~50  percent)  for  5  fxm-thick  films,  for  the  10^®-10l7  cm-3  carrier  concentration  range 
(Ref  2).  It  was  suggested  at  the  time  that  the  better  match  of  the  thermal  expansion 
coefficients  for  this  Af203  orientation  and  for  Si  may  have  been  responsible  for  the 
improvement. 

Analysis  of  the  data  obtained  during  the  third  six-month  period  (Ref  3)  further 
emphasized  the  strong  interrelationships  existing  among  the  various  parameters  selected 
as  most  meaningful  for  optimizing  Si  growth  on  insulators.  Evaluation  of  the  electrical 
properties  of  Si  films  on  Af203  orientations  identified  several  years  ago  as  providing 
reflective  Si  overgrowths  demonstrated  that  growth  conditions  (1)  must  be  optimized 
for  each  substrate  orientation  chosen;  (2)  differ  for  those  Af203  orientations  which 
lead  to  the  same  Si  orientation;  (3)  are  dependent  upon  reactor  geometry  and  gaseous 
atmosphere;  and  (4)  should  be  optimized  for  the  particular  film  thickness  desired. 

The  effects  of  annealing  during  film  growth  and  the  use  of  water-cooled  reactors 
were  investigated.  Preliminary  studies  of  Si  growth  by  SiH4  pyrolysis  at  reduced 
pressures  (1  to  10  torr)  were  carried  out  and  were  very  encouraging,  indicating  single¬ 
crystal  growth  can  be  obtained  over  a  fairly  wide  temperature  range,  when  conditions 
are  optimized,  on  both  Af203  and  MgAf204  substrates.  Brief  investigation  of  the 
growth  of  Si  films  on  both  A^203  ®nd  MgAf204  using  He  as  the  growth  atmosphere  and 
as  the  carrier  gas  showed  that  epitaxial  growth  could  be  achieved,  but  the  conditions 
for  good  quality  growth  were  not  established. 


14 


The  work  on  this  subtask  during  the  fourth  six-month  period,  described  herein, 
was  conc^entrated  on  continued  attempts  to  optimize  the  Si  deposition  process  for  growth 
on  (1120)-oriented  and  (0112)-oriented  A#20g  surfaces.  These  efforts  involved  (1) 
examination  of  the  effects  of  post-nucleation  annealing  on  ultimate  film  properties; 

(2)  evaluation  of  the  effects  on  film  growth  of  gas-phase  etching  of  the  AfoO,  surface 
prior  to  deposition;  (3)  investigation  of  the  effect  on  the  growth  process  of  using  a 
water-cooled  reactor  chamber;  and  (4)  preparation  of  a  variety  of  Si/AfoOo  film 
composites  for  the  study  of  the  anisotropy  of  electrical  properties  in  this  system  (see 
Subtask  6).  In  addition,  an  extensive  study  of  the  growth  and  properties  of  Si  films 
on  Czochralski-grown  stoichiometric  spinel  (MgAf204)  substrates  was  begun. 

It  should  be  noted  that  film  properties  and  CVD  growth  conditions  described  in 
the  preceding  semiannual  report  (Ref  3)  were  indicative  of  the  achievement  of  near- 
optimum  conditions  for  Si  growth  on  several  Af203  orientations,  using  the  vertical 
reactor  system.  There  was  considerable  spread  in  the  results  obtained  to  that  date, 
but  much  of  that  can  now  be  explained  based  on  the  discovery  and  subsequent  investiga¬ 
tion  of  strong  anisotropy  effects  in  the  Si/A#203  system.  (See  Subtasks  1  and  6.) 

a.  Effect  on  Film  Properties  of  Annealing  during  Growth 

To  examine  further  the  effects  of  annealing  Si  films  at  early  stages  of  growth  a 
group  of  Si/A(203  samples  was  prepared  with  the  deposition  process  interrupted  at 
various  stages  to  permit  a  20-minute  annealing  pieriod  followed  by  resumption  of  growth 
to  the  desired  final  thickness.  Electrical  measurements  were  then  made  on  the  Si  films 
to  evaluate  the  results  of  the  annealing  process.  Preliminary  studies,  described  in 
earlier  monthly  reports,  had  indicated  this  treatment  was  not  detrimental  to  a  2000  A- 
thlck  film,  but  its  effect  on  the  ~1000  A-thick  films  grown  in  tha  more  recent  experi¬ 
ments  was  not  known. 

The  data  obtained  are  presented  in  Table  III,  in  which  it  can  be  seen  that  a 
0.67  Km-thick  (111)  Si  film  annealed  at  800  C  exhibited  a  carrier  mobility  of  490  cm2/ 
V-sec,  considerably  better  than  the  mobility  of  190  cm2/V-sec  measured  for  another 
;,nn  annealed  at  1100  C.  The  data,  although  sparse,  seem  to  indicate  that 

1100  C  is  better  than  1050  C  as  the  growth  temperature  for  (111)  Si  on  ~(1120)  A<’203. 

^^®  effects  of  annealing  at  intermediate  stages  during  the  growth 
of  Si  films  on  AfoOg  involved  the  growth  and  anneal  of  ~2000  A  of  Si  at  a  growth  rate 
of  ~4  fim/min  followed  by  further  growth  at  the  same  rate.  The  results,  shown  in 
Table  IV,  confirm  that  the  annealing  of  veiy  thin  Si  films  on  Af20,  followed  by  additional 
filrn  growth  does  not  impair  the  ultimate  mobilities  of  films  so  grown.  The  best  films 
in  this  series  of  depositions,  annealed  or  unannealed,  were  grown  on  ~(1120)-oriented 
At203  substrates;  mobilities  >400  cm2/v-sec  w'ere  obtained  in  films  only  0.6-0. 7  fim 

IfllCK* 


However,  the  annealing  studies  to  date  have  not  demonstrated  a  significant 
improvement  in  the  overall  quality  of  Si  films  grown  by  this  deposit-anneal-deposit 
technique  realtive  to  those  prepared  by  the  conventional  single  -deposition  method  On 
this  basis,  it  was  determined  that  no  further  studies  of  this  technique  would  be  pursued 
during  the  second-year  program. 
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TABLE  III. 

ElectricaJ.  Properties  of  Si  Films  Annealed  for  20  min  after  Growth  of 
~1000  A  Thicknsaii  cn  Ai203  Substrates  (Growth  Rate  ~2tJ-m/min) 
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ibers  represent  single  determination  of  electrical  parameter  by  standard  techniques. 


Use  of  these  substrates  for  Si  growth  clearly  demonstrated  that,  just  as  for 
mechanically  polished  Af203,  the  appearance  of  the  surface  can  be  very  misleading. 
Some  surfaces  that  look  scratch-free  really  contain  scratches  that  are  filled  in  by 
polishing  powder  and  debris.  This  condition  was  revealed  during  attempts  to  clean 
MgA£204  surfaces  for  Si  growth;  it  was  found  that  after  about  10  to  30  min  in  hot  HF, 
most  of  the  "fill-in"  material  is  removed.  Figure  3  shows  adjacent  regioPiS  of  the 
surface  of  a  broken  slice  of  (111)  MgAf204;  the  right  side  was  "etched"  for  30  min  in 
hot  HF.  The  defects  showai  are  apparent  in  MgAf204  polished  by  either  of  the  two 
principal  vendors*  as  well  as  by  NR  laboratories. 

Just  as  was  the  case  for  A^203  substrate  surfaces,  it  was  not  evident  from  the 
quality  of  Si  films  grown  on  these  polished  but  defective  substrates  that  the  suifaces 
were  as  poor  as  later  scrutiny  showed  them  to  be.  The  explanation  for  the  high 
quality  of  Si  films  that  can  be  growm  on  relatively  poor  substrate  suifaces,  both  for 
Af203  and  for  MgAf204.  is  yet  to  be  found.  It  seems  clear  that  chemical  or  combined 
chemical-mechanical  polishing  processes,  probably  used  in  conjunction  with  various 
annealing  procedures,  should  be  studied  further  as  means  for  removing  or  minimizing 
surface  w’ork  damage  on  IVIgAf204  substrates  intended  for  epitaxy  use. 

It  is  important  to  note  that,  in  the  studies  on  this  program  to  date,  higher 
quality  Si  films  have  been  grown  on  MgAfgO^  substrates  than  on  AI202  substrates. 
Numerous  variables  remain,  however,  and  this  is  not  necessarily  regarded  as  a  final 
conclusion. 


a  b 

Figure  3.  (Ill)  Czochralski  MgAf204  (a)  before  Etching; 
(b)  after  Etching  in  Hot  HF  for  30  min 


*Union  Carbide  Corp.  and  Crystal  Technology,  Inc. 
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Electrical  Properties  of  Si  Filins  Annealed  for  20  min  after  Growth  of 
-2000  A  Thickness  on  A/2O3  Substrates  (Growth  Rate  ~4^m/min) 
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**No  intermediate  annealing  step  used. 


b.  Effect  of  Gas-Phase  Etching  of  Substrates  on  Film  Properties 

Gas-phase  etching  has  revealed  the  presence  of  considerable  surface  and 
subsurface  work  damage  on  A#202  that  is  caused  by  the  mechanical  processing  of 
the  substrate  in  preparing  a  surface  conducive  to  epitaxy.  Yet  it  is  striking  that  such 
damaged  surfaces  are  being  used  in  the  industry  for  the  epitaxial  growth  of  Si  and 
other  semiconductors. 

Experiments  were  carried  out  with  the  various  substrate  surfaces  shown  in 
Figures  6  and  7  of  the  Third  Semiannual  Report.  (Ref  3):  Si  films  were  grown  and 
their  electrical  properties  were  measured.  The  data  are  shown  in  Table  V  and  are 
somewhat  surprising.  Essentially  identical  carrier  mobilities  were  found  for  films 
grown  on  unetched  substrates  and  for  those  etched  for  as  long  as  100  min  ( -20  fim  of 
thickness  removed),  for  both  Czochralski-and  Verneuil-grown  Ai203.  This  was  taken 
to  indicate  that  (1)  the  work  damage  is  even  greater  than  20  |im  deep  in  these  sub¬ 
strates  (which  were  prepared  by  conimercial  vendors),  or  (2)  that  the  Si  overgrowth 
is  able  to  "overcome^'  much  of  the  substrate  surface  imperfection  during  its  growth, 
with  the  actual  mobility-limiting  factor  being  the  stress-producing  thermal  contraction 
and/or  lattice  mismatch  between  the  two  materials,  or  (3)  that  there  is  a  third  species 
present  on  the  substrate  surface  in  this  CVD  reaction  ai^  that  it  provides  the  real 
growth-influencing  sites. 

Further  experiments  are  needed  to  examine  this  important  question  in  more 
detail. 

c.  Effect  of  Use  of  Water- Jacketed  Reactors 

Si  growth  experiments  were  performed  in  a  water-jacketed  reactor  chamber, 
to  compare  the  electrical  properties  of  films  grown  in  a  water-cooled  reactor  with 
those  of  films  grown  without  the  reactor  walls  water-cooled.  The  results  are  given 
in  Table  VI  for  films  grown  both  with  and  without  water  cooling.  These  results  show 
that  there  is  no  apparent  advantage  in  the  use  of  a  water-jacketed  reactor  in  the 
vertical  reactor  system  used  in  these  experiments;  in  fact,  the  additional  handling 
necessary  to  detach  and  attach  the  reactor  from  and  to  the  system  was  found  to  be 
a  distinct  disadvantage. 

At  low  growth  rates  the  amount  of  reactor  wall  deposit  appeared  to  be  less  in 
the  case  of  the  cold— wall  reactor  than  for  the  uncooled  reactor.  However,  because 
of  the  additional  handling  required  and  in  view  of  the  fact  that  the  electrical  proper¬ 
ties  were  essentially  the  same  for  the  two  cases, no  further  studies  are  presently 
contemplated  with  the  cold-wall  reactor. 


d.  CVD  Growth  of  Si  on  Spinel 

Studies  of  the  growth  of  Si  on  spinel  were  continued  during  this  six-month 
period,  with  "G-D’s"  grade  He  as  the  major  carrier  gas.  Films  have  been  grown 
primarily  on  polished  w’afers  of  Czochralski  spinel* 


*Spinel  boules  and  pc  .ished  wafers  obtained  from  Union  Carbide  Corporation  and 
from  Crystal  Technology,  Inc. 
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TABLE  VI. 

Properties  of  Si  Films  Grown  in  Jacketed  Reactor 


Film 

Orientation 

Film  Thickness 
(fxm) 

Resistivity 
(ohm -cm) 

Carrier  Concentration  Mobility 

(cm“'l)  (cm2/V'-sec) 

With 

out  Water  Cooling 

(100) 

1.7 

0.71 

1.9x10^*’  400 

(111) 

1.8 

0.78 

1.3  xlO^*’  030 

Wi 

th  W'ater  Cool 

ing 

(100) 

2.0 

0.20 

5.0x10**’  490 

(111) 

1.9 

0..‘15 

2.9x10**’  020 

in  these  studies,  Pd-piirified  II9  was  added  to  the  He  carrier  gas  in  order  to 
reduce  the  wall  deposits  and  improvelhe  reaction  efficiency.  At  growth  temperatures 
of  1025  C.  for  example,  Si  growth  rates  in  He  (even  with  the  II2  addition)  were  reduced 
consideral)ly  with  respect  to  those  obtained  with  only  II2  as  the  carrier  gas.  For 
example,  a  Si  114  flowmeter  setting  of  .5  (25  ccpm  pure  Sill  t)  resulted  in  a  Si  growth  rate 
of  ~2  fim/min  in  Ho  I)Ut  only  0.4  fim/min  in  He  contatning  small  amounts  of  added  II2. 

The  effect  of  H2  additions  to  He  at  high  Silli  concentrations  is  shown  in  Figure  1. 
The  Sin.j  flow  was  kept  constant  at  150  ccpm  and  Ihe  Si  gi’owth  len’iperatiii’e  was 
1025  C.  The  data  show  that  es.scntially  a  linear  increase  in  growth  rate'  results  when 
112  is  added  to  the  He  carrier  f^as,  at  least  up  to  the  indicated  flow  rate  of  2000  ccpm 
(.‘iH  percent  of  the  total  carrier  gas).  Since  the  growth  rate  in  pure  II2  has  been  found 
to  be  ~15pm/min,  the  dependence  on  112  percentage  is  not  linear  throughout  the  range, 
leveling  off  at  higher  concentrations  of  ri2. 

It  was  also  estaUished  that  for  a  given  IIe-Il2  ratio  the  growth  rate  is  veiy 
temperature-dependent.  This  is  shown  in  Figure  2  for  a  constant  SilL  flow  (nominally 
150  ccpm)  and  0.8  percent  II2  in  lie.  Over  the  lOO-degree  temperature  span 
investigated  the  growth  rate  fs  seen  to  change  by  a  factor  of  4.  being  lower  at  the 
higher  temperature. 

It  is  difficult  to  understand  how'  higher  temperatures  cause  such  a  large 
decrease  in  the  growth  rate.  The  liehavior  is  very  different  from  that  obsen'ed  in 
a  pure  II2  carrier  gas  atmosphere;  distinct  mechanisms  appear  to  be  operative 
in  the  two  cases.  The  phenomenon  may  be  similar  to  that  observed  in  preliminan- 
experiments  with  growth  of  GaAs  on  GaAs  substrates  by  the  pyrolysis  of  trimetlnl- 
gallium  and  arsine_(AsIl3);  in  that  case  the  growth  rate  is  greater  on  the  (111)  "A" 
face  than  on  the  (III)  "B"  face  of  GaAs,  a  difference  presumably  caused  by  a  difference 
in  the  charge  on  the  surfaces  of  the  i  111)  faces. 

Since  it  w’as  earlier  found  that  the  net  carrier  concentration  in  As-dof)ed  Si 
films  is  growth-rate  dependent  when  II2  is  the  carrier  gas  (Figure  19.  Ref, 5),  it  w-as 
expected  that  similar  effects  would  be  obsened  in  the  Ile-I^  system.  The  results 
obtained  in  an  examination  of  this  question  are  shown  in  Table  VII.  The  data  indicate 
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Figure  1.  Changes  in  Si  Groulh  Rate  at  1025  C  Caused  by  II, 
Additions  to  He  Carrier  Gas  (Sill4  Flow  Rate  150  ccpm) 
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Figure  2.  Effect  of  Growth  Temperature  on  Si  Growth  Rate  for  SilU  Flow  Rate  of 

150  ccpm  and  Hg/He  -  0.  8  Percent 


TABLE  VII. 

Effect  of  Growth  Rate  on  Doping  Level  for  Si  Films  Grown  on  (111)  Spinel  at 
1025  C  in  Ile-IIg  Atmospheres  (Constant  Dopant  Gas  Flow  Rate) 


H„  in  He 

{%) 

Film 

Thickness 

(pm) 

Growth 

Rate 

(pm /min) 

Resistivity 
(ohm -cm) 

Carrier 

Concentration 

(cm“3) 

Mobility 

(cm^/V-sec) 

0.8 

1.90 

2.H 

0.33 

2.6  X  10^^’ 

747 

4.2 

1.87 

3.2 

0.41 

1.9  X  10^® 

792 

8.4 

2.05 

4.1 

0.46 

1.5  X  10^^ 

900 

16.8 

1.85 

5.8 

0.56 

1.6  X  10^^’ 

718 

33.3 

2.10 

9.7 

0.23 

3.0  X  10^® 

925 
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that  the  doping  level  does  not  change  appreciably  with  growth  rate  in  IIC-II2  atmos¬ 
pheres  at  1025  C.  The  dependence  of  doping  level  on  growth  temperature  has  not  yet 
been  determined. 

The  question  of  the  extent  of  autodoping  in  the  Si/spinel  system  has  been 
examined  in  preliminary  experiments.  Undoped  Si  films  were  grown  on  spinel  at  various 
deposition  temperatures  (1000  to  1100  C)  using  a  constant  amount  of  Sill^  and  H2-He 
mixture. 

As  shown  in  Table  VIII,  films  grown  in  the  1000  to  1050  C  range  were  n-type; 
at  1075  C  and  1100  C  they  were  p-type  before  annealing.  Annealing  in  O2  and  N2 
caused  the  films  to  become  more  strongly  n-type,  suggesting  the  removal  of  acceptor 
impurities  and/or  defects  by  the  annealing  process.  It  thus  appears  that  the  auto¬ 
doping  phenomenon  is  operative  for  Si/spinel  in  the  ne-H2  atmosphere  at  approximately 
the  same  temperatures  as  for  the  Si/Af203  system  (Ref  3). 


TABLE  VIII. 

Autodoping  in  Si/(111)  Spinel  Using  0.8  Percent  H2-in-He  as  Carrier  Gas 


Film 

Growih 

Film  Resistivity’ 

and  Type+* 

Deposition 
Temp  (C) 

Thickness 

(Km) 

Rate 

(fim/min) 

Before  Anneal 
(ohm -cm) 

After  Anneal 
(ohm-cm) 

1000 

2.1 

3.8 

4.2  X  10^  (n) 

1.3  (n) 

1025 

2.0 

3.1 

25.4  (n) 

0.4  (n) 

1050 

1.8 

2.1 

3.2  (n) 

0.28  (n) 

1075 

1.8 

1.2 

4. 8  X  10^  (p) 

76.  (n) 

1100 

4.4 

0.9 

1  X  10^  (p) 

1.  (n) 

(Reflective 
edge  reg  )n) 
6.6  X  105i  (p) 
(Gray  central 
region) 

Four-point -probe  measurement 
+  +  Thermoelectric  probe  determination 


Because  of  the  importance  of  the  phenomenon  it  is  planned  that  further 
investigation  of  autodoping  in  the  Si/spinel  system,  as  a  function  of  both  H2-He  ratio 
and  temperature,  will  be  carried  out  during  the  next  six  months. 
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3.  SUBTASK  3.  ANALYSIS  AND  PURIFICATION  OF  CVD  REACTANTS 


During  the  first  year,  techniques  of  gas  chromatography  were  developed  for 
analysis  of  the  reactants  used  for  semiconductor  heteroepitaxy  by  CVD.  A  general- 
purpose  gas-handling  system  was  constructed  for  the  highly  volatile  and  reactive 
gases  studied,  with  silicone  oil  and  polymer  colurtins  used  for  the  chromatography. 
Several  extraneous  impurity  peaks  were  observed  in  the  chromatograms  of  SiH4 
samples;  diborane  (B2ll(;)  was  tentatively  identified  as  a  significant  impurity  (-10  ppm), 
although  not  confirmed  by  mass  spectrometer  techniques.  Small  quantities  of  purified 
Sill^,  free  of  diborane,  were  prepared  by  successive  injections  in  the  chromatograph; 
quantities  were  too  small,  however,  fur  use  in  laboratory  CVD  experiments.  A  summary 
of  this  work,  which  was  carried  out  at  UCLA,  was  given  in  Appendix  A  of  the  Third 
Semiannual  Report  (Ref  3). 

In  the  first  six  months  of  the  second  year  the  emphasis  of  the  work  of  this  subtask 
was  changed.  Analyses  of  the  reactants  used  in  the  CVD  experiments  were  continued, 
primarily  by  sensitive  mass  spectrometric  techniques  with  the  assistance  of  external 
analytical  service  laboratories.  Samples  of  Sill4  and  of  trimethylgallium  (TMG)  used 
for  Si  and  GaAs  CVD  experiments  were  analyzed  for  impurity  content.  Disilane  and 
trimethylsilane,  together  with  several  other  impurities  of  less  concern,  were  found  in 
various  Sill4  samples.  The  analyses  of  TMG  samples  were  subject  to  various  uncer¬ 
tainties,  some  of  which  still  remain  unresolved. 

During  the  past  six  months  additional  reactant  analyses  have  been  carried  out, 
and  a  significant  new  study  of  the  chemistry  and  reaction  kinetics  of  the  CVD  processes 
used  for  depositing  Si  and  GaAs  films  in  heteroepitaxial  systems  has  been  initiated  in 
the  Department  of  Chemistry  at  California  State  University,  San  Diego.  These  two 
aspects  of  the  subtask  will  be  described  in  sequence. 

a.  Analysis  of  CVD  Reactants  for  Impurity  Content 

Four  different  tanks  of  "pure"  Sill4  have  been  used  in  the  studies  of  Si  deposition 
on  spinel  (MgAf204)  during  the  past  six  months.  Three  of  these  were  analyzed  for 
impurity  content  before  use  in  these  experiments  by  an  analytical  service  laboratory* 
using  sensitive  mass  spectrometric  techniques.  One  of  the  four  tanks  was  also 
reanalyzed  after  the  pressure  had  dropped  in  use  to  33  psi  from  an  original  value  of 
85  psi. 

These  data  are  given  in  Table  IX,  together  with  an  independent  analysis 
performed  by  another  laboratory '■+  and  supplied  with  one  of  the  SilL  tanks  by  the 
vendor.  First,  the  discrepancies  in  the  data  of  the  two  analyses  obtained  for  the 
tank  of  SilL  supplied  by  Vendor  B  are  significant.  The  analytical  procedures  used 
in  the  supplied  anal)'sis  apparently  are  not  able  to  detect,  for  example,  less  than 
3000  ppm  of  O2,  an  impurity  of  importance  in  SiIJ4  which  could  lead  to  the  formation 
of  Si02  (and  resulting  pinholes)  in  deposited  Si  films.  The  other  differ  i.Tces  in  the 
analytical  data  for  this  tank  indicate  insufficient  sensitivity  in  the  mass  spectrometer 
methods  used  in  the  analysis  to  properly  establish  the  concentrations  of  those 
impurities  which  could  be  detrimental  to  high-purity  Si  growth. 


*West  Coast  Technical  Service,  Inc.  (WCTS),  San  Gabriel,  CA. 
''♦Gollob  Analytical  Service,  Berkeley  Heights,  NJ 
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TABLE  DC. 

Mass  Spectrometric  Analyses!  of  SiH4  from  Three  Vendors 
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Note  also  the  relatively  high  conceit  rations  of  H2  and  disilane  in  the  tanks 
obtained  from  Vendors  A  and  B  compared  with  those  in  the  SiH4  prepared  by  Vendor  C. 
The  effects  of  these  impurities  on  Si  growth  seem  to  be  minimal  in  a  growth  system 
using  H2  as  the  carrier  gas.  However,  if  He  is  the  carrier  it  is  difficult  to  predict 
how  these  impurities  might  affect  the  nucleation  mechanism.  Another  source  of 
concern  is  the  CO2  content  in  the  tank  from  Vendor  B.  This  impurity  also  could  lead 
to  Si02  formation  in  a  H2  atm  ensphere  by  a  ’’water-gas"  reaction, 

C02  +  H2^C0  4.H20.  (1) 

the  water  then  reacting  with  SiH^  to  form  Si02.  The  impurity  levels  indicated  in  the 

table  preclude  further  use  of  SiH4  from  this  vendor  until  significant  improvement  in 
purity  is  achieved. 

The  remeasurement  of  the  impurity  levels  in  the  SiH4  from  Vendor  C  was 
prompted  by  a  sudden  inability  to  control  Si  film  growth  rates  within  reasonable 
limits  (~10  percent)  with  either  He  or  Hg  carrier  gases.  The  second  analysis, 
however,  did  not  provide  an  indication  of  the  cause  of  this  difficulty.  A  change  in 
SiH^  tanks  apparently  solved  the  problem  of  growth- rate  control,  although  the 
electrical  properties  of  the  films  grown  on  MgAf204  with  the  first  tank  of  Sin4 
were  quite  good;  mobilities  w-ere  between  700  and  925  cm2/V-sec  for  films  with 
carrier  concentrations  of  10li»-10l7/cm3, 

The  change  to  the  other  tank  of  SiH4  (supplied  by  a  fourth  vendor*  for  evaluation) 
was  coincident  with  the  growth  of  films  with  reduced  mobilities  -  lower  than  700  cm2/ 
V-sec,  With  the  change  of  SiH4  tanks  growth  and  dopant  control  had  to  be  re¬ 
established;  Si  films  w-ere  grown  on  (111)  MgAf204  with  mobilities  similar  to  those 
obtained  on  (0112)  AfgO^  (~400-500  cm2/V-sec)  but  inferior  to  those  obtained 
previously  on  (111)  lvrgAf204  (700-925  cm2/V-sec).  To  explore  the  situation  further 
the  condeasable  component  of  the  reactant  was  trapped  in  the  Si  114  line  at  -78  C  and 
transferred  to  a  specially  prepared  sampling  cylinder  for  mass  spectrometer 
analysis.  The  presence  of  H2O  and  tetrahydrofuran  (THF)  in  this  sample  was 
established  and  was  consistent  with  the  fact  (determined  subsequently)  that  the 
manufacturer  used  THF  as  a  solvent  in  the  preparation  of  the  SiH4.  The  low 
mobilities,  however,  were  not  consistent  with  the  fact  that  a  100  pm-thick  film  of 
Si  grown  on  p-type  Si  had  reasonable  electrical  properties  (P~  24.8  ohm-cm, 
n  =  1.98  X  10l-l/cm3,  and  P-  1270  cm2/v-sec).  It  appears  that  something  in  the 
S1H4  is  affecting  the  growth  quality  on  insulating  substrates.  The  results  also 
suggest  that  Sin4  and  HgO  can  coexist  in  a  closed  container  but  that  their  reaction  is 
catalyzed  by  exposure  to  air. 

This  impurity  problem  has  prevented  an  accurate  evaluation  of  the  effects  of 
deposition  parameters  on  Si  film  properties.  It  will  probably  be  necessary  to 
perform  such  studies  with  sufficient  quantities  of  starting  materials  to  complete  the 
entire  investigation,  with  the  hope  that  the  gas  purity  is  not  affected  by  changes  in 
total  tank  pressure.  In  view  of  this,  a  200-gm  quantity  of  Sill4  has  been  obtained 
from  Synthatron  Corporation  for  further  studies  of  Si  growth  on  MgAf204.  It  has  not 
yet  been  evaluated  for  impurity  content,  however. 


*  Liquid  Carbonic  Corporation 
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A  tank  of  AsHs-in-He  (nominally  10  percent  ASH3)  was  obtained  from  Airco 
for  GaAs  growth  studies.  In  evaluating  it  by  means  of  GaAs  growth  on  Af203,  using 
a  good  source  of  trimethylgallium ,  it  was  found  that  considerable  condensate  was 
obtained  at  -45  C.  The  tank  was  returned  to  the  vendor,  and  it  was  analyzed  by  West 
Coast  Technical  Service,  Inc.  ITie  gas  was  found  to  contain  the  following: 


ASH3 

12.2  mole  percent 

^«2 

4  ppm 

Water 

15  ppm 

H2S 

4  ppm 

^2 

56  ppm 

He 

balance 

Oxygen 

15  ppm 

Argon 

17  ppm 

The  high  concentration  of  water,  O2,  CO2,  and  H2S  could  be  expected  to  affect  the 
quality  of  the  GaAs  films.  A  replacement  tank  is  now  being  prepared  by  the  vendor. 

Another  tank  of  10  p'^rcent  AsH3-in-H2,  obtained  as  a  research  sample  from 
Phoenix  Research  Corporation  for  similar  evaluation  by  GaAs  deposition,  was 
described  in  the  vendor's  analytical  report  as  containing  no  impurity  greater  than 
1  ppm.  No  condensate  was  obtained  at  -45  C,  and  thick  undoped  GaAs  films  were  found 
to  be  p-type;  this  indicates  a  low  donor  impurity  content  in  the  reactant.  Several  GaAs 
films  were  grown  at  various  ASH3-TMG  ratios,  with  H2Se  added  as  dopant.  However, 
the  electrical  data  obtained  for  films  grown  on  Af203  using  this  AsH3-in-H2  were 
inconclusive;  it  was  found  that  the  dopant  levels  in  the  Se-doped  films  were  not  in  the 
10l5_io  cm“3  range,  desirable  for  permitting  mobilities  which  are  more  meaning¬ 
ful  for  evaluation  of  the  GaAs  film  quality.  Additional  GaAs  films  on  Af203  have  been 
grown  for  evaluation  of  this  particular  tank  of  AsH3-in-H2,  but  the  results  are  not  yet 
complete. 

b.  Study  of  Chemistry  and  Reaction  Kinetics  of  CVD  Processes 

During  the  past  six  months  arrangements  have  been  made  for  fundamental 
studies  of  the  chemistry  and  reaction  kinetics  of  the  CVD  processes  used  for  growing 
Si  and  GaAs  films  in  heteroepitaxial  systems  to  be  carried  out  by  scientists  in  the 
Department  of  Chemistry  at  California  State  University,  San  Diego  (CSUSD).  This 
group*  is  experienced  in  studies  of  chemical  reactions  of  the  Si  hydrides  and  related 
systems  and  is  in  a  position,  in  terms  of  present  work  and  past  experience,  to  make 
important  contributions  to  the  fundamental  studies  of  this  contract  program  and  to 
thereby  permit  further  improvements  in  the  manner  of  application  of  these  chemical 
reactions  to  the  practical  growth  of  semiconductor  thin  films  for  device  use.  For 
example,  the  long-standing  question  as  to  whether  the  release  of  Si  by  the  pyrolysis 
of  Sin4  occurs  in  the  gas  phase  near  a  substrate  surface  or  on  the  substrate  surface 
itself  can  probably  be  answered  by  the  experiments  to  be  undertaken  at  CSUSD.  This 
point  is  of  major  importance  to  the  planned  in  situ  CVD  experiments  in  the  electron 
microscope  (Subtask  5). 


*The  studies  will  be  carried  out  by  Prof.  M.  A.  Ring  and  associates. 
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In  these  studies,  an  investigation  will  be  initiated  into  the  basic  chemistry 
involved  in  the  formation  of  GaAs  from  the  metalorganic-hydride  reaction  of  tri- 
methylgallium  (Ga(CH3)3)  and  Asn3  and  in  the  deposition  of  Si  by  the  pyrolysis  of 
Si  14. 


The  flow  pyrolysis  of  trimethylgallium  (TMG)  with  ASH3  will  be  examined  over 
an  inert  surface  (Si  mirror)  to  determine  the  requirements  for  the  formation  of  pure 
GaAs.  The  parameters  to  be  varied  will  be  the  pyrolysis  temperature,  the  TMG- 
AsIIg  ratio,  and  the  partial  pressure  of  Hg.  In  addition,  the  volatile  products  will 
be  determined.  The  decomposition  modes  of  ASH3  and  TMG  will  be  examined  to 
determine  what  radicals  are  produced  in  the  pyrolyses  of  Asllg  and  of  TMG  and  in 
theco-pyrolysis.  The  reaction  of  these  radicals  with  TMG  ancf  AsHo  will  be 
investigated. 

The  pyrolysis  of  Sill^  and  Si2H6  over  Ai203  will  be  investigated.  The  AI2O3 
will  be  preheated  to  various  temperatures  (400  to  1100  C)  prior  to  each  pyrolysis 
to  determine  effects  due  to  heat  treatment.  The  pyrolysis  of  Si2ll6  over  Ai203 
will  be  carried  out  in  the  presence  of  CIl3SiD3  to  determine  whether  the  Af^OS  has 
changed  the  mechanism  of  decomposition  from  the  1,2-hydrogen  shift  observed  in 
the  homogenous  decomposition. 

- "^‘“2  ^  ^*”4  •  ■  -  (2) 


The  possibility  of  SiH^  decomposition  initiated  by  H  atoms  from  molecular  hydrogen 
decomposition  above  GOO  C  will  also  be  examined. 


All  of  these  reactions  will  be  carried  out  in  flow  systems  with  low-temperature 
traps  which  will  remove  the  heavier  less  stable  products  from  the  stream.  This  will 
allow  for  the  determination  of  the  effect  of  pyrolysis  temperature  on  the  yield  of 
Si2H(j  in  the  pyrolysis  of  8114; 


2  Silt,,— SijH,,  t  Ilj. 


(3) 


This  last  point  will  provide  some  information  concerning  the  species  which  are 
present  at  the  point  of  Si  deposition. 

The  first  experiments  to  be  undertaken  are  directed  toward  determining  the 
role  of  the  A£203  surface  in  catalyzing  the  decomposition  of  SiH4  and  other  Si 
hydrides.  Apparatus  has  been  fabricated  and  assembled, and  experimental  work 
is  in  progress. 

4.  SUBTASK  4.  PREPARATION  AND  CHARACTERIZATION  OF  SUBSTRATES 

In  the  first  year  of  the  program  it  was  demonstrated  that  A£203  surfaces 
prepared  by  mechanical  polishing  techniques  and  used  routinely  for  semiconductor 
heteroepitaxy  typically  have  severe  surface  and  subsurface  damage,  with  many 
scratches  often  several  microns  deep  and  rendered  invisible  to  close  inspection 
by  amorphous  or  fine-grained  debris  embedded  in  the  scratches  in  the  final  polishing 
stages.  Some  improvement  in  mechanical  polishing  procedures  was  achieved  in  terms 
of  the  density  and  depth  of  such  damage.  Gas-phase  etching/polishing  procedures 
using  SFg  and  various  fluorinated  halocarbons  in  the  1350  to  1500  C  temperature  range 
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produced  essentially  scratch-free  surfaces  on  (0112)  and  near-(1120)  Af20o  substrates. 
txteMiye  gas-phase  etch-rate  data  were  obtained  as  a  function  of  crystallogrt^ohic 
orientation  in  this  temperature  range. 

second  year  a  much  improved  technique  for  polishing 
K?  ■  j  was  developed,  and  excellent  surfaces  in  this  orientation  can  now  be 
Obtained.  Gas-phase  etching/polishing  techniques  were  further  developed  and 
exploited  for  (1)  thinning  .\f2^3  substrates  to  thicknesses  the  order  of  1  mil; 

(2)  evaluating  the  effects  of  prolonged  etching  on  (01i2),  (0001),  and  ~al2o’)  AfoO  • 

and  (3)  assessing  the  subsurface  damage  caused  by  various  mechanical  polishing  ^ 
procedures.  Routine  characterization  of  substrate  surfaces  at  various  stages  of 
preparatio^n  continued  using  techniques  of  X-ray  and  electron  diffraction  analysis  and 
optical  and  electron  microscopy. 

In  the  past  six  months  the  improved  polishing  procedure  (Ref  3)  for  (10l4)- 
oriented  AfgOo  has  been  verified  as  a  reproducible  process;  polishing  methods  for 
preparing  MgA<20^  surfaces  have  been  applied  to  wafers  used  for  Si  heteroepitaxy 
experiments;  considerable  effort  has  been  devoted  to  the  development  of  ion-beam 
sputterir^  techniques  for  thinning  Af203  wafers  for  use  in  the  in  situ  CVD  studies 
(Subtask  5);  gas-phase  etching  experiments  with  MgAf204  surfaces  were  initiated- 
experiments  with  Rutherford  backscattering  of  charged-particle  beams  as  a  mean’s 
of  characterizing  substrate  surface  conditions  were  resumed;  and  preliminary 
experiments  with  Auger  electron  spectroscopy  were  undertaken  to  evaluate  the 
techmque  as  a  means  for  characterizing  substrate  surface  impurities.  A  number 
of  the  investigations  planned  for  this  subtask  in  this  report  period  were  started 
late  because  of  increased  emphasis  on  other  subtasks  of  the  program;  only  the 
first  three  of  the  above  activities  will  be  discussed  further  below. 

a.  Mechanical  Polishing  of  (I014)-Oriented  AfgOg 

Si  deposition  experiments  using  (1014)  Af203  substrates  prepared  by  the 
improved  polishing  procedure  (Ref  3)  for  this  orientation  indicated  that  better 
films  can  now  be  obtained  on  these  improved  surfaces  than  had  been  realized 
previously.  This  is  consistent  with  the  observation  that  gas-phase  etching  in 
freon  at  1400  C  revealed  a  reduced  amount  of  surface  and  subsurface  damage  on  these 
substrates  relative  to  others  prepared  by  the  earlier  polishing  methods.  The 
improved  procedure  was  shown  to  be  a  reproducible  process;  several  groups  of  wafers 
in  this  orientati^on  were  prepared  in  the  same  v-ay,  and  Si  film  deposition  experiments 
indicated  that  the  surfaces  were  consistently  the  best  yet  achieved  on  this  orientation. 

b.  Polishing  Methods  for  MgAfgO^ 

Polishing  methods  being  developed  on  a  separately-funded  program*  for  a 
variety  of  crystallographic  orientations  of  stoichiometric  spinel  (MgA/„0.)  were 
applied  to  some  of  the  substrates  used  for  Si  heteroepitaxy  experiments  in  this 
contract  program. 


*  Subcontract  for  Crystal  Technology,  Inc. ,  Mountain  View,  CA.  as  part  of  AFML 
Contract  No.  F33615-70-C-1544. 
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c.  Ion- Beam  Sputtering  Techniques  for  Thinning  A<„0„  Substrates 

Ion-beam  sputtering,  also  called  ion-beam  machining,  ion  thinning,  ion 
sputtering,  or  ion  etching,  is  a  technique  whereby  material  is  removed  from  the 
surface  of  a  sample  by  controlled  sputtering  caused  by  an  incident  ion  beam.  The 
process  may  be  used  for  (1)  etching,  to  clean  substrates  prior  to  epitaxial  deposition; 
(2)  polishing,  in  which  a  scratched  or  work-damaged  surface  layer  typically  1  to  lOum 
deep  is  removed;  or  (3)  thinning,  in  which  bulk  samples  prepared  to  the  thickness 
limits  of  conventional  mechanical  polishing  are  thinned  to  ~200  A  as  required  for 
transmission  electron  microscopy.  Although  both  (1)  and  (2)  have  been  investigated, 
the  emphasis  to  date  Eias  been  placed  on  producing  thin  A<j>03  "wafers"  or  membranes 
for  use  as  substrates  for  the  in  situ  chemical  vapor  deposition  studies  of  this  contract 
(Subtask  5). 

The  ion  beam  sputtering  is  performed  in  an  Edwards  ion-beam  machining 
apparatus  (IBMA),  shown  in  the  photograph  of  Figure  4.  It  consists  of  a  conventional 
h^h-vacuum  system  containing  two  opposed  ion  guns  and  appropriate  feedthroughs. 

The  sample,  normally  0. 125  inch  in  diameter,  is  held  as  shown  in  the  schematic 
di&gratn  of  Figure  5,  so  that  each  side  is  stz*uck  by  an  ion  beam  at  a  glancing  angle 
of  5  to  30  deg.  Argon  ions  are  used,  with  energies  typically  5-10  kV.  The  sample 
is  rotated  in  its  plane  at  15  rpm.  Appropriate  shutters  are  used  to  measure  the 
ion-beam  current. 


Figure  4.  Ion  Beam  Machining  Apparatus  (IBMA) 
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The  IBMA  has  been  extensively  modified  to  overcome  certain  operational 
difficulties  and  to  improve  the  operating  convenience.  A  new  sample  holder  was 
built  to  hold  specimens  up  to  0. 75  in.  in  diameter  and  to  permit  tilting  the  sample 
for  a  beam  angle  down  to  5  deg  instead  of  the  standard  10  deg. 

A  new  gas  metering  system  was  built  by  replacing  the  original-equipment  leak 
valves  with  more  precise  vernier  metering  valves  for  control  of  the  low  gas  flow 
ranges  required.  Vacuum-tight  on-off  toggle  valves  w-ere  also  installed  so  the  gas 
flow  could  be  shut  off  without  disturbing  the  vernier  valve  setting. 

A  new  viewing  flange  has  been  built  and  installed  on  the  top  of  the  working 
chamber  (Figure  S).  This  and  other  modifications  permit  more  accurate  and 
convenient  operation  in  the  following  w'ays: 

1.  A  larger  recessed  window  permits  viewing  the  sample  with  a  stereo 
microscope  during  ion  bombardment  even  at  large  sample  tilt  angles. 
Considerable  time  is  saved  as  the  flange  no  longer  must  be  removed 
periodically  for  sample  examination  during  an  extended  thinning 
operation. 

2.  A  shield  protects  the  ion-current  measuring  shutters  from  all  but  the 
central  portion  of  the  beam.  Only  the  current  actually  used  for  ion-beam 
thinning  is  measured,  not  the  total  current. 

3.  A  movable  shutter  prevents  the  window  from  being  coated  with 
material  sputtered  from  the  sample  and  Its  holder. 

4.  A  larger  liquid-nitrogen  shield  was  installed,  increasing  the  shielding 
efficiency. 

The  first  step  in  preparing  Af203  for  transmission  electron  microscopy  is 
mechanical  lapping  and  polishing  of  the  wafer  to  as  small  a  thickness  as  possible, 
typically  0.002  to  0.003  inch.  Below  this  thickness  the  samples  fracture  during 
polishing  or  removal  from  the  polishing  jig  because  of  the  large  residual  stress 
present.  The  thinned  slices  are  then  cut  into  0. 125  inch  or  smaller  circular  discs, 
a  process  in  which  some  additional  samples  may  fracture.  Because  the  IBMA 
thinning  rate  is  so  low,  considerable  effort  has  been  expended  to  develop  a  procedure 
for  mechanically  polishing  At203  to  less  than  0.  002  inch,  but  the  yield  has  been  quite 
low.  A^203  wafers  eommereially  polished*  behave  similarly,  fracturing  during 
polishing  or  when  being  eut  into  disks. 

The  small  thinned  disks  are  installed  in  the  IBMA  sample  holder  sandwiched 
between  two  75-mesh  grids.  The  0.0007-inch  thick  grid  protects  the  Af203  below 
the  grid  bars  from  being  sputtered  until  the  grid  itself  is  sputtered  away.  The  Af203 
sample  after  thinning  has  a  "waffled"  surface,  with  the  ridges  providing  mechanical 
strength  so  that  very  thin  regions  will  not  break  during  handling.  The  grid  serves 
a  second  purpose  by  providing  a  number  of  individual  thinned  areas  rather  than  one 
larger  thinned  region,  thus  increasing  the  yield  by  a  factor  of  4  to  100.  Thinning 
usually  is  continued  until  a  perforation  is  noted  in  one  area;  then  the  sample  is 
examined  in  the  transmission  electron  microscope.  Ideally,  thinning  proceeds  until 
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a  suitable  electron-transparent  region  first  appears.  It  is  impossible  to  judge  with 
an  optical  stereo  microscope  the  precise  moment  to  stop  the  thinning  process,  so  in 
practice  it  is  continued  until  a  small  hole  is  visible  in  one  grid  square.  At  this  point 
some  grid  squares  will  have  only  thicker  regions,  while  other  grid  squares  will  have 
sizable  holes. 

'Hie  appearance  of  a  sample  thinned  beyond  the  customary  stopping  point  is 
shown  in  Figure  6a,  photographed  with  monochromatic  light  to  reveal  thickness 
contours  by  means  of  interference  fringes.  Since  this  sample  was  only  5.5 
thick  initially,  perforation  occurred  before  the  shielding  grid  was  completely 
sputtered  away. 

Three  different  orientations  of  have  been  successfully  thinned  by  ion- 

beam  sputtering  for  transmission  electron  microscopy  -  (0001),  (1014),  and  (0112). 
With  the  exception  of  samples  1,  2  and  8,  for  which  complete  thinning  was  not 
attempted,  each  sample  listed  in  Table  X  provided  thin  (~200  X)  areas  which  could 
be  examined  directly  and  which  were  at  least  marginally  suitable  for  use  as  sub¬ 
strates  for  the  in  situ  CVD  experiments.  The  electron  diffraction  spot  patterns  and 
Kikuchi  patterns  obtained  from  samples  in  these  three  orientations  are  shown  in 
Figures  7  through  11.  It  is  believed  that  these  are  the  first  such  patterns  obtained 
from  single-crystal  Af20Q.  The  spot  patterns  are  from  thin  regions  (<500  A  thick); 
the  Kikuchi  patterns  are  from  thicker  regions  (<10,000  A). 

Electron  micrographs  representative  of  the  samples  of  Table  X  are  shown  in 
Figures  12  through  20.  The  samples  are  of  relatively  high  quality  single-crystal 
A#20.j,  so  no  internal  defect  structure  is  visible.  The  ideal,  uniformly -thinned 
crystal  would  appear  featureless  on  an  electron  micrograph.  However,  three 
characteristics  are  clearly  visible  in  these  figures  -  extinction  contours,  surface 
sti’uctures,  and  a  thickness  gradient.  The  combination  of  these  factors  is  least 
noticeable  in  Figures  17,  18,  and  19,  which  accordingly  w'ould  be  the  most  desirable 
substrates  for  in  situ  CVD.  The  other  samples  would  be  marginally  suitable  as 
substrates.  The  importance  of  each  of  the  above  three  factors  will  be  discussed, 
and  then  other  features  of  specific  interest  in  certain  of  the  micrographs  will  be 
described. 

The  extinction  contours,  observable  in  all  of  the  micrographs  but  especially 
pronounced  in  Figures  12,  13,  14,  and  20,  are  the  result  of  constructive  and 
destructive  interference  of  the  electron  waves  passing  through  the  crystal.  Extinction 
contours  in  electron  microscopy  are  analogous  to  interference  fringes  in  optical 
microscopy,  but  unlike  the  latter  Iheir  spacing  is  not  simply  related  to  the  thickness 
of  the  crystal.  Slight  tilts  of  the  crj’sia!  will  dramatically  change  the  contour  spacing. 
Extinction  contours  are  quite  sensitive  to  crystal  orientation  and  thus  will  reflect 
minor  amounts  of  strain  or  buckling  in  the  crystal.  When  those  characteristics  are 
absent,  as  in  Figure  13,  the  extinction  contours  give  a  qualitative  indication  of 
local  thickness  variations.  The  combination  of  a  thickness  gradient  and  local  strains 
and  buckling  can  result  in  rather  interesting  combinations,  as  show'n  in  Figures  12 
and  14. 
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Figure  (>.  Optical  Micrograph  of  Thinncfl  f(lOOl)  Afot*;)  (Hample  No.  4, 
Sodium  l.ight)  fa)  25 X,  (b)  -l.lOX 

:)r, 


Samples  Thinned  by  Ion-Beam  Sputtering. 
(All  samples  are  Af-O  unless  noted  otherwise) 


Figure  7.  Electron  Diffraction  Spot  Pattern  from  (0001)  (Neg.  1362) 


Figures,  Electron  Diffraction  Kikuchi  Pattern  from  (0001)  Af  O  (Neg.  1360) 
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Figure  10.  Electron  Diffraction  Kikuchi  Pattern  from  (0112)  (Neg.  1510) 

3S 


Figure  12.  Electron  Micrograph  of  (0001)  A Sample  No.  3  (Neg.  1346) 
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Figure  14.  Electron  Micrograph  of  (1014)  A^20g,  Sample  No.  7  (Neg.  1450) 
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Figure  16.  Electron  Micrograph  of  (1014)  Sample  No.  9  (Neg.  1455) 
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Figure  18.  Electron  Micrograph  of  (10l4i  Af  O  ,  Sample  No.  9  (Neg.  1517) 
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Figure  19.  Electron  Micrograph  of  (0112)  Sample  No.  10  (Neg.  1514) 


Figure  20.  Electron  Micrograph  of  (0112)  Sample  No.  10  (Neg.  1515) 
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It  is  emphasized  that  the  circular  semisymmetrical  patterns  in  these  fiffures 
are  not  objects  in  the  crystal  but  rather  Intersections  of  the  extinction  contours.  This 
can  be  readily  demonstrated  by  tilting  the  sample  slightly  duiing  examination-  the 
circular  Intersections  then  move  about  the  crystal.  The  obscuring  presence  of  extinc¬ 
tion  contoura  IS  best  eliminated  by  producing  uniformly  thin,  strain-free,  unbuckled 
crystals  by  the  ion-beam  sputtering  process. 

Certain  isamples  exhibit  a  surface  structure  after  thinning;  this  is  particularly 
noticeable  in  Fibres  17  and  20.  This  semiregular,  wavy  or  terraced  structure  is 

orientations,  but  not  on  (0001)  samples,  which  suggests 
crysta  lographic  origin  However,  their  direction  in  Figure  17  does  not  appear  to  be 
crystallographic,  and  tte  structure  does  not  appear  on  all  samples  or  in  all  thin  areas 
of  the  same  sample.  The  surface  structure  is  more  pronounced  in  areas  where  the 
thickness  gradient  is  least,  suggesting  a  complex  interaction  between  the  crystalline 
sample  and  the  ion-sputtering  mechanism.  The  structure  does  not  appear  to  be  simply 

related  to  ion-be^  energy  or  incidence  angle.  The  origin  of  the  struchire  is  currently 
being  investigated.  ^ 


Many  regions  in  the  samples  examined  to  date  have  a  thickness  gradient  that  is 
larger  than  desirable  and  hence  are  marginally  suitable  as  substrates  for  in  situ  CVD 
Thin  regions  with  suitable  large  areas,  such  as  shown  in  Figures  17  through  19  are 
seldom  observed  more  than  once  per  sample.  ' 


The  magnitude  of  the  thicknesi  gradient  has  been  analyzed  by  optical 
interferometry.  The  central  thin  region  of  a  sample  appears  as  in  Figure  6b,  when 
examined  by  monochromatic  sodium  light,  The  asyisimetrical  shape  is  the  result  of 
the  bottotn  grid  being  slightly  misaligned  with  respect  to  the  top  grid  during  thinning, 
a®  In  Figure  6a.  The  thickness  variation  has  a  slope  of  1:50,  which  appears  to  persist 
down  to  the  thinnest  regions,  typified  in  Figures  12,  15  and  16.  The  thickness  will 
vary  from  zero  to  200  A  over  a  transverse  distance  of  ~l|jLm.  In  practice,  the 
thinnest  regions  (0  to  50  A)  near  a  hole  break  off  during  handling. 


P^ential  source  of  an  excessive  thickness  gradient  is  obvious  from  Figure  6a 
^e  grid  partially  sWelds  some  areas  of  the  sample,  the  exact  area  being  determined 
by  the  incidence  angle  and  the  thickness  of  the  grid.  If  a  hole  becomes  enlarged  and 
the  region  (^amin^  is  in  the  shadowed  area,  the  thickness  gradient  will  naturally  be 
greater.  This  problem  can  be  avoided  simply  by  limiting  electron  microscopic 
examination  to  the  unshadowed  center  of  the  grid  square. 

The  area  of  the  sample  always  exposed  to  the  ion  beam  within  a  grid  square 
during  the  tlunning  process  may  be  calculated  as  follows.  The  75-mesh  grid  has  a 
^peat  spaci^  of  338  nm,  with  open  areas  268  pm  wide  not  covered  by  the  grid  bars, 
^e  height  of  the  grid  (18  pm)  causes  some  of  the  sample  to  be  shadowed  a  portion  of* 
the  time.  At  a  beam  angle  of  30  deg,  the  area  never  in  shadow  is  205  pm  x  205  pm 
or  59  percent  of  the  open  area.  At  a  beam  angle  of  10  deg  the  corresponding  figure 
IS  5  percent  If  for  some  reason  the  grid  is  not  in  contact  against  the  sample  but  is 
instead  raised  another  18  pm,  for  example,  there  is  no  region  within  a  grid  square 

that  would  be  exposed  at  all  times  to  a  beam  incident  at  10  deg.  The  figures  for  other 
angles  are  listed  in  Table  XI.  b  b  umei 
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It  is  possible  that  a  limiting  but  relatively  steep  thickness  gradient  is  imposed 
on  the  sample  by  the  ion  sputtering  process  itself.  Thin  areas  are  more  rapidly  eroded 
than  thick  areas,  and  this  fact  alone  may  set  the  lower  Umit  on  the  gradient  attainable. 

A  ''thin”  area  would  be  one  in  which  the  sample  thickness  is  the  same  order  of  magni¬ 
tude  as  the  ion  penetration  depth.  The  liklihood  of  a  steep  limiting  gradient  from  the  ion 
sputtering  process  itself  is  less  because  adequately  thinned  regions  have  indeed  been 
obtained  (e.g. ,  Figures  17  ard  i9). 


One  feature  found  on  both  optical  and  electron  micrographs  of  Af203  samples 
has  been  conclusively  identified  as  an  artifact  of  the  ion  etching  technique,  and  that  is 
the  so-called  protuberances.  These  approximately  hemispherical  bumps  are  visible  in 
the  optical  photomicrographs  of  Figures  6  and  21  and  the  electron  micrograph  of 
Figure  13.  These  features  range  from  50  to  fractions  of  a  m  'ron  in  diameter,  and 
They  are  caused  by  some  form  of  contamination  shielding  the  sample 
surface  from  the  ion  beam,  thus  preventing  sputtering  in  that  location.  Figure  21 
clearly  shows  that  one  source  of  such  contamination  is  the  wire  grid  used  during 
thinning.  Protuberances  \wll  form  even  if  a  grid  is  not  used,  however,  for  the 
Sample  holder  is  an  unavoidable  source  of  contamination.  Protuberances  are  un- 
desirable  if  present  in  large  quantities  because  they  create  an  irregular  surface  on 
the  thinned  sample.  Their  incidence  can  be  reduced  but  not  entirely  eliminated  by 
thorough  cleaning,  appropriate  shields,  and  use  of  a  small  incidence  angle.  It  is 

emphasized  that  the  "objects"  visible  in  Figures  12  and  14  are  not  protuberances  but 
are  intersections  of  extinction  contours. 


Reproduced  from 
best  dveilable  copy. 


Figure  21.  Optical  Photomicrograph  of  Thinned  (10l4)  Ai„0„ 
(Sample  No.  6,  Sodium  Light,  IlOX)  ^  ^ 
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The  In  the  micrograph  of  aample  9  (Figure  18). 

in  the  4o3  aJl  a^o^  ph^ae  ClevClurfu"”  orlglnally''preaem 

resting  on  the  surface  is  determin^  bv  exam^n^fV.  Particulate  matter 

around  the  inclusion;  at  higher  mrenifk-adJ^!.  fh?  the  extinction  contours  in  the  ^203 

presence  la  revealed.  Thf  "fuz”y!?  ap^\“ree  tK-^ 

focus  but  rather  to  a  eoutamlnatlon  layS  Sch  hid  talU  u^l’K  "“'.‘'“.f 

prolonged  time  it  was  examined  Residual  ^  sample  during  the 

were  polymerised  by  the  electron  beam  Into  I  carbomTeSTpUr 

rates  from  0. 1  '^ariablM  TffSt-*'*<h*  **  “  PTnann®-  with  removal 

malerlal,  crystallograSc  oriental™  ^  ^™P'n 

respect  to  sample,  alignment  of  beam  and  samnhn 

limited  number  of  samples  thi^ed  fL.?f  ^  stability  of  the  ion  beam.  The 

preliminary  iesult.s  suggest  the  differences”  dllr  o'’th’e  toUhree''v  >>“< 

of  major  importance.  three  variables  are  not 

reduce?tMX'?a”I1,  leiralLles ‘‘rfll  a  definite  trend  touurd 

presence  of  profuberanles  oTthe^emnl^  'f  ■“  *■'  >t>inimize  the 

thls^has  the  practical  disadvantage  of  t..reasle;^:LI*^bTy^hTI^^^ 

(sample  2  In  Table°X)*'*  The  beam  has  c””"  tiptermlned  by  thinning  a  large  AfgO.  wafer 
along  the  axis  of  the  Ion  gun.  The  J  3  cTdf.*  “°®'  ‘"tense^reglon 

after  thinning,  with  a  3. 5-mm  dll  um  Zn  ’  ®®"®®t-®''aped  surface 

apering  to  10  pm  deep  near  the  edees^  In^h^ center, 

1.4  pm/hr,  while  at  L  edgis  the  fate  ™® 

IndlcaUnga  misalignment  bfthifaraount  of  renter  of  the  sample, 

the  axis  of  the  ion  beam.  For  most  samnlpa  rotation  of  the  sample  and 

for  samples  smaller  than  1  mm  dia  it  meins  the*sTmile^"TiT  “"^"’Portant;  hovwver, 
portion  of  the  ion  beam  for  a  fraction  of  the  tim#.  ^  ^®  intense 

thinning  rate  lower  than  normal.  uring  rotation,  making  the  effective 

determines  the  degree^tlTlwchlhe  mkc^iM  important  factor,  for  it 

sputtering  occurs  not  only  from  the  samole  but  As  thinning  proceeds, 

cathode  as  well,  and  the  U^red  material  is^  *1°^"®^.*  "®^'«  shields,  and 

In  particular,  material  from  the  cathorio  elsewhere  in  the  system, 

causing  a  short  circuit  and  automatic  shut-off  of^e  cathode, 

during  attended  operation  the  machine  pan  Ko  f  apparatus.  When  this  occurs 
5  to  120  min  before  again  arcil  °P®r®te  for 

lost.  Arcing  can  be  prevented  by  thoroughly  clea^ne^he^^o^^^'”"  operating  time  is 
requires  ~  1  hr  of  down-time  for  everv  20  to  In  ^  I!*®  T  this 

duty-cycle"  must  be  selected.  '  operating  hours,  an  economical  operating 
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Sputtering  of  the  cathode  occasionally  generates  a  more  serious  problem.  The 
original  cylindrical  beam  exit  canal  normally  becomes  conical  in  operation  due  to  the 
erosion  by  the  beam  but  it  may  become  skewed.  In  this  case  the  most  intense  portion 
of  the  ion  beam  does  not  strike  the  sample  and  the  thinning  rate  is  reduced.  The  reduced 
thinning  rate  of  sample  9  (Table  X)  is  partially  the  result  of  this  effect.  The  cathode 
canal  normally  increases  in  size  until  the  cathode  must  be  replaced  (typically  after 
100  to  200  hr  of  operation),  but  replacement  may  be  required  earlier  if  the  canal  becomes 
skewed. 

5.  SUBTASK  3.  STUDIES  OF  IN  SITU  FILM  GROWTH  IN  THE  ELECTRON 
MICROSCOPE 

In  the  first  year  of  the  program  many  of  the  modifications  required  in  the 
electron  microscope  for  iii  situ  observation  of  the  nucleation  and  early-stage  growth 
of  eVD  semiconductor  films  on  insulating  substrates  were  completed.  Provision  for 
motion  picture  recording  of  film  growth  was  assembled  and  tested,  and  the  heated 
specimen  stage  was  installed  and  tested.  The  first  in  si^  PVD  experiments  w-ere 
carried  out.  During  the  first  half  of  the  second  year  a  series  of  electron  microscope 
modifications  and  tests  was  completed,  culminating  in  a  series  of  successful  PVD 
experiments  inside  the  electron  microscope.  \t  was  deposited  onto  a  heated  carbon 
substrate  and  a  sequence  of  micrographs  taken  during  the  growth  process,  demon¬ 
strating  the  feasibility  of  performing  in  situ  nucleation  and  growth  studies  in  this 
equipment. 

A  transmission  phosphor  screen  (for  the  motion  picture  camera)  was  installed, 
permitting  motion  picture  photography  which  does  not  interfere  with  the  normal 
still  photography.  The  auxiliary  vacuum  pumping  system  for  the  specimen  chamber 
was  fabricated,  installed  and  tested.  The  basic  vacuum  system  of  the  microscope 
itself  was  improved  by  addition  of  a  cooled  baffle,  by  polishing  numerous  O-ring 
grooves,  and  by  thoroughly  cleaning  the  microscope  interior.  A  PVD  source  assembly 
was  fabricated,  installed,  and  used  in  conjunction  with  the  specimen  heater  to 
perform  the  above  PVD  experiments.  Calculations  and  design  for  the  CVD  micro¬ 
chamber  were  completed  and  fabrication  of  the  chamber  and  the  differential  pumping 
apertures  was  begun.  Work  was  also  begun.on  the  production  of  thin  (200  A) 
crystals  required  as  substrates  for  the  in  situ  CVD  studies.  ’’ 

Effort  oh  this  subtask  has  been  divided  betw’een  two  major  areas  during  the 
past  six  months:  (1)  design  and  construction  of  the  CVD  microchamber  and  associated 
hardware,  and  (2)  conducting  a  large  number  of  PVD  experiments  to  test 
and  demonstrate  various  experimental  situations  to  be  encountered  in  the  in  situ 
chemical  vapor  deposition  experiments.  A  new  design  for  the  CVD  microchamber  was 
completed,  superior  in  many  respects  to  that  described  in  the  previous  semiannual 
report  (Ref  .‘1).  Fabrication  and  assembly  of  the  microchamber  is  essentially  complete, 
with  testing  scheduled  for  the  first  month  of  the  third-year  program.  A  CVD  flange 
to  provide  the  gas  inlet  and  electrical  feedthroughs  for  the  microchamber  has  been 
fabricated  and  leak  tested.  i 

Numerous  in  situ  PVD  experiments  have  been  conducted  during  this  period, 
with  both  Af  and  Au  being  deposited  onto  amorphous  carbon  substrates.  These 
experiments  are  still  not  fully  satisfactory,  so  have  not  been  documented  on  tnotion 
picture  filr  i.  The  major  problems  have  been  the  result  of  shortcomings  in  the  PVD 
heating  stage,  coupled  with  a  contamination  problem  in  the  microscope  and  the  lack 
of  a  suitable  transmission  phosphor  screen. 
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a.  Design  and  Construction  of  in  situ  CVD  Microchamber 


The  present  design  of  the  CVD  microchamber  is  shovv'n  in  a  cutaway  drawing  in 
Figure  22  and  a  photograph  in  Figure  23.  The  substrate  is  held  between  two 
resistively  heated  grids  which  are  held  by  screws  in  contact  with  inner  and  outer 
electrodes.  The  sample  and  the  grids  are  totally  enclosed  in  a  thin  cylindrical  region, 
vacuum-tight  except  for  the  gas  inlet  and  two  apertures  (for  the  electron  beam 
entrance  and  exit).  The  Sill^  gas  enters  the  microchamber  through  a  coiled  gas  inlet 
tube,  surrounds  the  sample,  and  exits  via  the  two  apeitures.  The  microchamber  has  a 
gas  layer  thickness  of  l.G  mm  through  which  the  electron  beam  must  pass. 

Certain  constructional  details  (Figure  22)  of  the  CVD  microchamber  are  of 
interest.  The  center  electrode  contains  a  0.004-in.  drilled  hole  to  permit  the  electron 
beam  to  enter  and  a  limited  amount  of  gas  to  exit.  A  raised  projection  contains  a 
tapped  hole  for  the  No.  40  National  Miniature  (NM)  screw  (O.OlG-in.  dia,  254  threads 
per  in. )  which  holds  one  end  of  the  heating  grid.  The  center  electrode  is  brazed  into 
a  base  which  is  an  extensively  modified  AEI-EMG  electron  microscope  specimen  holder. 

One  end  of  a  coiled  copper  gas  inlet  tube  is  brazed  into  the  center  electrode  base, 
and  the  other  end  is  joined  to  the  CVD  flange  by  a  ’’Swaglok"  tubing  connector.  The 
flat  portion  of  the  base  is  slotted  to  permit  gas  to  flow  from  the  inlet  tube  to  the  annular 
space  betw-een  the  inner  and  outer  electrodes  and  then  to  the  specimen.  The  center 
electrode  is  beryllium  copper,  selected  for  its  superior  machining  characteristics 
over  those  of  pure  copper,  which  is  unsuitable  for  drilling  and  tapping  the  miniature 
holes  required.  As  anticipated,  some  difficulty  occurred  in  brazing  the  beryllium 
copper  to  the  pure  copper  base,  but  the  problem  has  been  adequately  solved’ 

l^e  outer  electrode  is  attached,  by  means  of  three  No.  80  NM  flat-head  screw'S 
(0.032  in.  dia,  127  threads  per  in. ),  to  the  lava  insulator  which  in  turn  is  fastened  to 
the  base  with  four  No.  80  N-Vl  fillister-head  screw’S.  The  lava  insulator  was  machined 
somew’hat  oversize  in  the  urfired,  soft  condition,  then  fired,  ground  flat,  and  lapped 
to  a  2-microinch  surface  finish.  The  bottom  of  the  outer  electrode  and  the  flat  portion 
of  the  center  electrode  base  were  also  lapped  to  a  2-microinch  surface  finish  so  that  the 
assembly  is  gas  tight  at  these  joints  when  assembled.  This  portion  of  the  microchamber 
need  not  be  disassembled  to  load  or  unload  specimens. 

The  outer  electrode  has  a  projection  containing  a  small  tapped  hole  which  accepts 
the  No.  40  NM  screw  holding  the  other  end  of  the  heating  grid.  A  threaded  cap  over 
the  outer  electrode  holds  a  thin  removable  copper  aperture  disc  in  place.  The  outer 
electrode  and  cap  are  made  of  Ti,  selected  for  its  machinability,  nonmagnetic 
properties,  and  low  vapor  pressure  at  elevated  temperatures. 

It  was  originally  anticipated  that  stainless-steel  screws  would  be  used  to  attach 
the  heating  grid  to  the  electrodes.  How-ever,  use  of  similar  screws  in  the  PVD 
experiments  has  revealed  a  small  trace  of  ferromagnetic  material  in  the  essentially 
nonmagnetic  screws;  this  creates  a  distorted  magnetic  field  in  the  vicinity  of  the 
specimen  and  prevents  attainment  of  maximum  resolution  of  the  electron  microscope 
above  20,000X.  Various  heat  treatments  have  failed  to  reduce  the  ferromagnetic 
content  of  the  stainless-steel  screws,  so  specially  machined  Ti  screws  have  been  made. 
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Figure  22.  In  Situ  Chemical  Vapor  Deposition  Microchamber. 

(A,  center  electrode  and  gas  limiting  aperture;  B,  lava  insulator;  C,  outer  electrode; 
D,  heatit^  grids;  E,  sample;  F,  removable  aperture;  G,  aperture  cap.) 
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Figure  23.  Microchamber  for  in  situ  CVD  Experiments 
b.  In  situ  PVD  Experiments  with  A< 

The  A£  PVD  experiments  initiated  in  the  first  six  months  of  the  second  year 
have  been  continued.  The  At  is  evaporated  from  a  specially  constructed  filament 
mounted  inside  the  microscope  onto  an  amorphous  carbon  film  heated  to  ~300  C  by 
a  resistively-heated  stainless-steel  support  grid.  A  sequence  of  electron  micro¬ 
graphs  taken  during  the  course  of  one  such  ^  situ  experiment  is  shown  in  Figure  24a 
through  24f.  The  first  photograph  in  the  series  is  that  of  the  original  substrate  prior 
to  nucleation;  the  particles  visible  are  dust  particles,  not  At  nuclei.  They  act  as 
reference  points  so  that  exactly  the  same  area  can  be  located  in  the  subsequent 
photographs. 

Several  important  aspects  of  thin-film  growth  can  be  observed  in  these  micro¬ 
graphs.  Locations  A  and  B  in  Figures  24b  and  24c  illustrate  the  formation  of  a  grain 
boundary  in  a  large  crystal  formed  from  several  smaller  crystals.  Various  contrasts 
visible  at  A  in  Figure  24b  indicate  that  a  number  of  differently  oriented  crystals  are 
present.  (The  loss  of  resolution  in  reproducing  the  original  photographs  here  may 
make  discrimination  of  individual  nuclei  difficult,  but  they  are  readily  discernible  in 
the  original  micrographs. )  Only  three  of  the  orientations  survived  the  coalescence 
process;  the  remainder  not  visible  in  Figure  24c  have  been  absorbed,  atom  by  atom, 
by  their  neighbors.  Note  that  at  A  both  orientations  have  apparently  grown  at 
approximately  the  same  rate. 
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Figure  24.  Sequence  of  Electron  Micrographs  Recorded  during  At  PVD  at  Times  nf 

fr%\  A  /U\  f  f  n _  _ _ 


At  B,  it  is  seen  that  the  original  crystals  grew  at  different  rates,  resulting  in 
one  orientation  becoming  dominant.  The  dark  portion  of  the  top  crystal  in  Figure  24c 
is  almost  completely  developed  in  Figure  24b,  and  no  further  growth  has  occurred 
except  in  the  light  portion.  The  bottom  crystal,  how’ever,  show's  the  reverse  behavior. 
(All  four  crystals  are  of  different  orientations.)  Note  at  location  C  that  not  all  crystals 
have  grown  at  the  same  rate.  At  location  D,  the  large  extent  to  which  rearrangement 
occurs  during  film  growth  is  visible. 

Some  areas  of  the  substrate  originally  covered  with  AS  may  become  bare  when 
several  islands  coalesce  and  secondary  nucleation  w'ill  then  occur  in  these  bare  areas. 
Such  an  area  can  be  seen  immediately  outside  and  to  the  upper  left  of  location  B  in 
Figures  24c  and  24d.  The  presence  of  secondary  nuclei  causes  some  of  the  apparent 
"clutter"  in  the  photographs;  the  remainder  is  caused  by  oxidation  of  the  A^. 

A I  deposition  w’as  stopped  after  190  sec  in  this  particular  experiment  but  the 
substrate  and  deposited  film  were  permitted  to  remain  at  an  elevated  temperature  for 
an  additional  10  min.  Almost  total  oxidation  of  the  AS  occurred,  completely  changing 
the  char'acter  of  the  film.  All  the  O2  required  for  the  reaction  w’as  supplied  by  the 
residual  O2  and  water  vapor  present  in  the  electron  microscope  vacuum  of  <1  x  10“''>torr. 
Much  of  the  water  vapor  originates  in  the  photographic  plates  used  for  still  pictures,  a 
problem  which  can  not  be  eliminated  until  the  external  motion  picture  camera  system 
has  been  made  to  operate  satisfactorily.  In  view  of  the  serious  oxidation  problem 
with  Al  deposits,  no  further  in  situ  PVD  experiments  with  this  material  are  being 
performed. 

In  situ  PVD  Experiments  with  Au 

WTien  Au  w’as  deposited  inside  the  microscope,  an  interesting  effect  was  noted. 

The  nuclei  were  extremely  small  in  that  area  of  the  substrate  exposed  to  the  electron 
beam.  The  resulting  Au  film  w-as  extremely  fine-grained  in  this  region,  but  consider¬ 
ably  larger- grained  outside  the  irradiated  area.  Various  repeats  of  the  earlier 
experiments  yielded  the  same  final  results,  an  example  of  w'hich  is  shown  in  Figure  25. 
Figure  25a  shows  the  region  visible  during  the  experiment.  Figure  25b,  taken  after 
completion  of  the  experiment,  show's  an  adjacent  area  which  was  outside  of  the  field 
of  view'  and  not  struck  by  the  electron  beam  during  the  deposition  experiment. 

Figure  25c  is  a  lower  magnification  view  taken  with  the  electron  beam  expanded  to 
include  both  of  the  areas  show'n  in  a  and  b. 

In  addition  to  the  restricted  size  of  the  Au  nuclei,  the  substrate  region  irradiated 
by  the  electron  beam  was  darkened  by  a  deposit  of  hydrocarbon  contamination.  The 
deposit  is  formed  from  residual  hydrocarbon  vapors  in  the  microscope  striking  the 
substrate  and  becoming  polymerized  by  the  incident  electron  bcvom.  The  carbonaceous 
deposit  thus  formed  gradually  builds  up,  eventually  obscuring  the  sample  being 
examined.  Hydrocarbon  vapors  which  strike  the  sample  outside  the  area  irradiated 
by  the  electron  beam  will  reevaporate  without  effect.  These  residual  hydrocarbon 
vapors,  present  to  varying  extent  in  all  electron  microscopes,  originate  from  the 
vacuum  oils  and  greases  used  in  the  pumping  system,  sealing  gaskets,  organic  insula¬ 
tion  on  interior  electrical  wires,  fingerprints, aixl  residues  of  cleaning  solvents. 

•The  first  hypothesis  presented  for  the  unusual  nucleation  behavior  was  that,  for 
some  unknown  reason,  the  microscope  contamination  rate  had  increased  greatly  between 
the  final  Ai  PVD  experiment  and  the  first  Au  PVD  experiment.  The  diffusion  of 
individual  Au  atoms  on  the  sample  surface  -  and  thus  the  rates  of  the  nucleation,  growth 
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Figure  25.  PVD  Au  Film  Grown  in  Electron  Microscope. 

(a)  Region  Observed  during  Growth;  (b)  Adjoining  Region  not  Irradiated 
by  Electron  Beam;  (c)  Lower  Magnification  View  of  Regions  of  (a)  and  (b) 
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and  coalescence  stages  -  would  then  be  seriously  affected.  However,  this  hypothesis 
was  discarded  when  the  contamination  rate  was  remeasured  and  found  to  be  140  A/min, 
only  slightly  larger  than  previously  measured  values. 

Since  the  contamination  rate  decreases  rapidly  with  increasing  temperature,  it 
was  logical  to  question  if  the  substrate  was  actually  being  heated  to  the  intended 
temperature  of  ~500C  (at  which  temperature  the  hydrocarbon  vapors  adsorbed  on  the 
substrate  are  thermally  reevaporated  before  the  electron  beam  polymerizes  them). 
Examination  of  the  substrate  heater  revealed  a  partial  short  circuit,  which  had  resulted 
in  reduced  heater  grid  current  and  correspondingly  lower  temperature.  Repair  of  the 
heater  stage  reduced  somewhat  the  abnormal  nucleation  behavior  but  it  did  not 
completely  eliminate  the  effect. 

The  heating  stage  was  again  repaired,  this  time  to  improve  the  electrical  contact 
between  the  heating  grid  and  the  electrodes.  The  original  cap,  which  made  pressure 
contact  to  the  grid,  was  removed  and  No.  70  NM  screw  threads  were  tapped  directly  into 
the  electrodes.  Testing  in  the  electron  microscope  indicated  a  loss  of  resolution  which 
prohibited  usage  above  ~20,000X.  As  mentioned  earlier,  the  stainless-steel  screws 
apparently  contained  a  trace  of  ferromagnetic  material  which  distorted  the  electro¬ 
magnetic  field  of  the  objective  lens  around  the  specimen.  Various  heat-treating 
procedures  failed  to  reduce  the  ferromagnetic  content  appreciably,  so  nonmagnetic  Ti 
screws  are  being  used. 

Experiments  have  continued  with  the  same  results;  greatly  altered  Au  nucleation 
and  growth  behavior  still  occurs  in  regions  Irradiated  by  the  electron  beam.  Two 
hypotheses  are  currently  being  considered:  (1)  Au  nucleation  and  growth  are  intrinsically 
more  susceptible  to  contamination  effects  thanfar  Af ,  and  (2)  the  substrate  is  still  not 
making  adequate  thermal  contact  with  the  heating  grid.  The  first  suggestion  is 
intuitively  difficult  to  accept  and  spears  contrary  to  the  experimental  evidence  of 
related  nucleation  and  growth  studies.  Furthermore,  it  would  be  indicated  only  upon 
elimination  of  all  other  possible  causes.  Experiments  are  currently  underway  to 
test  the  second  hypothesis,  with  no  results  to  date. 

d.  Characteristics  of  Transmission  Phosphor  Screens 

The  transmission  phosphor  screen  is  an  integral  component  of  the  motion 
picture  recording  set-up  on  the  modified  AEI  EM6  electron  microscope.  It  is  a  glass 
p  ate  with  a  coating  of  an  electron-activated  phosphor  material  on  the  side  internal  to 
the  electron  microscope,  e.g. ,  the  vacuum  side.  To  date,  five  screens  from  three 
manufacturers  have  been  tested;  one  of  these  is  marginally  adequate  but  the  others 
have  been  unsuitable. 

The  first  four  screens  tested  used  a  "Pll"  phosphor  coated  with  a  thin  layer 
of  At  to  prevent  charging  by  the  electron  beam.  The  Pll  phosphor  was  chosen 
because  its  maximum  spectral  efficiency  occurred  in  the  blue  at  0.46  pm,  nearly 
the  wavelength  to  which  Kodak  RAR  2479  film  is  most  sensitive.  The  intensity  available 
on  the  screen,  even  when  struck  with  a  maximum-intensity  electron  beam,  was 
inadequate  to  obtain  an  image  on  the  film  using  a  Swltar  fl.  6  lens  and  a  shutter  speed  of 
1/30  sec.  Presumably  the  phosphor  layer  and/or  the  Af  coating  was  too  thick.  In 
addition,  the  resolution  was  less  than  desired,  due  either  to  a  large  phosphor  grain  size 
or  to  structure  in  the  Af  coating. 
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The  best  screen  tested  to  date  is  of  S2  green  F2  phosphor  material  coated  with 
a  conducting  layer  of  amorphous  carbon.  The  resolution  is  more  than  adequate.  The 
intensity  is  better  than  that  obtained  with  previous  screens  and  is  just  sufficient  to 
record  images  on  the  film  when  the  developing  process  is  "forced, "  i.e. ,  the  developing 
time  extended. 


The  camera  lens  and  film  being  utilized  are  the  best  commercially  available,  so 
solution  of  the  problem  appears  to  be  in  increasing  the  exposure  time  of  a  single  frame. 
This  has  the  unattractive  feature  of  requiring  a  more  elaborate  and  more  costly  (but 
commercially  available)  motion  picture  camera.  This  alternative  is  now  being  pursued. 

().  SUBTASK  ().  EVALUATION  OF  FILM  PROPERTIES 

During  the  first  year,  the  routine  evaluation  of  film  properties  was  earned  out 
by  established  methods  of  X-ray  and  electron  diffraction  analysis,  metallographie 
analysis,  and  electrical  measurements  of  transport  properties.  In  addition,  a  new 
technique  for  evaluating  the  characteristics  of  the  interfacial  region  of  hetcroepitaxial 
films  was  developed,  involving  measurement  of  the  photoelectron  emission  from 
monoehromatieally-illuminatcd  films  in  the  MIS  configuration  on  insulating  (viz. , 
A22O3)  substrates.  Relatively  large  photoeurrents  due  to  electron  transport  through 
thick  (-10  mils)  single -crystal  A^gO^  substrates  were  measured  as  a  function  of 
photon  energy.  Photoelectric  threshold  energies,  escape  length  (mean  free  path)  of 
excited  electrons,  and  band  bending  in  the  semiconductor  film  adjoining  the  interface 
were  determined  in  the  Si/A^2^'3  ^ind  GaAs/A^20o  systems.  Observation  of  the 
energy  spectrum  of  back-scattered  proton  or  alpha-partiele  beams  injected  in  chan¬ 
neling  directions  In  hetcroepitaxial  semieonduetor  films  was  also  investigated  as  a 
means  of  measuring  the  density  and  the  location  of  structural  defects  in  the  films. 
Experiments  indicated  that  Si/insulator  films  have  less  imperfect  interfacial  regions 
than  do  GaAs/insulator  films.  The  best  structures  of  those  examined  to  date  were 
found  in  GOO)  SI  films  on  (0112)  A^20g  substrates  and  in  (111)  Si  films  grown  on 
near-(l  120)  A22O.5  substrates. 


In  the  first  six  months  of  the  second  year  the  study  of  the  effects  of  changes  in 
deposition  parameters  on  Si/AggOy  film  properties  continued.  These  studies 
led  to  considerable  insight  into  the  factors  which  most  profoundly  influence  film 
quality,  and  identification  of  conditions  for  optimized  film  growth  on  the  various  AfgOo 
orientations  progressed  well.  The  iniportanee  of  reactor  geometry  was  reeognizeu 
and  demonstrated.  The  extent  of  A2  autodoping  from  the  substrate  was  established, 
and  appropriate  annealing  procedures  for  minimizing  these  effects  were  determined. 
The  use  of  new  A220g  orientations  for  Si  growth  led  to  substrates  which  yielded 
epitaxial  films  of  as  good  or  better  quality  than  those  previously  obtained.  The 
measurements  of  photoemission  of  electrons  from  hetcroepitaxial  semieonduetor  films 
and  transport  of  electrons  in  A^203  was  carried  further;  work  functions  of  additional 
metals  were  determined,  and  the  mechanism  of  electron  transport  through  the 
insulator  was  studied  further.  Measurements  of  high-field  transport  properties  of 
Si  and  GaAs  hetcroepitaxial  films  were  also  initiated  in  that  period. 

During  the  piist  six  months  a  major  part  of  the  contract  effort  has  been  expended 
on  this  subtask.  The  speeifie  areas  investigjited,  described  in  further  detail  below, 
include  (1)  a  study  of  the  variation  in  the  electrical  properties  of  Si/Af20g  with 
temperature;  (2)  an  extensive  study  of  the  anisotropy  of  the  electrical  properties  of 
Si  films  on  A220g,  involving  both  theoretical  (see  Subtask  1)  and  experimental  investi¬ 
gations;  (3)  an  evaluation  of  the  eleetrieal  properties  of  Si  films  on  MgA2204 
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substrates;  (4)  further  evaluations  of  the  eleetrical  properties  of  Si/A^20o  to  establish 
relationships  among  temperature,  growth  rate,  and  substrate'orientation; 

(5)  further  study  of  the  photoeleetrie  effeets  observed  in  the  Si/A^oOo  system,  particu- 
lany  in  the  AII2O2  substrates;  and  (6)  investigation  of  the  high-field  transport  properties 
of  heteroepitaxial  films. 


a. 


Variation  of  Eleetrical  Properties  of  Si/A£203  with  Temperature 


Ho  11  variability  of  a  computer  program  for  the  analysis  of  low  temperature 

Hall-effeet  data  has  provided  a  technique  for  deducing  information  about  electricallv 
and  defects  in  Si/A«203  films.  This  information,  coupled  with 
knowledge  of  the  carrier  scattering  mechanisms  in  Si  films  that  can  be  gained  from 
low -temperature  studies,  has  suggested  a  study  of  the  feasibility  of  using  low- 
temperature  data  as  a  measure  of  film  quality  and  a  means  to  assist  in  optimizing 

Si/A?  n detailed  examination  of  two  n-type 
S1/A2203  frlms,  has  been  completed  during  this  reporting  period.  The  electrical 

a^d's'orr'  Si/(01l2)A£2O3  films  having ^hickSesscT^^^^^ 

and  3. 9  pm  were  measured  between  room  temperature  and  ~40K.  The  growth 
f  substrates  were  identical  for  the  two  samples  (with  the  e^eption  of 

iauW  be  mTde  It  iif  comparison  of  thin  and  thick  Si  films 

could  be  made.  It  vas  expected  that  any  differences  in  electrical  behavior  in  the  two 
samples  could  be  attributed  to  the  difference  in  strain  and/or  defect  stmeyre  in  ^ 
tw  o  samples.  After  growth,  both  samples  were  annealed  for  a  total  of  1  5  hr  in  Oo 
followed  by  a  2-hr  anneal  in  N2  at  llOOC.  (As  determined  previouX  this  anilLuL 
sequence  assures  the  electrical  deactivation  of  any  At  impurities  prLcnt  in  the  film^.) 

The  eleetrical  characteristics  of  the  two  samples  at  room  temoeraturp  arp 

Utemture^Olif  13?  ^  Th®  ®‘"Sle-crystal  Si  obtained  from  the 

literature  (Ref  13).  The  major  differences  to  be  noted  are  those  in  the  mobilities  of 

-HOotmvT  s'e^'i^'^Air  cm2/V-sec),  and  bulk  Si 

(  1400  em  /V-sec).  Although  the  room -temperature  mobilities  of  the  two  Si/A£oO 

3lt'  wUh’f  ^  in  thTvrrraUon  oT^ 

initHMv  t«"iperature.  As  can  be  seen  from  Figure  26,  the  thicker  sample 
imtidlly  shows  an  increasing  mobility  with  decreasing  temperature  and  becomes  nearlv 

constant  from  300  to  lOOK,  but  decreases  with  decreasing  temperature  from  100  down 
o  40K.  These  data  suggest  that  the  ratio  of  the  low-temperature  to  the  room- 
in^irVh”*^^^  mobility  may  be  a  good  measure  of  the  extent  of  defect  structure.  Aecord- 

thP^Hh^r^f  I''  computed  for  77  and  300K  (temperatures  easily  obtained  in 

the  laboratory)  and  has  been  included  in  Table  XII. 

=  "f*  carrier  concentration  with  temperature  appears  to  be 

lightly  different  from  that  obtained  for  bulk  samples;  smooth  experimental  curves 

fn  oh  9pm-thick  sample  and  for  bulk  Si.  The  difference 

eters'!;sed1rfu"thG  r  manifested  in  the  different  values  for  the  param- 

?  r,  theoretically  by  the  computer  analysis  discussed  in  the 
Second  Semiannual  Report  (Ref  2);  these  parameters  are  given  in  Table  XII.  The 

found  for  the  As  donor  in  the  Si/A£20.,  samples  compared 
with  that  for  bu  k  material  is  primarily  a  result  of  the  difference  in  the  slopes  of 
the  curves  in  Figure  27  at  low  temperatures.  Another  significant  difference  between 
the  properliee  of  the  bulk  material  and  those  of  the  epitaxll  samples  ijirthe 
effective  mass,  where  there  is  an  order  of  magnitude  difference.  This  suggests  that 
the  mass  may  also  be  useful  as  a  parameter  for  eharaeterizing  film  properties. 
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Electrical  Charateristics  of  Two  As-doped  Si/(0112)Af2O 
Films  Used  for  Low-temperature  Measurements 
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Figure  26.  Electron  Mobility  as  Function  of  Temperature  for  Two  Samples  of 
Arsenic-Doped  n-Type  Si/AfgOg  Grown  under  Identical  Conditions. 
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Figure  27.  Net  Currier  Concentration  as  Function  of  Temperature  for  3.9pni-Thick 
Si/A^2^3  Film  and  for  Bulk  Single-Crystal  Si 
(Smooth  curves  drawn  through  experimental  data  points) 


It  was  originally  intended  that  these  measurements  of- earner  eoneentration 
might  yield  a  value  of  minority  earner  eoneentration  (N^  for  these  n-type  films)  that 
eould  be  related  to  defeet  strueture.  The  slight  differenee  in  this  parameter  for  the 
two  films  (thiek  and  thin)  shows  this  not  to  bo  the  case.  The  negative  result  here  is 
probably  due  to  one  of  the  following  charaeteri sties  of  the  defect  centers  involved:(l) 
they  do  not  act  as  acceptor  centers;  (2)  they  can  act  as  both  acceptor  and  donor  centers; 
(3)  they  act  as  deep  donors;  or  (4)  they  are  electrically  neutral.  Additional  measure¬ 
ments  on  p-type  samples  should  indicate  which  of  the  first  three  is  the  more  likely 
possibility.  The  fourth  alternative  is  considered  unlikely  since  defects  appear  to  be 
too  effective  as  scattering  centers.  It  should  be  noted  that  if  (2)  is  the  correct  ex¬ 
planation  then  the  model  used  in  the  analysis  of  the  data  is  not  correct,  and  the 
resulting  parameters,  although  possibly  useful,  arc  not  physically  meaningful. 

It  is  imj,iortant  to  attem^it  to  iindetatuivl  jrt'UiMjti  for  ll.-u  tllst'rcpiti'.ty  Ih  llwJ 
activation  energy  deduced  from  the  data.  At  the  present  time  it  appears  that  this 
result  may  be  a  consequence  of  inhomogeneous  strain  in  the  film.  It  is  well  known 
VJS  tiv-ir  bai^a  hrtniurtr.ii  and  Biuuiiu-siaLu  ci.cigy 

of  an  impurity.  For  large  homogeneous  strains  these  effects  tend  to  lead  to  an  increase 
in  measured  ionization  energy;  however,  if  the  film  has  a  high  locallv-inhomogeneous 
stiaui  field,  as  migHi  be  expected  with  a  large  number  of  defects,  the  situation  would 
be  different.  If,  for  example,  the  strain  "smears  out"  the  donor  energy  level  and 
"smears  out"  or  lowers  the  conduction-band  minimum,  the  energy  measured  might 
then  be  expected  to  be  lower  than  in  an  unstrained  crystal  or  in  a  homogeneously 
strained  ciystal,  and  would  be  approximately  equal  to  the  amount  the  conduction -band 
minimum  is  lowered.  Estimates  of  this  shift  have  been  given  as  -0.015  eV  (Ref  14) 
and  from  0.  040  to  0.  070  eV  (Ref  12),  either  value  of  which  could  easily  account  for 
the  difference  in  As  level  of  -0,014  eV  that  is  observed  between  the  bulk  Si  value  and 
that  for  the  thinner  film.  Further  experiments  with  films  of  different  thicknesses 
should  provide  additional  insight  into  this  problem.  The  continuation  of  these  studies 
has  been  temporarily  postponed  due  to  the  importance  of  the  anisotropy  studies  dis¬ 
cussed  below.  It  is  anticipated  that  these  studies,  including  measurements  on  p-type 
Si/Ag2f^;j.  will  be  resumed  during  the  next  reporting  period. 

^>7  hi  tlVv.  Fi opciTics  oT  5i/  At,,t7„  Turns 
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A  detailed  examination  of  the  variation  of  the  electrical  properties  in  the  plane 
of  Si  heteroepitaxial  films  has  been  undertaken  during  the  last  report  periorl,  and  has 
ibceivcd  the  major  portion  of  the  effort  aevoteci  to  Sublask  0.  The  results  are  potentially 
very  significant  to  the  preparation  of  heteroepitaxial  semiconductor  films  for  device 
usn*  ItilMllMir  sTUll))  crihfW  im  bv-Wlns  tA  imRftna  m 

work  involved  and  the  importance  of  these  results  to  the  entire  program. 

The  study  to  date  has  basically  been  in  two  areas:  first,  a  study  of  anisotropy  in 
a  relatively  low-symmetry  Si  plane,  represented  by  the  (221)Si/(1122)A^2<^3  system; 
and  second,  a  similar  study  of  a  high-symmetry  plane,  (100)Si/(0112)A^2^\3*  'He  two 
Si  planes  are  basically  different  in  the  sense  that  the  anisotropy  in  the  (22l')  plane  can 
be  expected  for  any  Si  heteroepitaxial  system,  while  the  anisotropy  seen  in  the  (lOO)Si 
plane  is  a  result  of  the  anisotropic  thermal  expansion  of  A^2t^3  J>nd  would  not  in  general 
be  present  for  (lOO)Si  growth  on  other  substrates  -  MgA^O^,  for  example.  The  results 
obtained  in  these  two  areas  are  described  in  the  following  paragraphs. 
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The  data  reported  below  are  from  measurements  taken  on  Si  grown  on  an  Ai^O^ 
orientation  ~3  deg  from  the  (1122)  plane.  Previous  orientation  studies  have  shown  that 
(221)Si//(ll22)A403  and  the  [IlOlSiMnooJAfgOg  in  the  plane  of  the  film.  For  this 
experiment,  As-doped  films  having  thicknesses  from  1.9  to  2.5  pm  and  net  donor 
coneentrations  of  2-7x1016  em  *  were  grown  at  llOOC  at  a  growth  rate  of  approximately 
2pm/min.  After  growth,  the  samples  were  annealed  in  O.,  for  30  min  followed  by  a  N 
anneal  for  2  hr  in  order  to  stabilize  film  properties  and  electrically  neutralize  any  Af** 
impurities  present  in  the  film. 


Electrical  measurements  were  earned  out  on  a  specially  designed  double  Hall- 
bridge  pattern  etehed  in  the  Si  film  using  photolithographic  techniques.  This  bridge 
pattern  is  shown  in  Figure  28.  The  legs  in  eaeh  bridge  are  spaeed  72  deg  apart  and 
the  two  bridges  are  rotated  with  respect  to  each  other  by  18  deg.  The  bridge  thus 
allows  an  independent  measurement  of  Hall  mobility  (p„)  every  18  deg  in  the  plaue  of 
the  film.  “ 

Plots  of  Hall  mobility  (p)  vs  angle  (0)  measured  from  the  [lT0]Si  direetion  are 
shown  for  representative  samples  in  Figure  29.  For  eaeh  sample,  a  least-squares 
fit  IS  obtained  from  a  eomputer  analysis  by  fitting  the  data  to  the  equation  (Pi,)”l  = 
a^beos20  (see  below).  An  anisotropy  factor  A  can  be  defined  according  to 

^  ~  ^^max  ~  ^min^^A^  ^  ~~~a  ’  where  and  are  the  maximum  and 

minimum  mobilities  determined  from  the  least-squares  fit,  and  pa  =  1/2  (p  _  +  p  •  ). 
Table  XIII  shows  the  values  of  these  parameters  for  a  number  of  samples,  observed 
anisotropy  is  seen  to  be  appreciable,  as  evidenced  by  an  average  room  temperature 
anisotropy  factor  A  of  approximately  39  percent. 

In  order  to  explain  this  behavior  the  reduction  in  mobility  due  to  the  piezoresistance 
effect  resulting  from  the  difference  in  thermal  expansion  coefficients  of  Si  and  A£_0„ 
has  been  calculated.  For  the  sake  of  completeness  the  difference  in  thermal  ^  ^ 
expansion  coefficients  for  A^^Og  parallel  and  perpendicular  to  the  (OOOlj  axis  has  also 
been  taken  into  account.  Using  representative  values  for  the  piezoresistance 
coefficients  =  -102,  =  54^  =  _i4^  x  10"^2  em2/dyne,  values  of  elastic 

constants  Cji  =  16.3,  C^g  “  6.4,  =  7.93,  x  10  dynes/cm2,  thermal  expansion 

coefficients  ogj  =  3. 9  x  10-6/deg  C,  2O3  =  9.  30  x  10-6 /deg  C,  ^  (Le)  = 

8.31  X  10  Vdeg  C,  and  a  temperature  difference  of  llOOC  over  which  the  thermal  con¬ 
traction  takes  place,  it  is  found  that  (p/tio)"^  -  1.06375-0.21891  cos  20,  where  p  is 
the  mobility  in  the  absence  of  stress  and  0  is  defined  as  above.  The  theory  predicts 
an  anisotropy  factor  of  41  percent,  in  good  agreement  with  the  experimental  values 
given  in  Table  XHI. 

An  interesting  feature  of  the  piezore  si  stance  in  this  plane  is  that  the  mobility 
is  predicted  to  be  larger  than  the  unstressed  mobility  for  9  deg  <  0  <  37  deg  (if  only 
the  first  quadrant  is  considered).  At  37  deg,  theory  predicts  that  the  mobility  re¬ 
duction  due  to  thermal  stress  is  zero,  so  that  p(37  deg)  p  .  At  this  angle  the  ex¬ 
perimental  value  of  mobility  can  be  used  to  estimate  the  mobility  degradation  due  to 
causes  other  than  thermally  induced  stress. 
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Figure  28.  Double  Hall- Bridge  Pattern  for  Studying  Anisotropic  Electrical  I^rameters 


TABLE  XIII 

Anisotropy  Parameters  of  Several  (221)Si/(1122)  AfgOg  Films 


Temp 

Sample  No.  ("K) 


-'max  ;^min 

(cm2/v-sec)  (cm^/V-sec)  (cm^/V-sec) 


RMS  Error 

{%) 


108.6 


i 
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e (Degrees) 


Figure  29.  Variation  in  Hall  Mobility  as  Function  of  Direction  in  (221)Si/(1122)  A£  O 
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If  it  is  assumed  that  the  film  mobility  at  37  deg  (i.e.,  =  (500  cm  /V-sec  for  the 

films  used  in  this  study)  is  given  by  (I/^q)  =  (l/|J^n)  +  (1/l^D),  where  Mg  is  the  bulk 
mobility  (=1000  em^/V-see),  then  =1500  cni2/v-see  represents  the  limiting  value 
of  mobility  due  to  all  mechanisms  other  than  thermally  induced  stress.  However,  pg 
probably  is  related  almost  entirely  to  defect  structure  in  the  film.  The  fact  that  the*^ 
value  of  mobility  |Xq  is  much  smaller  than  the  bulk  mobility  strongly  suggests  that 
thermally  induced  stress  is  not  the  dominant  mechanism  in  lowering  the  mobility. 

The  piezoresistance  theory  also  predicts,  in  the  absence  of  any  extexmally  applied 
magnetic  field,  an  electric  field  in  the  plane  of  the  film  but  orthogonal  to  the  current 
direction.  This  transverse  field  Eq.  is  given  theoretically  by  EqVjpQ  0.21227  sin  20, 
where  j  is  the  current  density  and  Pq  is  the  resistivity  in  the  absenee  of  stress. 
Aecordingly,  this  transverse  field  has  been  measured  (aeross  the  same  terminals  used 
to  measure  the  Hall  voltage);  the  quantity  EqVjpQ  is  shown  plotted  for  two  samples  in 
Figure  30.  The  quantity  p  is  found  from  Pq  -  ^H  the  Hall  constant 

and  Ho  as  before  is  the  mobility  at  0-  37  deg.  A  least-squares  fit  of  the  data  to  the 
predicted  equation  yields  an  average  value  for  the  coefficient  of  sin  20  of  0.2,  in  good 
agreement  with  theory. 

The  temperature  dependence  of  the  anisotropy  has  not  been  examined  in  detail, 
although  measurements  on  selected  samples  at  77K  show  a  much  larger  anisotropy,  as 
evidenced  by  the  data  in  Table  XIH.  Such  an  increase  is  consistent  with  what  might  be 
expected  for  the  variation  of  piezoresi stance  with  temperature.  Although  detailed 
measurements  of  all  of  the  piezoresistance  coefficients  are  not  available  in  the 
literature,  ttjj  appears  to  be  the  order  of  three  times  greater  at  77K  than  at  300K. 

If  it  is  assumed  that  all  of  the  coefficients  vary  by  the  same  amount,  it  can  bo  shown 
that  the  anisotropy  factor  A  would  be  expected  to  be  the  order  of  2. 7  times  greater  at 
77K.  The  data  on  two  samples  show  an  average  increase  of  a  factor  of  2.4,  consistent 
with  that  predicted. 

The  enhanced  mobility  predicted  for  some  current  directions  suggests  that  the 
(221)  plane  of  Si/A^20g  may  yield  superior  electron  mobilities.  In  comparison  with 
-ind  (lll)Si  grt-wd  uj  U/Ih  lrxboruit>jy ,  thi'  pKnii  Iad  likiocii  hlf^K>r 

mobilities.  It  would  be  expected  that  some  types  of  devices  would  show  improved 
performance  if  oriented  in  such  a  way  as  to  take  advantage  of  this  improved  mobility. 

sitvss  hit  Jm;  Ctrrrenl  .,Tncl;t;  m  ib 

the  (221)  plane,  insight  may  be  gained  into  the  scattering  mechanisms  in  Si/AP20q  by 
judicious  studies  of  the  mobility  versus  temperature.  For  example,  the  comparison 
of  such  data  at  various  angles  may  make  u  possible  lo  separate  tne  conlributions  of 
various  scattering  mechimisms  to  reduction  in  the  mobility. 

The  films  used  for  the  study  of  the  (001)Si/0li2)Af2Og  system  were  grown  in  two 
(iifferent  reactor  systems.  One  of  these  is  the  vertical  sysiem  normally  used  in  this 
program,  and  the  other  is  a  horizontai  system  of  a  type  used  for  commercial  prepara¬ 
tion  of  epitaxial  Si.  As  in  the  study  described  above,  the  films  were  formed  from  the 
thermal  decomposition  of  Sill^,  with  Ho  used  as  a  carrier  gas.  Intentional  doping  to 
concentrations  from  1  to  G  x  10l6  cm-^  was  achieved  with  the  use  of  AsHg  as  the 
.lopanl  gbs.  tetiiperaiwrfcs  -Hinged  titnn  in  Wje  biuw.uinai  oysicm  lu 

in  the  vertical  system,  and  film  thicknesses  varied  from  1.5  to  1.  8pm. 
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Electrical  data  were  taken  on  eight  samples  of  (001)Si/(01l2)A22O3.  For  every 
sample,  the  orientation  of  each  leg  of  the  Hall  pattern  was  determined  with  respect  to 
the  <100>Si  direction  which  is  parallel  to  the  (2110]A^2O3  direction  in  the  plane  of  the 
substrate.  The  mobility  in  each  leg  could  then  be  plotted  vs  the  angle  0  between  the 
current  direction  and  the  <100>  Si  axis.  Typical  mobility  data  are  shown  plotted  in 
Figure  31  for  two  samples.  The  mobility  can  be  seen  to  vary  with  angle  and  appears 
to  approach  a  maximum  at  0  deg. 

As  in  the  case  of  the  (221)Si  plane  discussed  above,  theoretical  treatment  of 
piezoresistance  predicts  a  mobility  variation  of  the  form  p/tJ^  =  (a+bcos  20)"^,  consist¬ 
ent  with  the  qualitative  behavior  observed.  Accordingly,  a  least-squares  fit  of  the  data 
to  the  above  equation  is  sho’wi  by  the  solid  curve  in  the  figure.  The  calculated  maximum 
and  minimum  mobilities  along  with  the  anisotropy  parameter  A  are  tabulated  in 
Table  XIV. 

It  can  be  seen  that  the  anisotropy  varies  between  approximately  7  and  11  percent; 
the  average  value  is  9.  2  percent.  Also  shown  is  the  RMS  error  (in  percent)  between 
the  experimental  points  and  theoretical  curve  for  each  sample. 

Measurements  were  taken  at  77K  on  several  selected  samples,  and  representa¬ 
tive  data  for  one  sample  arc  shown  in  Figure  32;  the  parameters  deduced  from  the 
curve  fitting  procedure  are  also  shown  in  Table  XIV.  The  anisotropy  is  found  to 
increase  in  going  from  room  temperature  to  77K  by  roughly  a  factor  of  three  for  those 
samples  measured.  This  is  consistent  with  the  variation  of  the  peizoresistance 
coefficients  in  Si  with  temperature. 

The  theoretical  piezorcsistance  in  the  (OOl)Si  plane  is  found  to  be 

^  =  K  +  K'  Jb^g^ogxc  “  "A^gOgZ/c] 

where  K  and  K'  are  constants.  The  anisotropy  in  this  plane  is  thus  proportional  to  the 
difference  in  thermal  expansion  coefficients  of  A£20q  parallel  and  perpendicular  to  the 
c-axis.  Such  effects  would  therefore  be  expected  to  be  absent  for  (OOl)Si  on  MgA^oO.. 
From  the  above  equation  it  follows  that,  numerically,  ^ 

=  (1.44192  -  0.04449  cos  20)"^ 

For  each  sample  measured,  has  been  calculated  from  this  equation,  and  these 
values  are  also  given  in  Table  XIV.  Comparison  with  Table  XIII  shows  that  the  values 
of  (jLq  are  very  similar  for  both  the  (001)  and  (221)Si  planes. 

The  theory  also  predicts  a  transverse  electric  field  effect  in  the  (OOl)Si  plane, 
given  by  Ej/jp^  =  0. 04449  sin  20.  Attempts  have  been  made  to  measure  this  voltage; 
however,  the  small  voltage  involved  makes  the  determination  extremely  difficult.  Repre¬ 
sentative  data  on  one  sample  are  shown  in  Figure  33.  A  least-squares  fit  to  the  data 
gives  a  coefficient  of  =0.  08,  which  is  somewhat  larger  than  predicted;  however,  the  fit 
is  not  good.  The  value  of  Pq  used  in  this  expression  is  found  from  the  value  of  Pq  deter¬ 
mined  from  the  longitudinal  mobility  measurements  described  above.  This  transverse 
voltage  also  increases  at  77K  and  is  also  shown  in  Figure  33.  The  magnitude  of  the 
increase  is  a  factor  of  between  2  and  3,  similar  to  that  for  the  longitudinal  effect. 

Note  that  the  fit  to  theory  is  much  better  at  77K,  due  to  the  larger  signal. 
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Figure  31.  Variation  of  Hall  Mobility  with  Direction  in  (001)Si/(01l2)Aje  O 
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TABLE  XIV 

Anisotropy  Data  for  (001)Si/(01]2)Ai?2n.j 


Sample 

u 

^max 

min 

A  (', ) 

BMS  Error  Cfi) 

o 

Boom  Temperatun'  Data 

102 

44(» 

40H 

427 

0.7 

1.0 

015 

lH-1 

432 

397 

415 

H.  0 

1.4 

598 

If) 

472 

425 

449 

10.5 

2.5 

048 

3412-2 

313 

472 

493 

8.2 

1.9 

711 

3411-1 

510 

47H 

497 

7.0 

2.  1 

717 

lH-2 

422 

3H7 

405 

8.0 

2.5 

584 

31 

451 

401 

420 

11.0 

2.2 

014 

13 

443 

394 

419 

11.7 

2.9 

004 

(Average) 

(402) 

(420) 

(441) 

(9.2) 

(030) 

Li(|ui(l-N^  Teniporalure  Data 


18-1 

1014 

797 

900 

24 

5.7 

1305 

31 

1003 

832 

948 

24 

7.0 

1307 

MOBILITr  fcni^/V-secl 
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c.  Electrical  Properties  of  Si  Grown  on  MgAfi204 

Additional  studies  of  the  Si/MgA<204  system  have  recently  been  undertaken,  with 
emphasis  on  the  growth  of  (lll)Si/(lll)MgA/204  in  a  mixed  H2-He  atmosphere  (Sub¬ 
task  2).  Preliminary  data  indicate  that  good  quality  n-type  films  can  be  grown  on 
(lll)MgA<’204  with  carrier  mobilities  higher  than  those  achieved  in  the  Si/AfoOo 
system.  Mobilities  as  high  as  925  cm^/V-sec  have  been  measured  in  ~2|^m-thick  films 
which  approaches  bulk  Si  mobilities.  The  variation  of  electrical  properties  with  growth 
conditions  for  these  films  is  discussed  in  the  following  paragraphs. 

Since  growth  temperature  has  been  shown  to  be  an  important  parameter  in  deter¬ 
mining  film  quality,  the  electrical  properties  of  a  group  of  (lll)Si/MgAf  2O4  samples 
having  net  donor  concentrations  in  the  10l6cm“3  range  were  examined  as  a  function  of 
growth  temperature  between  950  and  1050  C.  Low  growth  rates  (<l|jLm/min)  were  used 
for  this  study,  in  a  mixed  0.8  percent  H2-Hc  atmosphere.  The  data  are  shown  in 
Table  XV,  and  indicate  that  equivalent  quality  films  can  be  grown  over  a  wide  tempera¬ 
ture  range.  No  optimum  growth  temperature  is  suggested  by  the  data;  however,  the 
use  of  more  desirable  growth  rates  may  lead  to  a  more  definitive  answer  to  this  question. 

The  variation  in  electron  mobility  with  growth  rate  for  a  growth  temperature  of 
1025  C  is  shown  in  Figure  34.  Growth  rates  were  controlled  bv  changing  SiH4  flows 
and  also  by  varying  the  concentration  (see  Subtask  2)  from  approximately  0.8  percent 
to  4  percent.  It  is  apparent  from  this  data  that  excellent  quality  Si  can  be  grown  at  the 
higher  growth  rates  (>2H.m/min). 

Additional  mobility  data  are  shown  in  Figure  35,  which  illustrates  a  limited 
amount  of  data  at  higher  growth  rates.  This  figure  also  demonstrates  the  effect  of 
different  tanks  of  SiH4  on  filni  quality;  the  open  and  solid  circles  are  data  obtained  from 
samples  grown  with  different  sources  of  SiH4.  The  purity  of  the  SiH4  was  evaluated  by 
growing  a  thick  (~100  pm)  film  of  Si  on  p-type  Si  substrates  and  measuring  the  film 
characteristics.  The  open  circles  were  grown  with  a  tank  of  SiH4  which  yielded  films 
having  the  following  characteristics:  P  =  25  ohm-cm,  n  =  2  x  IQl^cm-S,  and  p  =  1300 
cm^/V-sec.  The  closed  circles  were  grown  with  SiH4  which  produced  films  with 
P  =  47  ohm-cm,  n  =  1.  0  x  10l4  cm-3,  and  p  =  1300  cmVv-sec.  The  purity  of  the  two 
tanks  thus  appears  to  be  equivalent  based  on  Si  growth  on  Si,  but  the  quality  of  Si  films 
produced  on  MgA<’204  is  markedly  different. 

An  Af203  substrate  was  included  in  several  of  the  Si/MgAf204  growth  runs  in 
order  to  compare  the  quality  of  the  resulting  films.  However,  it  should  be  kept  in 
mind  that  (lll)Si  may  tend  to  have  a  higher  mobility  since  stress  in  the  film  produces 
a  less  deleterious  effect  on  the  mobility  for  (lll)Si  growth.  This  data  is  recorded  in 
Table  XVI  for  two  growth  temperatures.  At  1025  C  the  growth  on  (lll)MgAf204  was 
consistently  superior  to  that  on  (0ll2)A i'203*  however,  the  lower  mobilities  on  rngOo 
are  consistent  with  those  found  previously  in  the  H2  system  at  1025C.  At  1050C,  where 

the  growth  on  (0112)A/2O„  improves  in  the  H„  system,  the  films  are  comparable  in 
quality. 
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TABLE  XV 

Variation  of  Electrical  Properties  with  Temperature 
for  (lll)Si/(lll)MgA£204  for  Low  Growth  Rates 


Growth 

Temperature 

(C) 

Substrate 

Growth  Rate 
(pm/min) 

Thickness 

(fim) 

Net  Donor 
Carrier 
Concentration 
(cm"^) 

Mobility 

(cm^/V-sec) 

950 

LO 

0.23 

2.3 

2.2  X  lO^'^ 

490 

955 

LO 

0.67 

2.0 

2.0  X  10^^ 

510 

LO 

0.44 

2.2 

1.4  X  10l7 

550 

LOT 

0.42 

2.1 

9.6  X  10l6 

460 

975 

'  LO 

0.33 

2.3 

7.3  X  10l6 

640 

LOT 

0.34 

2.4 

7.3  X  10l« 

660 

LOT 

0.33 

2.3 

7.0  X  10lf> 

520 

LOT 

0.36 

2.5 

5.5  X  10l6 

590 

1000 

LO 

0.29 

2.0 

2.  7  X  10l6 

480 

LOT 

0.31 

2.2 

1.4  X  10l6 

390 

1025 

LO 

0.23 

1.6 

3.4  X  lOlC 

560 

LOT 

0.21 

1.5 

3.2  X  10^6 

420 

LOT 

0.23 

1.6 

2.9  X  10l6 

300 

LO 

0.33 

2.3 

5.  7  X  10l(> 

660 

LOT 

0.27 

1.9 

1.5  X  1016 

480 

1050 

LO 

0.33 

2.0 

2.3  X  10l6 

630 

*LO:  Substrate  purchased  from  Union  Carbide  Corp. 

LOT:  Substrate  supplied  by  Crystal  Technology,  Inc.,  for  evaluation;  surface 
finish  poor  on  these  substrates. 
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TABLE  XVI 

Comparison  of  Properties  of  Si/AggO^  and  SiAlgAf204 
Grown  at  Two  Different  Temperatures 


Growth 

Temp 

(C) 

Percent 

Ho  in  He 

Growth 

Rate 

(pm/min) 

Thickness 

(pm) 

Resistivity 
(ohm -cm) 

Net  Donor 

Carrier 

Concentration 

(cm"^) 

Mobility 

(cm^/V-sec) 

Substrate 

1025 

0.8 

0.40 

1.95 

0.20 

7.6  X  10^® 

410 

(0n2)A#2O3 

0.41 

1.93 

0.25 

3.8  X  10^® 

660 

(lll)MgA/204 

1025 

4.2 

2.2 

1.6 

1.2 

1.3  X  10^® 

411 

(0n2)A/2Og 

2.4 

1.8 

0.77 

1.1  X  10^® 

744 

(lll)MgAf204 

4.2 

2.9 

2.2 

0.10 

1.4  X  lO^"^ 

440 

(0li2)A;2O3 

3.1 

2.3 

0.11 

9.4  X  10^® 

598 

(lll).MgAf204 

1050 

0.8 

0.27 

1.60 

0.35 

2.7  X  10^® 

650 

(0li2)A/2O3 

0.33 

1.96 

0.43 

2.3  X  10^6 

630 

(lll)MgAf204 

d.  Electrical  Properties  of  Si/A£203 


Additional  experiments  to  evaluate  further  the  parametric  relationships  among 
growth  temperature,  growth  rate,  and  substrate  orientation  have  been  carried  out  for 
the  Si/Af203  system.  For  the  experiments  examining  the  effects  of  growth  rate,  the 
amount  of  dopant  gas  added  to  the  gas  stream  was  adjusted  for  each  growth  rate  so  that 
the  carrier  concentration  in  the  film  was  nearly  constant  for  the  film  orientation  being 
grown. 

The  mobilities  of  ~2  pm-thick,  n-type  (100)Si/(0112)Af2O3  films  grown  at  various 
temperatures  are  plotted  versus  growth  rate  in  Figure  36.  The  data  show  the  increasing 
mobility  with  increasing  growth  rate  observed  previously  (Third  Semiannual  Report, 

Ref  3),  although  the  maximum  mobility  values  are  not  as  high  as  previously  reported. 
The  explanation  for  these  smaller  values  is  not  known,  but  may  relate  to  the  successive 
reuse  of  the  substrates  employed  for  these  studies.  As  seen  from  Figure  36,  there 
appear  to  be  no  significant  differences  for  growth  in  the  range  1050  to  1100  C. 

The  effects  of  growth  rate  for  (lOO)Si  at  higher  and  lower  temperatures  are 
recorded  in  Table  XVII.  These  data  do  not  suggest  appreciable  differences  in  mobility 
for  the  various  growth  rates.  The  generally  lower  values  of  mobility  with  respect  to 
those  shown  in  Figure  36,  however,  indicate  that  the  growth  temperature  is  out  of  the 
optimum  range  for  growth  of  (lOO)Si,  The  fact  that  thinner  samples  were  grown  at 
1125  C,  however,  also  leads  to  lower  mobilities  in  those  instances. 

The  effect  of  growth  rate  for  (lll)Si  is  shown  in  Table  XVIII.  As  in  the  case  of 
(lOO)Si,  the  higher  growth  rates  lead  to  improved  films.  The  high  mobilities  obtained 
in  these  films  with  thicknesses  less  than  2  pm  are  evidence  of  the  desirability  of  high 
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Figure  36.  Hall  Mobilities  of  (100)Si/(0ll2)Al2O3  as  Function 
of  Growth  Rate  (Grown  at  Three  Different  Temperatures) 


growth  temperatures  for  (lll)Si  growth.  Although  the  amount  of  Ai  autodoping  from 
the  substrate  IS  severe  at  this  growth  temperature  (Ret  3),  the  A(  does  not  appear  to 
affeet  severely  the  mobility  in  annealed  n-type  films. 


TABLE  XVn 

Effeet  of  Growth  Rate  on  the  Eleetrieal  Properties  of 
(100)Si/(0li2)A£2O3  Films  Grown  at  1025  and  1125  C 


Growth 

Temp 

CC) 

Growth  Rate 
(fim/min) 

Film 

Thiekness 

(pm) 

Resistivity 
(ohm -cm) 

Carrier 

Concentration 

(em"^) 

Mobility 

(cm^/V-sec) 

1025 

0.36 

1.8 

0. 12 

1.5  X  10^'^ 

360 

1025 

1.7 

1.7 

0.36 

3.4  X  10l6 

450 

1025 

1.8 

1.8 

0.19 

7. 9  X  10l6 

420 

1025 

2.8 

1.9 

0. 18 

8.6  X  1016 

400 

1025 

3.9 

2.0 

0.19 

7.2  X  10l6 

460 

1025 

4.9 

2.0 

0.20 

7.  7  X  10l6 

400 

1025 

5.  6 

1.9 

0.25 

5, 9  X  10l6 

420 

1125 

0.27 

1.4 

0.33 

4.7  X  10l6 

400 

1125 

1.5 

1.5 

1.6 

9. 5  X  lOl® 

400 

1125 

1.9 

1.3 

1.5 

1.0  X  10l6 

410 

1125 

2.3 

1.2 

0.82 

1.9  X  10l6 

400 

1125 

2.9 

1.3 

1.6 

9.4  X  10l5 

420 

T.iBLE  XVni 

Eifect  of  Growth  Rate  on  the  Electrical  Properties  of 
(llllSi/A^gOg  Films  Grown  at  1125  C 


Growth  Rate 
(pm/min) 

Thiekness 

(pm) 

Resistivity 

(ohm-em) 

Carrier  Coneentration 
(em-3) 

Mobility 

(cm2/v-see) 

1.0 

0.21 

5.3  X  10^® 

560 

0.31 

1.3 

2.3  X  10^® 

450 

1.7 

1.7 

1.0 

1. 1  X  10l6 

560 

1.7 

1.7 

1.2  X  10l6 

600 

2.2 

1.5 

1.1 

9.  7  X  10l5 

610 

3.8 

1.3 

0.34 

2.5  X  10^® 

730 

4.5 

1.5 

0.  75 

1.2  X  10l6 

690 
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e.  Measurement  of  Photoelectric  Effects  in  Heteroepitaxial  Films  on  AfgOg  Substrates 

The  experiments  on  photoinjection  of  electrons  from  metal  films  into  A/oOo 
substrates  have  been  continued  at  UCLA  during  this  reporting  period.  Further 
experiments  were  carried  out  to  ensure  that  the  current  observed  in  the  photoelectric 
experiments  (Ref  3)  is  indeed  due  to  photoinjection  of  electrons  from  the  metal  films 
into  the  insulator. 

In  one  such  experiment,  six  Al  film  stripes,  each  ~300A  thick  and  1  mm  wide, 
were  depostited  on  one  side  of  the  substrate,  with  a  spacing  of  approximately  1  mm 
between  adjacent  stripes  (Figure  37).  A  Au  film  was  deposited  uniformly  on  the  other 
side  of  the  substrate.  A  sharply  focussed  light  beam  was  scanned  at  a  slow  rate  across 
the  face  of  the  substrate  having  the  semitransparent  At  stripes,  and  the  photoelectic 
current  was  recorded  by  an  X-Y  recorder  as  a  function  of  the  position  of  the  light 
beam. 

Six  peaks,  with  well-defined  minima  between,  were  observed  in  the  photoelectric 
current  for  both  polarities  of  the  bias  voltage,  asshowm  in  Figures  38  and  39.  The  peaks 
correspond  to  incidence  of  light  on  a  stripe  and  the  minima  correspond  to  incidence  of 
light  in  the  spaces  between  stripes.  Observation  of  the  six  peaks  even  when  the  M 
films  v/ere  biased  positive  clearly  demonstrates  that  (1)  electrons  were  injected  into 
Af203  from  the  Au  film,  and  (2)  such  injection  occurs  only  in  those  regions  of  Al^On 
where  an  electric  field  exists,  viz.  ,  the  regions  under  the  Af  stripes. 


Figure  37.  Configuration  of  Semitransparent  Al  Stripes  Deposited  on  Af203  Substrate 
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POSITION  OF  LIGHT  BEAM 
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Figure  38 . 


Photocurrent  as  Function  of  Position  of  Scanning  Light  Beam, 
with  At  Biased  Negative 


PHOTOCURRENT  10  AMP 
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Figure  39.  Photocurrent  as  Function  of  Position  of  Light  Beam, 

with  At  Biased  Positive 


For  a  given  bias  voltage,  as  soon  as  the  light  beam  was  directed  onto  the  metal 
electi’ode,  the  photocurrent  reached  a  peak  value  in  a  time  interval  the  order  of  a 
few  milliseconds  and  then  decayed  as  a  function  of  time,  reaching  a  steady-state  value 
after  several  minutes,  thereby  indicating  that  traps  were  being  filled.  The  higher  the 
bias  voltage,  the  larger  was  the  number  of  traps  that  were  filled. 

For  the  case  of  electron  emission  from  the  Al  electrode  to  the  A^gOg,  the  steady- 
state  current  values  were  measured  for  various  photon  energies  with  the  bias  voltage 
kept  constant.  The  square  root  of  the  photoelectric  current  was  then  plotted  as  a 
function  of  the  photon  energy  to  obtain  a  Fowler  plot,  as  shown  in  Figure  40,  Two  linear 
regions  corresponding  to  two  threshold  energy  values  of  3. 15  and  4.3  eV,  respectively, 
were  observed.  Taking  the  lower  value  as  that  corresponding  to  the  height  of  the  inter¬ 
face  barrier  between  A2  and  A^2^3»  the  energy-band  diagram  shown  in  Figure  41  was 
constructed.  The  work  function  of  A£  was  taken  to  be  4. 1  eV  (Ref  15),  and  this  gives  a 
value  of  1.0  eV  for  the  electron  affinity  of  A^g^^g*  which  agrees  w’ell  with  the  value 
quoted  by  Pong  (Ref  Ifi)  for  amorphous  A^203.  The  value  of  3. 15  eV  for  the  interface 
barrier  height  agrees  well  with  the  value  3.05  eV  measured  by  Szydlo  and  Poirier 
(Ref  17)  for  the  polycrystalline  A^203-M  interface. 

When  the  photoelectric  current  corresponding  to  emission  of  electrons  from  other 
metals,  such  as  Au,  Ag,  and  Cu,  was  measured  as  a  function  of  photon  energy,  the 
resulting  Fowler  plots  were  again  characterized  by  two  linear  regions.  The  lower 
threshold  energy  values,  again  taken  as  the  interface  barrier  heights,  when  combined 
with  the  value  of  1. 0  eV  for  the  electron  affinity  of  A^2®3  gave  work  function  values 
(Ref  18)  agreeing  with  previously  published  values  for  the  metals  involved. 

In  all  these  cases,  the  threshold  energy  value  corresponding  to  the  second 
linear  region  was  1  to  1. 1  eV  higher  than  the  interface  barrier  height.  This  suggests 
the  possibility  that  the  higher  threshold  energy  might  correspond  to  vacuum  emission 
of  electrons  from  the  metal  into  pinholes  in  the  metal  film.  Further  experiment  is 
needed  to  verify  this. 

The  fact  that  a  photoelectric  current  is  observed  in  A£203  layers  several  mils 
thick  clearly  establishes  that  the  range  of  photoinjected  electrons  is  much  larger 
than  the  value  of  130  A  reported  for  amorphous  A£203  layers  (Ref  16).  It  is  possible 
to  measure  the  range  of  the  injected  electrons  by  carrying  out  measurements  with 
A22C)3  substrates  of  several  different  thicknesses. 

The  magnitude  of  the  current  as  a  function  of  the  thickness  of  A£  films  deposited 
on  the  same  A^2^3  substrate  was  also  measured,  w'ith  the  light  incident  on  the  A^ 
films.  A  peak  was  observed  around  2500  A  thickness,  and  the  current  was  found  to  be 
an  appreciable  fraction  of  the  peak  even  when  the  thickness  of  the  A2  layer  was  as 
large  as  5000  A.  This  result  is  in  contrast  to  the  mean  range  for  electrons  in  A£  as 
measured  by  other  workers  (Ref  19). 

Using  this  technique,  it  is  possible  to  measure  the  transport  properties  of 
electrons  injected  into  ordered  insulators  and  to  determine  the  range  of  electrons 
in  metal  films.  All  the  measurements  described  above  were  conducted  at  room 
temperature;  it  is  planned  that  similar  experiments  be  carried  out  at  various  other 
temperatures. 


82 


3.0  4.0  6.0 

ENERGY  (ev) 

Figure  40.  Fowler  Plot  of  Photocurrent  au  Function  of  Photon  Energy 

for  Af-A£20g  Interface 
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Figure  41.  Energj'-Band  Diagram  for  A2-A220^  Interface 

The  effect  of  annealing  on  the  photoinjection  properties  of  metal  films  was  also 
investigated.  A  metal-insulator-metal  san^ich  ccnsisting  of  semitransparent  Au  films 
(300  to  400  A  thick)  deposited  on  the  two  opposite  faces  of  an  (0ll2)  A£203  substrate 
0.016  in.  thick  was  fabricated.  Photoelectric  measurements  were  made  on  this 
sample  under  four  different  experimental  conditions,  corresponding  to  opposite 
polarities  of  bias  voltage  first  with  light  incident  on  one  metal  film  and  then  with  light 
incident  on  the  other  metal  film.  The  samples  were  then  annealed  for  30  min  at  SOOC 
in  vacuum,  after  which  the  photoelectric  measurements  were  repeated. 

It  was  found  that  the  photoelectric  characteristics  had  changed  considerably 
as  a  result  of  the  annealing.  For  example,  the  plot  of  the  square  root  of  the  photo¬ 
current  as  a  function  of  photon  energy  app>eared  much  smoother  (i.e. ,  devoid  of  fine 
structure)  after  annealing,  for  electron  injection  from  the  illuminated  film. 

Some  experiments  were  carried  out  in  an  attempt  to  determine  the  mobility  of 
injected  electrons  in  the  insulator  by  using  a  pulsed  N2  gas  laser  for  transient  internal 
photoemission  studies.  Due  to  a  high  noise  pick-up  in  the  measuring  circuit  from  the 
switching  circuit  of  the  laser  the  signal  could  not  be  observed.  The  experiments  have 
not  yet  been  repeated. 
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f.  Measurement  of  High-Field  Transport  Properties  of  Heteroepitaxial  Films 

Experiments  to  determine  the  high-field  transport  properties  of  heteroepitaxial 
semiconductor  films  on  insulating  substrates  have  been  continued  at  UCLA  during 
this  reporting  period.  The  experimental  arrangement  was  shown  in  Figure  27  of  the 
Third  Semiannual  Report  (Ref  3).  Measurements  made  in  the  first  half  of  the  second 
year  had  not  yet  revealed  the  presence  of  high-field  effects. 

Considerable  difficulty  was  experienced  initially  in  fabricating  the  device 
structures  required  for  these  measurements.  However,  a  few  test  devices  were 
successfully  made  and  were  tested.  It  was  found  that  the  metal  pads  deposited  on  the 
semiconductor  films  did  not  provide  a  good  ohmic  contact.  Finally,  after  considerable 
additional  experimentation,  good  ohmic  contacts  were  obtained  on  GaAs  films  by 
evaporating  an  alloy  of  Ag,  In  and  Ge  (90  percent  -  5  percent  -  5  percent)  onto  the  GaAs 
film  and  heating  the  sample  for  t^\o  minutes  at  70dC  in  He  gas  ambient. 

Some  resultsofthe  high-field  transport  measurements  are  given  in  Figure  42, 

It  can  be  seen  that  the  current  deviates  from  linear  behavior  for  electric  fields  beyond 
-1.4  kV/cm  and  exhibits  a  peak  around  3  kV/cm.  However,  because  of  heating  effects 
(a  low- resistivity  sample  was  used),  the  current  tends  to  increase  again  at  fields 
higher  than  3.4  kV/cm.  In  addition,  significant  sparking  was  observed  at  field 
strengths  of  6  to  8  kV/cm. 
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Figure  42.  High- Field  Transport  Effects  in  Heteroepitaxial  GaAs/AfgOg 
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In  the  case  of  Si  films  on  A£203,  electrodes  of  Au  doped  with  1  percent  Sb  were  deposited 
on  the  Si  films  and  heated  to  700C  for  two  minutes  in  a  No  ambient.  Ohmic  contacts 
were  still  not  obtained,  however,  and  this  prevented  the  observation  of  high-field 
characteristics.  Further  experimentation  is  progressing  directed  toward  achieving 
good  ohmic  contact.  It  was  noted  that  the  nonohmic  contacts  did  becjome  ohmic  after 
a  few  occurrences  of  the  sparking  referred  to  above.  If  necessary,  such  a  "forming^ 
process  will  be  systematically  employed  for  preparing  the  required  ohmic  contacts. 

7.  SUBTASK  7.  DESIGN  AND  FABRICATION  OF  DEVICES 

During  the  first  year  apparatus  for  determining  carrier  lifetime  by  C-V 
measurement  in  MOS  structures  was  designed,  constructed  and  tested.  A  special 
MOS  structure  was  designed  for  measurement  of  channel  conductance,  high-  and 
low-field  transport  properties,  and  various  interface  characteristics  of  heteroepitaxial 
films.  Initial  attempts  to  fabricate  Schottky-barrier  diodes  in  Si/A£203  films  as  a 
means  of  evaluating  their  electrical  properties  were  not  successful. 

In  the  first  six  months  of  the  second  year  preliminary  design  of  a  Schottky-barrier 
type  of  FET  was  completed  for  use  in  making  experimental  FET  structures  in  GaAs/. 
insulator  films  for  operation  at  1  GHz.  Preliminary  results  on  carrier  lifetime  in 
Si/A^203  films  were  obtained,  after  initial  development  of  the  measurement  technique 
on  bulk  single-crystal  Si  samples  was  completed  and  after  impurity  contamination 
problems  encountered  in  the  oxide  growth  process  were  solved. 


The  main  device-oriented  effort  has  centered  about  the  carrier  lifetime  studies 
using  the  MOS  pulsed-capacitance  technique  and  the  attempts  to  fabricate  a  Schottky- 
barrier  FET  in  GaAs/A2203.  Both  activities  are  being  carried  out  at  UCLA. 

The  GaAs  FET  work  was  hampered  seriously  earlier  in  this  reporting  period 
by  difficulties  encountered  in  processing  the  heteroepitaxial  samples  through  the 
photolithography  steps  using  the  new  device  masks,  fabrication  of  which  was  completed 
shortly  after  preparation  of  the  Third  Semiannual  Report.  It  is  believed  that  these 
processing  problems  are  now  successfully  resolved,  after  several  sessions  in  which 
the  procedures  being  used  at  UCLA  were  critically  reviewed  and  empirically  modified 
by  NR  personnel  in  NR  laboratories.  The  device  fabrication  work  has  been  resumed, 
although  some  further  delay  is  being  encountered  while  heteroepitaxial  samples  are 
being  prepared  at  NR  as  part  of  the  Subtask  2  effort. 

Lifetime  studies  using  the  technique  of  pulsed  MOS  capacitance  have  been 
continued  during  this  report  period.  Many  of  the  Si/A^203  samples  received  earlier 
were  found  to  have  a  large  number  of  pinholes,  and  this  gave  rise  to  device  fabrication 
difficulties.  The  samples  prepared  more  recently  did  not  have  many  pinholes, 
containing  typically  about  20  pinholes  of  diameter  larger  than  ~1  fJLm  on  each  wafer. 

The  number  of  pinholes  was  found  to  be  larger  in  the  (111)  Si  orientation  than  in  the 
(100)  orientation.  The  lower-resistivity  samples  had  fewer  pinholes  and  the  number  of 
pinholes  seemed  to  progressively  decrease  as  the  thickness  of  the  Si  film  was 
increased  -  not  a  surprising  effect. 
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In  the  oxidation  process  used  for  making  the  MOS  structure,  the  impurity 
distribution  is  changed  in  the  bulk  of  the  semiconductor  in  the  region  near  the  oxide 
interface.  An  analysis  was  made  to  consider  the  effect  of  impurity  redistribution  on 
the  interpretation  of  lifetime  measurement  data  using  the  pulsed  MOS  capacitance 
technique,  and  this  is  described  in  subsequent  paragraphs. 


The  use  of  the  pulsed  MOS  capacitance  technique  to  determine  the  bulk  generation 
lifetime  Tg  and  the  surface  recombination  velocity  s  of  both  films  and  wafers  of  Si 
as  well  as  GaAs  has  been  reported  by  several  authors  (Refs  20  and  21).  In  this 
technique  the  capacitance  C  of  the  MOS  structure  is  monitored  as  a  function  of  lime 
after  the  application  of  a  depleting  voltage  step.  Zerbst  (Ref  20)  has  shown  that  for  an 
ideal  MOS  structure  having  an  oxide  thickness  Xq,  corresponding  to  an  oxide  capacitance 


C  ,  a  plot  of 


ields  a  straight  line  whose  slope  is  related  to 


and  whose  intercept  on  the  - 
capacitance  at  time  1  co. 


dt 


axis  is  related  to  s  (Ref  3). 


is  the  MOS 


This  technique  has  been  suggested  for  in-process  measurement  of  Xg  and  s  during 
device  fabrication,  and  has  recently  been  extended  (Ref  22)  for  the  case  of  schottky- 
barrier  diodes  to  determine  the  lifetime,  energy  level,  and  density  of  both  deep  and 
shallow  impurity  levels.  In  a  study  of  the  impurity  redistribution  in  Si  near  the  SiOn- 
Si  interface  that  occurs  during  thermal  oxidation.  Deal  et  ^  (Ref  23)  used  an  "effective" 
surface  doping  concentration  in  examining  the  C-V  characteristic  of  an  MOS  capacitor. 
However,  this  approach  is  unsuitable  in  pulsed  measurements  to  determine  lifetime 
since  the  depletion  layer  width  W,  and  hence  capacitance,  changes  with  time  t. 


It  is  possible  to  illustrate  quantitatively  the  effect  of  impurity  redistribution  on 
the  interpretation  of  bulk  generation  lifetime  using  the  pulsed  MOS  capacitance 
technique,  and  this  is  done  as  follows. 


Grove  ct  ^  (Refs  23  and  24)  have  showai  that  if  a  wafer  of  Si,  initially  uniformly 
doped  with  boron  to  a  concentration  Np,  is  sibjected  to  a  thermal  oxidation  process,  the 
resulting  concentration  Ng  of  boron  in  the  Si  at  the  Si02-Si  interface  will  be  lower  than 
the  bulk  value  by  a  factor  typically  from  2  to  7,  depending  on  oxidation  conditions.  For 
a  phosphorus-doped  Si  wafer  Ng  will  be  larger  than  Ng  by  a  factor  of  from  1. 1  to  4. 

An  exponential  doping  profile  of  the  form 


N(x)  =  Ng  +  (Ng 


-  V 


-x/x, 


(4) 


was  compared  with  the  exact  solution  of  Grove  et  al  and  was  found  to  be  in  very  close 
agreement  (i  5  percent).  The  parameter  Xj  is  used  to  fit  the  exponential  doping 
profile  to  the  exact  profile  at  Grove  et  ai  (Ref  24)  and  is  plotted  in  Figure  43  against 


N 

N 


s 

B 


for  both  boron  and  phosphorus  under  various  oxidation  conditions. 


Using  the 


exponential  doping  profile  in  a  double  integration  of  Poisson’s  equation,  an  expression 
for  the  voltage  across  the  space  charge  region  is  obtained. 
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Figure  43.  vs  x.  (in  fim)  for  Phosphorus- and  Boron-Doped  Si  Under  Various 

Nb  1 

Oxidation  Conditions.  (Compiled  with  the  aid  of  data  from  Ref  24. ) 
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Following  the  analysis  of  Zerbst  |Ref  20)  it  can  be  shown  that  the  rate  of  increase 
of  charge  Qg  at  the  Si02-Si  iiterface  is  rekted  to  the  MOS  capacitance  C{t)  by 


where  tg  is  the  permittivity  of  the  Si  and  q  is  the  electi-onic  charge.  It  has  been 
assumed  that  the  area  A  of  the  metal  gate  electrode  is  large  enough  to  neglect  surface 
effects  and  that  the  total  capacitance  may  be  represented  by  C  in  series  with  the 

^s  A 

depletion  layer  capacitance  .  'fhe  first  term  on  the  right-hand  side  of 

Eq  5  is  the  result  due  to  Zerbst  assuming  a  uniformly-doped  sample,  and  the  second 
term  is  the  modification  due  to  the  exponential  doping  profile. 


Setting  Eq  5  equal  to  the  bulk  and  surface  generation  current  terms  gives 


Nt,  « 

B  s 


2C 


(«) 


where  ijj  is  the  intrinsic  carrier  concentration  in  the  semiconductor  and  the  function 
B(C)  is  given  by 


B(C)  = 


(7) 


w1>  ’re  Cl  -  . 

*  /C  \ 

A  plot  of  the  left  side  of  Eq  6  as  ordinate  versus  l-^  "  ^  /  3  straight 

line  whose  slope  is  related  to  and  whose  intercept  on  the  vertical  axis  is  reicted 
to  s.  A  Zerbst  plot  assuming  umform  doping  yields  a  lifetime  value  ^  due 

to  the  function  B(C).  If  \V  <  xj  then  the  entire  space-charge  region  is  located  within 
a  region  of  impurity  redistribution,  and  B(C)  departs  from  unity. 


Figure  44  showc  the  v-"riation  of  B(C)  during  the  capacitance  transient  for  various 
oxidation  conditions.  From  Eq  G  and  7,  and  Figures  43  and  44  n  can  be  seen  that  for 
boron  the  function  B{C)  <  1  and  ^  >  t  ».  For  phosphorous  the  function  B{C)  >  1  and 

^  _  f  6  6 
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Figure  44.  Impurity  Redistribution  Function  B(C)  vs  Normalized  MOS  Capacitance 
C/C  for  Phosphorus-  and  Boron-Doped  Si  under  Various  Oxidation  Conditions 


For  W  »xi  then  B(C)  =  1.0  and  Xg*  =  Xg  for  the  entire  transient  response. 

,«,m  doping  profile  used  in  this  analysis  illustrates  the  effect  of 

o'!  ®‘^‘bution  on  the  interpretation  of  litetime  measurements  from  pulsed 

^e?ter?r"!r  temperature  processing  such  as  annealing, 

fion  fS  ?  oxidation  step  may  cause  departures  in  the  impurity  redistribu- 

^  ^  °  apparent  changes  in  lifetime  unless  the  actual  doping 

prolile  IS  used  in  the  analysis.  ^  ” 

mpf  ^  obtained  from  the  C-V  characteristic 

til  p  v’  a  computer  calculation  is  usually  used  to  extract  the  doping  profile  from 
Irnfn"'^  characteristic,  it  is  convement  simply  to  incorporate  this  measured  doping 
profile  into  the  computer  program  analyzing  the  capacitance-time  response,  thereby 

SnuJU"®  Tf®  P*’®'-'*®®  “rements  of  Xg  and  s.  Alternatively,  thermal  oxide 
growth  conditions  may  be  ;hosen  from  Figures  43  and  44  that  result  in  negligible 

;oTT/ed!*t"l:“SrS^ 

r  lifetime  measurements  on  both  Si  and  GaAs  films  on  substrates  primarily 

of  are  continuing.  Detailed  results  as  functions  of  film  thickness,  carrier 

concentration  and  crj-stallographic  orientation  are  being  obtained  and  wall  be 
described  in  the  next  semiannual  report. 
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SECTION  III 

WORK  PLANNED  FOR  NEXT  SIX  MONTHS 


The  planned  effort  during  the  next  six  months  in  the  seven  subtasks  will  include 
^  the  following.  Deviations  from  these  activities  will  be  introduced  when  and  as  the 

i''  progress  of  the  total  program  retiuires  it. 

j  SURTASK  1.  theory  OF  EPITAXY  AND  IIETEROEPITAXIAL  INTERFACES 

The  investigation  of  surface  reconstruction  of  the  A^203  lattice  using  the 
■:  previously  determined  potentials  will  be  eontinued  and  eompleted.  An  existing  Morse 

J  1°*'  Si  will  be  investigated  for  suitability  in  the  heteroepitaxial  problem,  and 

the  adatom  potential  interaction  between  Si  and  the  will  be  parameterized, 

f^f^aiputer  simulation  of  the  growth  of  Si  on  A^203  using  previously  determined 
potentials,  surfaee  reeonstruetion,  and  parameterized  adatom  interaetions  will  then 
be  initiated;  it  is  intended  that  the  simulation  work  will  begin  late  in  the  eoniing  six— 
montli  period,  although  it  may  be  delayed  until  early  in  the  final  six  months  of  the 
program. 

The  extensive  theoreUcal  study  of  the  anisotropy  in  mobility  due  to  differential 
thermal  eontraetion  at  the  interfaee  between  film  and  substrate  will  be  eompleted  and 
a  detailed  paper  will  be  prepared  and  submitted  for  publication.  Further  work  in  this 
area  may  be  performed  as  required  to  eomplement  and  further  advanee  the  experi- 
niental  studies  of  film  mobility  and  other  eleetrical  properties. 

2.  SUBTASK  2.  DEPOSITION  STUDIES  AND  FILM  PRE PARATION 

Beeause  of  the  very  eneouraging  eleetrical  results  obtained  to  date  for  Si  grown 
on  MgA£204  substrates  in  H2-fIe  mixtures  this  system  will  be  explored  further,  not 
only  for  MgA2204  but  also  for  Aif203.  Ihis  study  will  eneomipass  various  substrate 
orientations,  and  eonditions  will  be  established  for  optimized  growth  to  permit 
identification  of  the  preferred  substrate.  Further  investigation  of  the  autodoping 
phenomenon  as  a  fmtetion  of  growth  atmosphere  (H2-He  ratio),  temperature,  substrate 
,  orientation,  and  growth  rate  will  be  earned  out  with  both  A^2^3  and  MgA2204 

;  substrates.  The  nature  of  the  autodoping  meehanism  will  be  eii^lored  using  bulk  Si 

wafers  and  A^203  and  MgA^2t>4  substrates  as  sourees  of  dopant,  in  various  gaseous 
atmospheres.  Some  further  attention  will  be  given  to  the  effeets  on  the  properties  of 
Si/MgOg  films  of  gas  flow  patterns,  reaetor  geometries,  and  temperature  and  gas 
■  flow-rate  eontrol.  Studies  of  film  nueleation  phenomena  in  the  Si/A£203  system  will 

be  resumed  and  additional  deposition  experiments  at  low  pressures  (<0. 1  torr)  will 
be  earried  out  espeeially  to  eomplement  the  in  situ  CVD  experiments  in  the  electron 
microscope  (Subtask  5).  A  p-type  doping  capability  for  Si  films  on  insulators  will  be 
developed  so  as  to  study  multilayer  film  interfaces  and  permit  fabrication  of  multi¬ 
layer  device  structures.  Methods  for  optimizing  the  controlled  doping  of  thin  epitaxial 
layers  of  GaAs  on  insulators  will  be  developed  in  order  to  prepare  eomporite  samples 
p  for  other  subtasks  of  the  program,  particularly  those  which  involve  device  fabrication 

and  evaluation. 
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3.  SUBTASK  3.  ANALYSIS  AND  PURIFICATION  OF  CVD  REACTANTS 


Analysis  of  CVD  reactants  for  impurity  content  will  be  continued,  primarily  by 
mass  spectronietric  techniques.  The  efforts  to  obtain  reactants  of  improved  purity 
will  continue,  depending  heavily  upon  cooperation  from  the  vendors  involved.  The 
previously  planned  analysis  of  the  impurity  content  of  a  scries  of  Si  and  GaAs  films 
as  a  function  of  distance  from  the  film-substrate  interface,  as  an  indirect  measure 
of  the  purity  of  the  reactants  used,  will  be  carried  out  during  this  six-month  period. 
The  investigations  of  the  basic  chemisti’y  involved  in  the  formation  of  Si  by  SiH4 
pyrolysis  and  of  GaAs  by  the  trimethylgallmm  (TMG)  -  Asllg  reaction  at  CSU  San  Diego 
will  emphasize  a  study  of  the  influenee  of  A  £2^3  surfaces  on  the  mechanism  of 
decomposition  of  SiH4  and  Sioll^j  as  a  lunction  of  temperature  and  an  examination  of 
the  decomposition  1,  odes  and  reaction  products  of  'I’MG  and  AsILj  as  a  function  of 
temperature,  concentration  ratio,  and  II2  partial  pressure  to  determine  the  require¬ 
ments  for  formation  of  optimum-quality  GaAs. 

4.  SUBTASK  4.  PREPARATION  AND  CHARACTERIZAT  ION  OF  SUBSTRATES 

Substrate  surfaces  of  both  A^203  and  MgA^204  will  be  prejiarcd  for  experimental 
use  in  the  program  using  the  best  polishing  and/or  etching  techniques  developed  to 
date.  Additional  effort  will  be  devoted  to  improving  the  surface  preparation  techniques 
for  both  materials,  with  emphasis  on  chemical-mechanical  polishing,  liquid-phase 
etching  procedures,  and  the  use  of  annealing  cycles  in  conjunction  with  liquid-  and 
gas-phase  etching  procedures,  for  both  Ae203  and  MgAi204.  The  effects  of  gas-phase 
etching  of  MgA^204  will  be  further  examined  and  correlated  with  Si  film  growth 
characteristics  and  resultant  film  properties.  Ion-beam  sputtering  will  also  be 
investigated  as  a  means  of  removing  damaged  surface  layers  and  contaminants  from 
substrates  prior  to  film  deposition.  The  experiments  with  Rutherford  backscattering 
of  charged-particle  beams,  reflection  electron  diffraction,  replica  electron  micro¬ 
scopy,  and  Auger  electron  spectroscopy  will  be  continued  in  the  attempt  to  correlate 
the  properties  of  CVD  Si  films  on  insulator  substrates  with  the  pre-deposition 
properties  of  the  substrate  surfaces. 

The  development  of  the  ion-boam  sputtering  technique  for  preparing  thin 
(~200A)  substrates  of  A£203  the  in  si^  CVD  experiments  of  Subtask  5  will  continue. 
The  emphasis  will  be  on  producing  more  uniform  thickness  of  the  Af203  in  the  thinned 
regions  to  be  used  in  the  deposition  experiments  and  on  developing  mechanical  or 
other  polishing  or  etching  procedures  which  will  have  higher  yield  for  the  preparation 
of  the  pre-ion-thinning  substrates  of  thickness  the  order  of  0.  001  in.  or  less. 

5.  SUBTASK  5.  STUDIES  OF  D4  SITU  FILM  GROWTH  IN  THE  ELECTRON 

MICROSCOPE 

The  CVD  microchamber  will  be  installed  and  tested  early  in  the  coming  six 
months;  upon  completion  of  the  focus,  resolution,  alignment,  heating  and  gas-flow 
tests,  the  initial  CVD  experiments  will  begin.  The  first  in  ^  work  will  concentrate 
on  the  temperature  and  substrate  orientation  variables  of  the  Si  CVD  process  based 
on  SiH4  pyrolysis;  the  effect  of  Sillq  gas  pressure  will  be  examined  later.  The 
nucleatlon  and  early-stage  growth  processes  will  be  studied  to  determine  the  mecha¬ 
nism  whereby  structural  defects  are  incorporated  into  the  growing  film.  Improve¬ 
ments  in  the  CVD  microchamber  and  the  motion  picture  recording  system  will  be  made 
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as  required.  It  is  anticipated  that  numerous  changes  will  be  made  in  the  experimental 
apparatus  and  the  procedures  used  as  the  work  progresses;  this  period  of  exploratory 
work  will  probably  occupy  the  entire  six-month  period. 

6.  SUBTASK  6.  EVALUATION  OF  FILM  PROPERTIES 

The  routine  evaluation  of  the  structural  and  electrical  properties  of  Si  and  GaAs 
heteroepitaxial  films  will  be  continued  by  means  of  various  methods  of  X-ray  and 
electron  diffraction  analysis,  metal lographic  analysis,  and  measurement  of  both 
fundamental  and  device-related  electrical  properties.  Rutherford  backscattering 
analysis  using  charged  particle  beams  will  be  continued  for  establishing  the  depth  and 
density  of  defects  in  the  films  and  the  results  correlated  with  characteristics  of  the 
substrates  used  to  grow  the  films. 

Studies  of  the  anisotropy  in  electrical  properties  in  the  Si/A;f203  system  will  be 
continued  with  emphasis  on  Si  films  grown  on  or  near  the  (111)  orientation.  In 
addition,  further  studies  will  be  made  of  films  having  orientations  near  the  (100)  Si 
plane.  Low -temperature  studies  of  the  electrical  properties  of  Si/A£203  will  be 
resumed  and  compared  with  similar  studies  to  be  undertaken  on  the  Si/MgA/204  system, 
using  both  n-  and  p-type  films.  Attempts  will  be  made  to  identify  experimentally  the 
stress  conditions  in  these  films  using  X-ray  techniques.  The  comparison  of  the 
properties  of  Si  films  on  Aie203  and  MgA£204  substrates  will  continue  so  that  the 
goal  of  selecting  the  preferred  substrate  species  for  optimized  properties  in  the 
Si/lnsulator  system  can  be  realized. 

Measurements  of  the  escape  length  in  Si  and  GaAs  films  are  to  be  completed. 
Studies  of  transient  internal  photoemission  and  the  charge  transport  properties  in 
^^2^3  various  temperatures  will  be  continued.  Observations  will  be  made  of  the 
effect  of  various  annealing  treatments  on  the  photoinjection  properties  measured  in 
A£203  -  based  systems.  In  addition,  the  photoelectric  current  in  these  same  configu¬ 
rations  will  be  measured  at  very  low  electric  fields.  Finally,  the  high-field  transport 
properties  in  both  Si  and  GaAs  will  be  measured  as  a  function  of  temperature,  film 
resistivity  and  thickness,  and  crystallographic  orientation. 

7.  SUBTASK  7.  DESIGN  AND  FABRICATION  OF  DEVICES 

Determination  of  carrier  lifetimes  in  Si  and  GaAs  films  on  both  AI2O2  and 
MgA^204  substrates  by  the  MOS  pulsed  C-V  techniques  will  be  continued.  Systematic 
measurements  will  be  made  of  lifetime  as  a  function  of  film  thickness,  resistivity, 
and  crystallographic  orientation;  attempts  will  also  be  made  to  observe  lifetime 
dependence  on  temperature  of  measurement.  The  variation  of  the  properties  of  simple 
junction  diodes,  MOS  structures,  and  FET's  with  conditions  of  growth  of  the  films 
will  be  used  as  a  measure  of  film  quality.  Comparison  device  structures  will  be  made 
and  evaluated  in  the  Sl/A£203  and  Sl/MgA£204  systems  to  help  select  the  preferred 
substrate  for  optimized  Si  film  properties.  The  Schottky-barrler  FET  fabrication 
experiments  will  continue  in  GaAs/A?203  and GaAs/MgA£204  structures;  following 
successful  fabrication  and  appropriate  evaluation  of  these  devices  at  -IGHz  (design 
frequency),  attempts  will  be  made  to  extend  the  operation  Into  the  microwave  region. 
Some  attention  will  also  be  given  to  the  fabrication  of  other  special  experimental 
device  structures  to  be  used  to  characterize  heteroepitaxial  film  quality  and  utilize 
the  unique  properties  of  these  films. 
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SECTION  IV 

PROGRAM  SUMMARY  TO  DATE 


The  overall  objeetive  of  the  program,  for  whieh  this  is  the  Fourth  Semiannual 
Report,  is  to  carry  out  a  fundamental  study  of  thenueleaticmand  film  growth  mechanisms 
in  heteroepitaxial  semiconductor  thin  films,  leading  to  new  knowledge  and  understanding 
of  these  processes,  and  then  to  apply  these  results  to  the  preparation  of  improved 
semiconductor  thin  films  and  thin-film  devices  on  insulating  substrates. 


The  specific  technical  objectives  of  the  three-year  program  are  the  following: 

1.  Investigation  of  the  many  aspects  of  the  mechanisms  of  heteroepitaxial  film 
growth,  to  establish  (through  accumulation  of  basic  knowledge)  sets  of  tech¬ 
nical  guidelines  for  the  preparation  of  better  films  whieh  ean  then  be  applied 
to  real  situations. 

2.  Preparation  of  improved,  high-quality,  deviee-grade  heteroepitaxial  films 
of  Si  and  GaAs  on  insulating  substrates  by  chemical  vapor  deposition  (CVD) 
methods . 

3.  Development  of  methods  of  characterizing  heteroepitaxial  films  as  to  their 
suitability  for  subsequent  deviee  fabrieation. 

4.  Design  and  fabrication  of  selected  thin-film  devices  whieh  take  advantage 
of  the  unique  properties  of  such  films. 

The  general  plan  for  accomplishing  these  objectives  involves  the  study  of  the  funda¬ 
mentals  of  heteroepitaxial  semieondu  iitor  film  growth  on  insulating  substrates,  with 
specialized  device  fabrication  used  as  a  means  of  evaluating  certain  film  properties 
(and  thus  as  a  measure  of  film  qualify  as  the  program  progresses)  and  as  a  means  of 
exploiting  those  properties  unique  to  heteroepitaxial  semiconductor-insulator  systems. 

The  determination  of  the  fundamental  mechanisms,  properties,  and  processes 
to  be  investigated  is  based  on  extensive  background  luiowledge  of  epitaxy  and  thin-film 
device  difficulties  encountered  for  many  years  in  many  laboratories.  The  problems 
under  study  are  not  restricted  to  those  identified  a  priori;  experimental  and  theoretical 
attention  is  shifted  as  needed  as  the  program  progresses. 


The  program  involves  both  theoretical  and  experimental  investigation  of  the 
nucleation  and  growth  mechanisms  of  heteroepitaxial  films  in  semieonduetor-insulator 
systems.  The  theoretical  studies  consist  of  two  types  of  activity:  (1)  direct  interaction 
with  the  experimental  program  involving  data  analyses,  suggestion  of  definitive  experi¬ 
ments,  and  postulation  of  specific  models  to  explain  experimental  observations,  and 
(2)  development  of  original  contributions  to  the  theory  of  heteroepitaxial  growth,  with 
the  goal  of  generating  significant  advances  in  fundamental  epitaxy  theory. 

The  experimental  investigations  are  also  of  two  types:  (1)  fundamental  explorations 
to  delineate  mechanisms  and  general  empirical  principles  of  the  heteroepitaxial  growth 
process,  and  (2)  practical  studies  that  accompany  the  fundamental  investigations  so  that 
any  new  developments  can  be  immediately  applied  to  the  improvement  of  semiconductor 
films  and  thin-film  devices  on  insulating  substrates. 
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The  work  has  emphasized  the  CVD  method  of  growing  semiconductor  thin  films 
because  of  its  importance  in  the  semiconductor  industry.  Specifically,  the  program 
emphasis  is  on  films  of  Si  and  GaAs  and  substrates  of  sapphire  (Af203)  and  spinel 
(MgAf204)*  The  initial  emphasis  has  been  on  the  Si-on-Af203  system,  with  increasing 
attention  being  given  to  the  Si-on-MgA#204  and  GaAs-on-A£203  systems.  Si  and  GaAs 
have  been  chosen  because  of  the  preeminence  of  the  former  in  the  semiconductor 
industry  and  the  high-frequency  and  high-temperature  attributes  of  the  latter. 

The  program  is  carried  on  primarily  at  the  facilities  of  the  North  American 
Rockwell  (NR)  Electronics  Group  in  Anaheim.  Parts  of  three  of  the  subtasks  have 
been  performed  by  personnel  of  the  University  of  California  at  Los  Angeles  (UCLA), 
in  the  Department  of  Electrical  Sciences  and  Engineering  and  in  the  Chemistry  Depart¬ 
ment,  and  by  personnel  of  the  California  State  University,  Sau  Diego  (CSUSD),  by 
means  of  subcontracts  from  NR. 

The  accomplishments  of  the  contract  program  to  date  are  summarized 
subtask  as  follows. 

1.  SUBTASK  1.  THEORY  OF  EPITAXY  AND  HETEROEPITAXIAL  INTERFACES 

During  the  first  year  a  formal  theoretical  method  of  replacing  overgrowth  atoms 
on  a  substrate  with  Gaussian  mass  distributions  was  further  developed  for  those  cases 
where  the  effective  interatomic  potential  is  known.  The  technique,  applicable  to 
irregular-shaped  islands  or  films  of  finite  extent,  was  applied  to  a  simplified  model 
to  determine  preferred  orientation  relationships  from  calculated  film -substrate  inter¬ 
action  energies.  The  method  was  not  pursued  further,  however,  because  it  was  not 
sufficiently  adaptable  to  real  systems.  Several  other  possible  approaches  to  the 
theoretical  modeling  of  heteroepitaxial  systems  were  critically  reviewed,  including 
the  Frank-Van  der  Merwe  model,  a  Green's  function/Wannier-function  approach,  a 
contrived  potential -energy  model,  and  the  two-body  interactomic  potential  method.  It 
was  concluded  that  most  existing  theories  are  inadequate  for  application  to  real  systems. 

The  feasibility  of  a  molecular  orbital  development  of  the  heteroepitaxial  interface 
was  then  investigated.  However,  it  was  determined  infeasible  to  apply  this  technique 
in  a  manner  directly  relevant  to  heteroepitaxjr,  so  the  effort  was  terminated.  The  inter¬ 
atomic  potential  approach  to  heteroepitajgr  was  then  reinstigated,  with  the  goal  being  the 
computer  simulation  of  growth  of  Si  on  AfgOs*  Mechanical  stability  conditions  for  an 
AI2O3  lattice  modeled  with  two-body  potentials  were  investigated  and  determined  to  the 
depth  required  for  the  applications.  Computer  programming  of  the  AI2O3  lattice 
energy  and  elastic  constants  was  begun  for  use  in  determining  appropriate  empirical 
potentials  required  for  modeling  the  AI2O3  lattice,  a  major  requirement  for  modeling 
Si  growth  on  Ai203.  In  addition,  the  application  of  the  electron-on-network  theory  to 
the  problem  of  determing  surface  configurations  and  interfacial  binding  energies  in 
heteroepitaxial  systems  where  the  surface  structure  is  allowed  to  relax  was  investigated 
and  for  a  time  appeared  promising  for  the  real  systems  of  interest. 

During  the  past  six  months  the  theoretical  studies  have  been  devoted  to  three 
main  areas:  (1)  modeling  of  the  Sl/Ai203  system  by  means  of  interactomic  potentials 
and  computer  simulation;  (2)  use  of  the  electron-on-network  technique  to  calculate 
work  fiinctlons  and  surface  double-layer  potentials  of  monovalent  metals  (although  this 
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nK-tlu)(|  will  not  1,(.  pursued  lurther  for  |3ossible  uppliention  to  heteroepitaxiul  svstenisi- 
.  m  (.1  euleulutions  relating  the  mobility  of  Si  films  on  Af2()>t  to  stress  effects  arisinu  ’ 

alul  substrale  ‘'-m..!  contraction  of  fil^ 

The  modeling  of  the  A«2();]  lattice  has  been  earned  out  in  terms  of  Morse 
kn  lals;  only  anion-anion  and  cation-anion  potentials  have  been  employed  These 
.ve  been  chosen  to  conform  to  specific  constraints  related  t<.  the  physic'al  iattice 
linestipition  ol  surface  reconstruction  in  basal-plane  AfoO-,  has  also^been  initiated 

orkiJtld- ---'‘^tin.  si  nim 

ihn  fii,**'-  calcinated  changes  in  mobility  and  resistivity  in  the  plane  of 

n.si<^  ?  'r  system  caused  by  thermally  induced  stresses  and  the  pic/o- 

i  tsist.incr.  Clfect  hayc  been  found  in  excellent  aKrccincnt  with  experimental  results 

MgAf^O^  substiates  and  lor  a  variety  of  crystallograpliic  orientations.  ^ 

2.  SUllTASK  2.  l)Klk)Srn()N  STI’DIK.S  AND  FILM  PHKPARA  I’lON 

In  the  first  year  it  was  determined  that  tlu‘  electrical  propc.rties  of  undoni-d 
heteroejntaxial  Si  films  grown  on  various  orientations  of  AfoOo  (and  also 
Mg,A  20,)  l,y  the  CVl)  method  of  pyrolysis  of  Si.,,  are  doniinaUaf  b^y 

aoo  -md  /; nr’""*''’;  em-.2_or  Ixdow.  Kssentially  e(,uilalenl 

(  0  and  (llD-oriented  Si  films  were  grown  on  (0112)  and  (1014)  AfaOa  substrates  at 

the  nS  fiitodoplng  range  (~lor,()  C).  AC2()3  o’l’icntations'near 

strates  Ibr  I  Id’  Si ’  uwm' “’"a  ^  studies,  wert  utilized  as  sub- 

V  en/.  f,  1)  fc>i  hotel oepitaxy;  this  resulted  in  elect ron  mobilities  of  (500-700  cm^ 

\  -sec  for  carrier  concent ratimis  of  lolo  -  lt)17  cm-‘5.  which  exceeds  mobilities  ' 
obtained  on  either (0112)  or  (1014)  AM);}  substrates. 

y  ;vi,  I’C'calcd  the  strong  interrelationshi|)s  that  exist  among  the 

\aii oils  parameters  considered  as  nio.st  meaningful  for  optimizing  Si  growth  on  ^ 
insulatois.  Fvaluation  of  the  electrical  properties  of  Si  films  on  thosi'  A/.a)..  orienta- 
t  ons  iccogni zed  years  ago  as  providing  reflective  Si  overgrowths  has  deinonit rated 
111  "’“«t  be  maximized  for  the  substrate  orientation  chosen- 

(2)  differ  for  those  Af20;}  orientations  which  lead  to  the  same  Si  orientation;  (;i)  are 
l  epcndent  upon  reactor  geometry  and  gaseous  atmosphere;  and  (4)  should  he  optimized 
or  ic  particular  film  thickness  desired.  Preliminai-y  studies  of  Si  growth  by  Sill  i 
pyrolysis  n  roiluccil  pressures  (1  to  lOlorr,  luive  beei  very  eueouruKinr  udie^  1 
sine  e-eryslul  srowlh  cun  be  obluineil  over  a  fairly  iviilo  liunperature  rauie  ivhelr 

I  f.'  ’  subsirales.  InvesUenlion  of  Hie 

e  wll,  of  bi  films  on  AfjOa  and  MgA(2(l4  using  lie  as  Hie  einivth  alniospliere  11^  Hu 
carrier  gas  showed  that  epitaxial  growth  could  be  achieved,  but  the  conditions  for  good 
quality  growth  remain  to  be  established.  ^ 

MHr  offort  on  this  subtask  in  the  past  six  months  was  eoneeiitrated  on  eontinuina 

r  ttempts  to  optimize  the  Si  deposition  process  for  growth  on~(1120)  and  (t)lT2)  AM)^ 
surf,, CCS.  in  the  couse  of  this  work  the  effects  of  |»st-„ueleation  anilealllig  'lurin/lfi. 
deposition  process  on  ultimate  ftitii  properties  luive  Ix'en  examined,  but  no  sigiiifiealit 
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improvement  in  overall  (juality  of  Si  films  resulting  from  these  proeedures  has  been 
demonstrated.  The  effects  on  film  growth  of  gas-phase  etching  of  Af203  surfaces 
prior  to  deposition  have  been  evaluated  further;  there  is  some  indication  that  surface 
damage  may  not  be  the  primary  factor  in  determining  the  quality  of  Si  overgrowths, 
and  this  question  is  being  examined  iu  more  detail.  The  effect  of  cooled  reactor 
chamber  walls  on  the  Si  growth  process  has  also  been  investigated  further,  but  no 
significant  advantage  of  cooled  walls  in  the  vertical  reactor  systems  used  in  this  work 
has  been  observed.  Considerable  additional  study  has  been  made  of  the  growth  and 
properties  of  Si  films  on  Czochralski-grown  stoichiometric  spinel  (MgAf204).  This 
work  has  indicated  that  Si  films  with  electrical  properties  at  least  as  good  as  those 
grown  on  Af203  can  be  obtained  when  Ile-Il2  gas  mixtures  are  used  for  the  growth 
environment.  It  has  also  been  determined  that  autodoping  is  operative  in  the  Si/ 
MgAf204  system  (in  Ile-Il2  atmospheres)  at  approximately  the  same  temperatures 
as  for  Si/AfoO^. 

M  tj 

SUBTASK  3.  ANALYSIS  AND  PURIFICATION  OF  CVD  REACTANTS 

During  the  first  year,  techniques  of  gas  chromatography  were  developed  for 
analysis  of  the  reactants  used  for  semiconductor  heteroepitaxy  by  CVD.  A  general- 
purpose  gas-handling  system  was  constructed  for  the  highly  volatile  and  reactive  gases 
studied,  with  silicone  oil  and  polymer  columns  used  for  the  chromatography.  Several 
extraneous  impurity  peaks  were  observed  in  the  chromatograms  of  Sillq  samples; 
diborane  (B2II6)  was  tentatively  identified  as  a  significant  impurity  (~10  ppm),  although 
not  confirmed  by  mass  spectrometer  teclmiques.  Small  quantities  of  purified  Sillq, 
free  of  diborane,  were  prepared  by  successive  injections  in  the  chromatograph; 
quantities  were  too  small,  however,  for  use  in  laboratory  CVD  experiments. 

Beginning  in  the  second  year  of  the  program  samples  of  SiH4  of  trimethyl- 
gallium  (TMG)  used  for  Si  and  GaAs  CVD  expert  Jents  were  analyzed  for  impurity 
content  by  sensitive  mass  spectrometric  techniques.  Disilance  and  trimethylsilane, 
together  with  several  other  impurities  of  less  concern,  were  found  in  the  SiH,  samples. 
The  analyses  of  TMG  left  some  uncertainties  regarding  the  actual  impurity  levels,  and 
these  are  not  yet  fully  resolved. 

In  the  final  half  of  the  second  year  additional  reactant  analyses  have  been  carried 
out  by  mass  speetrometrie  techniques.  Significant  impurity  concentrations  in  some  of 
the  reactants  (especially  SiH4)  have  severely  limited  the  accuracy  of  the  study  of  the 
effects  of  deposition  parameters  on  Si  film  properties  on  several  occasions  during 
this  period.  Cooperative  efforts  with  vendors  for  preparation  and  analysis  of  improved- 
purity  reactants  have  continued. 

A  significant  new  study  of  the  chemistry  and  reaction  kinetics  of  CVD  processes 
used  for  growing  Si  and  GaAs  films  in  heteroepitaxial  systems  has  been  initiated; 
the  first  experiments  undertaken  are  directed  toward  determining  the  role  of  the 
Af203  surface  in  catalyzing  the  pyrolysis  of  SiH4.  Other  experiments  with  the 
trimethylgallium  (TMG)-arsine  (ASH3)  reaction  used  for  GaAs  growth  will  also  be 
performed. 
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4.  SUBTASK  4.  PREPARATION  AND  CHARACTERIZATION  OF  SUBSTRATES 


In  the  first  year  of  the  program  it  was  demonstrated  that  A£2(\‘5  surfaees  prepared 
by  mechanieal  polishing  teehniques  and  used  routinely  for  semieonduetor  heteroepitaxy 
typieally  have  severe  surfaee  and  subsurfaee  damage,  with  many  seralehes  often 
several  mierons  deep  yet  often  rendered  invisible  to  close  inspection  because  of  amor¬ 
phous  or  fine-grained  debris  embedded  in  the  scratches  in  the  final  polishing  stages. 
Some  improvement  in  mechanical  polishing  procedures  was  aehieved  in  terms  of  the 
density  and  depth  of  such  damage.  Gas-phase  elehing/polishing  procedures  using  SF,j 
and  various  fluorinated  haloearbons  in  Urn  13.'>0  to  1500  C_  temperature  range  produced 
essentially  serateh-free  surfaees  on  (0112)  and  ncar-(1120)  A£2f^<^  substrates.  Exten- 

st.v  guo— plitiisv  late  dattt  utVc  .lo  .i  fuiiCltwii  of  sliiil Ogl 

in  this  temperature  range. 


ry  sliiil  ogT  ic  i  i  ..n 


During  the  first  part  of  the  second  year  a  much  improved  technique  for  polishing 
iFdi  ij  A*  *ihl  wtfrf-iverf:  -mi  drfiL  iurtv'  loc 

obtained.  Gas-phase  ctehing/polishing  techniques  were  further  developed  and  exploited 
for  (1)  thinning  A^qO.^  substrates  to  thicknesses  the  order  of  1  mil;  (2)  evaluating  the 
effecis  of  prolonged  etching  on  (uili:),  ^uuul>,  and  (>0  assessing  die 

subsurface  damage  caused  by  various  mechanieal  polishing  procedures.  Routine 
characterization  of  substrate  surfaces  at  various  stages  of  preparation  continued  with 
techniques  of  X-ray  and  electron  diffraction  analysis,  optical  and  electron  microscopy, 
and  charged-particle  back-scattering  measurements. 

fliiP  fTKrul  TTi'ml  Hi'i*  iliblaiil"  h*i  h  »( kjr  iJ  tl»K  impfO-'rJ 

polishing  procedure  for  (1014)  ^  reprodiieible  proeess.  Evaluation  of  polishing 

metliods  for  MgA^204  surfaees  has  indieated  that  surfaee  fill-in  oeeurs  for  this  material 
Jnst  as  far  A42iJ3.  Gas-phasc  etching  eJtpcriinv-fils  tilth  MgAi^O^  iMipfaees  have  be^fi 
initiated.  Ion-beam  sputtering  techniques  have  been  developeafor  preparing  ultra-thin 
(-200A)  A?203  wafers  for  use  as  substrates  i;  the  in  situ  CVD  experiments  with  Si. 
tVafers  succtssfuliy  fnimied  lo -0.002  in.  oy  inechanical  po'iishlng  lecliniqaes  have  been 
subsequently  thinned  by  ion  etching  to  the  point  of  perforation  in  some  areas,  which 
iicaallji  in  ndjuining  j'U|^on]i  'of  IhiekD/'flfleft  nuiEnbU*  Uir  IrunandBaiKi  ptectreri  rrirrp>- 
scopy  as  applied  in  the  m  situ  experiments.  Typical  thinning  rates  are  0.1  to  2.0  pmVhr. 
Three  different  A£203  orientations  have  been  successfully  thinned  by  this  method  - 
(0001),  (1014),  and  (0112).  Considerable  study  of  properties  of  the  resulting  thinned 
substrates  has  been  carried  out  and  is  discussed  in  detail. 

5.  SUBTASK  5.  STUDIES  OF  IN  SITU  FILM  GROWTH  IN  THE  ELECTRON 
MICROSCOPE 


In  the  first  year  of  the  program  many  of  the  modifications  required  in  the  electron 
microscope  for  in  situ  observation  of  the  nucleation  and  early-stage  growth  of  CVD 
semiconductor  films  on  insulating  substrates  were  completed.  Provision  for  motion- 
picture  recording  of  film  growth  was  assembled  and  tested,  and  the  heated  specimen 
stage  was  installed  and  tested.  The  first  in  situ  PVD  e.xperimcnts  wore  also  carried 
out  near  the  end  of  the  first  year. 

Di  ring  the  first  six  months  of  the  second  year  a  series  of  electron  microscope 
modifications  and  tests  was  completed,  culminating  in  a  series  of  successful  PVD 

ItMlOr  UU  clrftftit'  lbir'rt>rT;|:x  Allfrt^Untai  via  lAh  n  tir  -lAfvf 
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carbon  substrate  and  a  sequence  of  micrographs  takefi  during  the  growth  process 

demonstrating  the  feasibility  of  performing  m  situ  nucleation  and  growth  studies  in  this 
equipment. 

A  transmission  phosphor  screen  (for  the  motion  picture  camera)  was  installed 
permitting  motion  picture  photography  which  does  not  interfere  with  the  normal  still* 
photography.  The  auxiliary  vacuum  pumping  system  for  the  specimen  chamber  was 
fabricated,  installed  and  tested.  The  basic  vacuum  system  of  the  microscope  itself 
addition  of  a  cooled  baffle,  by  polishing  numerous  O-ring  grooves 
and  by  thoroughly  cleaning  the  microscope  interior.  A  PVD  source  assembly  was 
fabricated,  installed  and  used  in  conjunction  with  the  specimen  heater  to  perform  the 
above  PVD  experiments.  Calculations  and  design  for  the  CVD  microchamber  were 
also  completed,  and  the  fabrication  of  the  microchamber  and  the  differential  pumping 
apertures  was  begun.  Work  was  also  begun  on  the  production  of  thin  (200  A)  A^oOo 
crystals  required  as  substrates  for  the  in  situ  CVD  studies.  ^ 

In  the  past  six  monlns  the  modified  design  of  the  CVD  microchamber  and 
associated  hardware  has  been  completed  and  fabrication  is  essentially  complete,  with 
testing  and  the  first  CVD  in  situ  experiments  now  planned  for  early  in  the  third-year 
program.  Numerous  in  situ  PVD  experiments  have  been  carried  out,  with  both  Ai  and 
Au  deposited  onto  amorphous  carbon  substrates  to  delineate  further  the  necessary 
techniques  and  experimental  problems  to  be  encountered  in  the  CVD  experiments. 
Kxtensive  oxidation  occurred  with  the  A£  films  when  grown  at  elevated  temperatures, 
and  unusual  nucleation  effects  were  observed  in  the  portions  of  the  Au  films  exposed  to 
the  electron  beam;  further  experiments  arc  in  progress. 

().  SUBTASK  (5.  EVALUATION  OF  FILM  PROPERTIES 

During  the  first  year,  the  routine  evaluation  of  film  properties  was  carried  out 
by  establishing  methods  of  X-ray  and  electron  diffraction  analysis,  metallographic 
analysis,  and  electrical  measurements  of  transport  properties.  In  addition,  a  new 
technique  for  evaluating  the  characteristics  of  the  interfacial  region  of  heteroepitaxial 
films  was  developed,  involving  measurement  of  the  photoelectron  emission  from 
monochromatically-illuminated  films  in  the  MIS  configuration  on  insulating  (viz. , 

A^gOg)  substrates.  Relatively  large  photocurrents  due  to  electron  transport  through 
thick  (-10  mils)  single-crystal  A£20g  substrates  were  measured  aj  a  function  of  photon 
energy.  Photoelectric  threshold  energies,  escape  length  (mean  free  path)  of  excited 
electrons,  and  band  bending  in  the  semiconductor  film  adjoining  the  interface  were 
determined  in  the  Si/AfgU'i  and  GaAs/AfgOg  systems.  Observation  of  the  energy 
spectrum  of  back-scattereJ  proton  or  alpha-particle  beams  injected  in  channeling 
directions  in  heteroepitaxial  semiconductor  films  w'as  also  investigated  as  a  means  of 
mepuring  the  density  and  the  location  of  structural  defects  in  the  films.  Experiments 
indicated  that  Si/insulator  films  have  less  imperfect  interfacial  regions  than  do  GaAs/ 
insulator  films.  Ihe  best  structures  of  those  examined  to  date  were  found  in  (100)  Si 
films  on  (0112)  Af20„  substrates  and  in  (lll)Si  films  grown  on  near-(1120)  AfoOo 
substrates.  ^  3 

In  the  first  part  01  the  second  year  the  study  of  the  effects  of  changes  in  deposition 
parameters  on  Si/AfgOg  film  properties  continued.  These  studies  have  led  to  consider¬ 
able  insight  into  the  factors  which  most  strongly  influence  film  quality;  identification  of 
conditions  for  optimized  film  growth  on  various  Af203  orientations  appeared  nearly 
complete.  The  importance  of  reactor  geometry  was  recognized  and  demonstrated. 
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The  extent  of  AH  autodoping  from  the  substrate  was  established,  and  appropriate 
annealing  procedures  for  minimizing  these  effects  w'ere  determined.  The  use  of  new 
substrate  orientations  for  Si  growth  led  to  epitaxial  films  of  as  good  or  better 
quality  than  those  previously  obtained. 

The  measurements  of  photoemission  of  electrons  from  heteroepitaxial  semiconductor 
films  and  of  the  transport  of  electrons  in  were  carried  further.  Work  functions 

of  additional  metals  were  determined,  and  the  mechanism  of  electron  transport  through 
the  insulator  was  studied  further.  Measurements  of  high-field  transport  properties  of 
Si  and  GaAs  heteroepitaxial  films  were  initiated. 

The  evaluation  work  in  the  past  six  months  has  included  (1)  study  of  the  variation 
of  the  electrical  properties  of  Si/A£2^3  v\’ith  temperature,  in  which  some  effects 
attributed  to  high  defect  densities  or  inhomogeneous  strains  were  observed;  (2)  an 
extensive  experimental  study  of  the  anisotropy  of  the  electrical  properties  in  Si/A£2^^3; 

(2)  evaluation  of  the  electrical  properties  of  Si  films  on  MgA£20j  ,  with  evidence  that 
n-type  films  with  mobilities  higher  than  those  obtained  in  the  Si/A^^O^  system  can  be 
obtained;  (4)  additional  measurements  of  the  electrical  properties  oT  Si/A£2<\3  fo 
establish  parametric  relationships  among  temperature,  growth  rate,  and  substrate 
orientation;  (5)  further  study  of  the  photoelectric  effects  observed  in  the  Si/A?2<)o 
system,  particularly  in  the  substrate  material,  including  verification  that  the 
observed  phenomena  do  result  from  photoinjection  of  electrons  from  metal  films  into 
the  A?203  f'nd  determination  of  the  work  functions  of  several  metals  and  the  heights  of 
the  metal-A£2^\3  interface  barriers,  and  ((5)  measurement  of  the  high-field  transport 
properties  of  neteroepitaxial  films. 

The  anisotropy  studies  have  been  concentrated  in  the  (221)Si/(l  122)  AH 2^^^  system 
and  in  the  (100)  Si/(0112)  Afi203  system.  The  two  Si  planes  are  basically  different  in 
that  the  anisotropy  in  the  (221)  plane  can  be  expected  for  any  Si  heteroepitaxial  system, 
while  that  in  the  (100)  Si  plane  results  from  the  anisotropic  thermal  contraction  of 
A^203  and  would  not  be  present,  e.g. ,  in  the  Si/MgA£2^^4  system.  The  mobility 
anisotropy  factor  A,  defined  as  the  ratio  of  the  difference  between  the  maximum  and 
minimum  values  of  carrier  mobility  in  the  plane  of  the  film  to  the  average  value  of  the 
mobility  in  tliat  plane,  was  found  to  be  about  39  percent  for  the  (221)  plane  and  about 
9  percent  for  the  (100)  plane.  Results  of  calculations  based  on  the  piezoresistance 
effect  in  Si  resulting  from  the  difference  in  the  thermal  expansion  coefficients  for  Si 
and  A^203  nre  given  and  agree  well  with  the  experimental  data.  The  calculations  and 
the  experimental  results  also  indicate  that  (221)  Si  probably  exhibits  higher  electron 
mobilities  man  omer  more  commonly  used  orieiuauons.  Tlie  lusults  of  the  anisotropy 
studies  are  highly  significant  to  preparation  of  device  structures  in  heteroepitaxial 
semiconductor  films. 

In  the  first  year  of  the  contract  apparatus  for  determining  carrier  lifetime  by 
C-V  measurement  in  MOS  structures  was  designed  and  constructed  and  tests  were 
begun.  A  special  MOS  structure  was  designed  for  measurement  of  channel  conductance, 
high-  and  low-field  transport  properties,  and  various  interface  characteristics  of 
heteroepitaxial  films.  Initial  attempts  to  fabricate  Schottky-barrier  diodes  in  Si/A£203 
films  as  a  means  of  evaluating  their  electrical  properties  were  not  successful. 
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tvoe  ofVET  Preliminary  design  of  a  Schottky-barrier 

t^e  ot  l-E  r  was  complelad,  with  experimental  FET  structures  to  be  fabricated  in 

in  SiVA^^O  n  W  ^  Preliminary  results  on  carrier  lifetime 

nn  h,\\\r  2.^3  Were  obtained,  after  initial  development  of  the  measurement  techm-’-ie 
on  bulk  single-crystal  Si  samples  was  completed  and  after  impurity  contamination 
problems  encountered  in  the  oxide  growth  process  were  solved. 

onf  mon^s  of  the  program  the  device-oriented  effort  has 

ce  tered  about  the  determination  of  carrier  lifetimes  using  the  MOS  pulsed-capacitance 

fabricate  a  Schottky-barrier  FET  in  GaAs/AfioOo.  ProcessJnc 
defects  have  been  encountered  in  both  of  these  investigations,  most  of  which  are  now 
tifp  analysis  was  made  to  evaluate  the  effect  of  the  impurity  redistribution  in 

the  Si  in  the  region  near  the  oxide  interface  (resulting  from  the  oxidation  process  used 
for  making  MOS  structures)  on  the  interpretation  of  lifetime  data  obtained  by  this^ 
technique.  Lifetime  measurements  in  both  Si  and  GaAs  films  as  a  function  of  film 
ickness,  carrier  concentration,  and  crystallographic  orientation  are  being  obtained. 
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APPENDIX  I 


WORK  FUNCTIONS  AND  SURFACE  DOUBLE  LAYER  POTENTIALS  OF 
MONOVALENT  METALS  FROM  A  NETWORK  MODEL* 

ABSTRACT 


The  model  described  in  this  paper  uses  an  electronic  wave  function  which  is 
defined  to  be  nonzero  only  along  the  linos  connecting  first  nearest  neighbors  in  the 
metallic  lattice.  The  electrons  are  assumed  to  move  freely  along  the  lines  between 
nearest  neighbors.  No  electron-electron  or  electron-nucleus  force  is  included  in  the 
model  calculations  (except  for  forces  arising  from  the  Pauli  exclusion  principle).  The 
work  function  is  defined  as  the  amount  of  energy  required  to  move  an  elootron  from  a 
point  slightly  inside  the  crystal  to  a  point  slightly  outside.  The  contribution  of  the 
electronic  double  layer  is  included  in  the  calculation  of  the  work  function  as  well  as 
the  dependence  of  the  double  layer  potential  on  the  surface  geometry.  Surface  states, 
where  the  electron  is  localized  in  the  neighborhood  of  the  face  of  the  crystal,  are 
found  to  have  energies  sufficiently  above  the  Fermi  level  to  eliminate  the  possibility 
that  they  make  any  contribution  to  the  double  layer  potential  for  the  case  of  the  (100) 
crystal  plane.  Consequently,  surface  states  have  been  ignored  in  all  the  calculations. 
The  surface  double  layer  is  assumed  to  be  caused  by  the  presence  of  a  finite  potential 
barrier  at  the  surface  of  the  crystal.  Bulk  electronic  wave  functions  can  penetrate 
this  barrier  and  decay  exponentially  outside  the  crystal.  The  only  parameters 
required  by  the  model  are  the  nearest  neighbor  distance  for  the  lattice  and  the  height 
of  the  potential  barrier  at  the  surface.  The  former  quantity  is  fixed  by  the  lattice 
structure  (body  centered  cubic  for  the  alkali  metals)  and  by  the  density,  while  the 
latter  quantity  can  be  adjusted  to  give  the  best  agreement  between  the  model  calcula¬ 
tions  and  experiment.  For  the  alkali  metals,  lithium  through  sodium,  the  best  value 
of  the  barrier  height  is  about  50  percent  of  the  sum  of  the  ionization  potential  energy, 
the  heat  of  vaporization,  and  the  calculated  Fermi  level  for  the  corresponding  metal. 
In  addition,  the  value  of  the  double  layer  potential  for  sodium  agrees  very  well  with  a 
more  sophisticated  calculation  by  Bardeen  and  is  reasonably  close  to  the  experimental 
measurement. 


♦Preprint  of  a  paiper  by  P.  K.  Rawlings  and  II.  Reiss,  Department  of  Chemistry, 
UCLA,  submitted  for  publication.  This  work  supported  at  UCIJ\  in  part  by  ARPA 
Order  No.  1585,  Contract  No.  DAAH01-70-C-1311,  through  North  American 
Rockwell  Corporation,  and  in  part  by  NSF  Grant  No.  GP-28722X. 
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1.  INTRODUCTION 


The  free  eleetron  molecular  orbital  method  (FEMO)  has  been  successfully 
applied  by  Kuhn  (Ref  1)  and  by  Rudenberg  and  Scherr  (Ref  2)  to  the  treatment  of 
n -eleetron  systems  in  conjugated  molecular  systems.  As  in  the  case  of  Huckel 
theory  (Ref  3),  the  free  electron  orbital  method  is  directed  towards  the  development 
of  a  "relative"  internally  consistent  theory  for  these  complicated  systems.  In  the 
case  of  Ruckel  theory,  parameters  such  as  overlap  integrals  are  determined  by 
experiment  in  one  molecule  and  are  assumed  to  have  the  same  magnitude  in  others. 
Strictly  speaking,  in  the  FEMO  method  there  are  no  adjustable  parameters  (unless 
the  distances  between  atoms  are  allowed  to  depart  from  the  values  actually  observed 
in  the  molecule).  Nevertheless,  both  methods  provide  a  reasonable  semi-quantitative 
picture  of  the  electronic  state  of  affairs  within  r -systems,  and  have  been  useful  for 
making  semi -quantitative  estimates  concerning  the  chemical  and  spectroscopic 
properties  of  ir-systems. 

Recently  Montroll  and  his  coworkers  (Ref  4)  have  expanded  the  FEMO  method 
to  include  the  treatment  of  local  potentials  due  to  individual  atoms.  The  theory  of 
course  remains  highly  schematic,  as  well  as  parameterized,  but  it  does  allow 
relatively  simple  analytical  solutions  to  be  obtained  for  previously  complicated 
problems.  In  this  way.  Montroll  and  his  collaborators  have  been  able  to  explore  the 
qualitative  features  of  many  intei'csting  solid  state  phenomena  including  the  electronic 
properties  of  defects  and  surfaces.  Since  the  original  FEMO  method  proved  useful 
in  the  semi-quantitative  discussion  of  the  electronic  properties  of  molecules,  it  is 
interesting  to  see  whether  both  the  FEMO  and  Montroll's  modification  of  it  can  be 
used  in  the  semi -quantitative  study  of  solids.  We  have  in  mind  not  only  the  examina¬ 
tion  of  schematic  relationships  (Mon troll  and  his  collaborators  have  already  done 
that),  but  also  the  development  of  an  internally  consistent  parameterized  theory  in 
the  same  spirit  as  Iliickel  theory. 

In  order  to  study  this  question,  we  have  elected  to  apply  the  network  theory  to 
the  estimation  of  the  work  functions  of  monovalent  metals,  especially  the  alkalis. 

Our  principal  objective  has  been  to  determine  the  work  function  and  its  variation 
with  crystal  orientation.  The  approach  is  as  simple  as  possible,  involving  a  single 
adjustable  parameter,  namely  the  difference  between  the  energy  of  a  free  electron 
in  vacuum  and  the  lowest  level  of  the  conduction  band.  This  parameter  is  invariant 
to  crystal  orientation  and  could  in  principle  be  determined  from  measurements  on 
a  single  crystal  plane.  All  refinements,  such  as  exchange  and  correlation  effects, 
are  explicitly  ignored  but  should  be  implicitly  involved  within  the  barrier  height 
which  constitutes  the  adjustable  parameter  mentioned  iibovc.  The  detailed  explicit 
treatment  of  these  additional  effects  would  be  inconsistent  with  the  severe  approxima¬ 
tions  already  contained  within  the  theory.  One  would  hope,  however,  that  the  network 
model  would  give  a  reasonable  semi-quantitative  account  of  the  variation  of  work 
function  with  crystal  orientation,  once  the  adjustable  parameter  is  determined,  by 
using  the  experimental  data  available  for  one  crystal  plane.  Furthermore,  one  would 
hope  that  there  would  be  an  orderly  relationship  between  the  predictions  of  the  theory 
for  one  alkali  metal  as  compared  to  another. 
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2.  THE  NETWORK  MODEL 

Cur  treatment  actually  makes  use  of  the  FEMO  method  coupled  with  the  Green's 
function  approach  developed  by  Montroll  for  the  treatment  of  defects  and  surfaces  in 
connection  with  his  more  elaborate  network,  model.  The  one  electron  wave  function  is 
chosen  to  be  nonzero  only  along  the  lines  connecting  the  first  nearest  neighbors  in  the 
metallic  lattice.  This  is  equivalent  to  requiring  that  the  potential  be  infinite  every¬ 
where  within  the  crystal  except  on  these  lines.  The  electrons  are  assumed  to  move 
freely,  i.e. ,  wc  do  not  associate  potential  wells  (as  in  the  Montroll  modification)  with 
the  individual  atoms.  To  this  extent,  our  model  is  more  like  the  original  FEMO 
method.  However,  we  shall  have  to  make  use  of  the  Bloch  theorem  (Ref  5)  in  order 
to  deal  with  the  crystal  and  the  Green's  function  approach  will  have  to  be  used  to 
treat  the  extended  defect  represented  by  the  surface.  The  crystal  symmetry  enters 
the  problem  through  the  nodal  conditions,  Eq  (11)  and  (12)  below,  and  through  the 
Bloch  theorem. 

The  work  function  is  defined  as  the  amount  of  energy  required  to  move  an 
electron  from  the  point  slightly  inside  the  crystal  to  a  point  slightly  outside,  and  is 
effected  by  the  electronic  double  layer  (Ref  6)  which  in  turn  depends  on  crystal  surface. 
One  of  our  main  tasks  will  involve  calculating  the  strength  of  this  double  layer. 

As  long  as  an  electron  remains  on  a  line  joining  two  nearest  neighbor  atoms 
inside  the  lattice,  it  has  no  force  acting  on  it,  and  its  potential  energy  may  be  set 
equal  to  zero.  Wc  assume  that  an  abrupt  change  in  potential  occurs  as  the  electron 
moves  on  to  a  dangling  surface  bond  where  it  experiences  a  constant  potential  Vg 
(Vr  >0),  The  surface  barrier  Vg  is  the  only  parameter  in  the  theory  and  will  have 
to  be  determined  through  reference  to  some  single  experiment.  In  fact,  the  form  of 
the  potential  barrier  in  the  real  crystal  is  not  abrupt,  but  since  we  arc  searching  for 
the  simplest  possible  parameterized  theory  (and  also  because  of  the  severe  approxima¬ 
tion  already  contained  within  the  network  model  itself),  we  assume  that  it  is  abrupt. 

The  dangling  bonds  for  surface  atoms  have  the  same  direction  as  if  they  were  joined 
to  another  lattice  site;  but  instead  of  being  of  length  R  (where  R  is  the  nearest 
neighbor  distance),  they  extend  all  the  way  to  infinity.  The  wave  equation  for  this 
system  may  be  expressed  in  the  following  form; 


-  {a  -k^)4>=  0. 


(8) 


in  which  ,  k^  =  and  4^  is  the  electronic  wave  function,  m  the  mass  of 

the  electron,  E  the  energy  of  the  state  described  by  4^,  V  the  potential,  and  h  is 
Planck's  constant  divided  by  2Tr.  Clearly  V  (and  a)  equals  zero  for  positions  along 
bonds  in  the  bulk,  while  V  =  Vg  on  the  dangling  bonds.  The  solutions  of  Eq  (8  )  for 
the  surface  and  bulk  bonds,  respectively,  are; 


4jg  =  a^  exp  { -  (a^  -k^)^^^x  )  , 


(9) 


=  a  cos  {kx+6n)}. 
S  B 


(10) 
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Each  is  defined  only  on  a  particular  bond  and  is  zero  elsewhere.  Where  several 
bonds  meet  at  a  lattice  point,  the  wave  fimctions  associated  with  them  must  be  equal. 
In  addition,  there  are  requirements  associated  with  the  conservation  of  mcvnentum. 
These  restrictions  may  be  expressed  as  follows:  (Ref  2) 

"  4^3  (o)  ^  •  -■  4^^  (o)  . 


(12) 


In  these  equations  n  is  the  number  of  bonds  meeting  at  a  particular  lattice  site. 

If  the  quantity  4)  (2.m,n)  is  the  value  of  the  wave  function  at  the  lattice  site 
(2,m,n),  Eq  (11-12)  require,  for  a  body  centered  cubic  lattice,  the  following  relationship 


8cos(kR)  4)  (2,m,n)  =  4)( 2+1, m+l,n+l)  +  4> (2+l,ni-l,n+l)  +  4>(2-l,m+l,n+l) 
+  4)  (2-l,m-l,n+l)  +4>(^+l,m+l,n-l)  +  4>(^+l,ni-l,n-l)  +  4>  (2-l,m+l,n-l) 


+  4>  (^-1*  ni-l,n-l)  .  (13) 

In  addition,  for  the  bulk  crystal  we  make  use  of  periodic  boundary  conditions: 

4>(2m,n)  =  4>(2,N,m,n)  =  4»(2,m+N,n)  =  4)(2,m,n+N).  (14) 

A  possible  solution  to  Eq  (13)  and  (14),  which  also  satisfies  Bloch's  theorem,  is; 

4)(2,n,n)  "  A  expl^^  (^Sj +ms2 +nsg)}  ,  (15) 

where  Sj,  s„,  Sg  =  1,2,3, . .  .N.  Substitution  of  Eq  (15)  into  (13),  with  appropriate 
collection  ofterms,  yields: 

2tts-  ^^®3 

cos  (kR)  =  cos  (— — )  cos  — )  cos  (— jj”)  =  E.  (16) 


Since  the  energy  of  an  electron  depends  on  k,  and  Sj,  S2,  and  Sg  are  indices 
denoting  possible  states  of  the  system,  Eq  (16)  determines  the  density  of  states 
for  network.  The  density  of  sets  (s^,  Sg,  Sg),  satisfying  Eq  (16)  in  the  range  dE, 
has  been  derived,  in  another  context,  by  Jelitto  (Ref  7)  (see  Figure  45).  The  density^ 
of-states  function  is  symmetric  about  cos  (kR)  =  0,  and  it  follows  that  the  Fermi  level 
of  the  network  corresponds  to  an  energy  level  such  that  kR  =  tt/2.  Therefore,  the 
Fermi  level  is  given  by: 


^F 


h 


2 


32mR^  • 


(17) 
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Figure  45.  The  Det 
as  a  Functio 
Momentum 


In  order  to  deal  with  surface  states,  we  assume  that  the  network  is  infinite  in 
both  the  X  and  y  directions,  but  has  a  finite  thickness  in  the  z  direction.  It  is  possible 
to  write  a  general  relationship  (analogous  to  Eq  (13) )  which  applies  to  either  bulk  or 
surface  sites.  Thus,  we  have: 

8cos  (kR)  (}>(£, m,n)  -  <t>(f+l,m+l,n+l)  -  <j>(£+l, m-l,n+l)  -  (t>(£-l,m+l,n+l) 

-  <j>(£-l,m-l,n+l)  -  <t>(2+l,m+l,n-l)  -  <|>(£+l,m-l,n-l)  -  4>(^-l,m+l,n-l) 


-  4)(£-l,m-l,n-l)  =  H  (£,m,n)  ,  (18) 

(t>(£+N,m,n)  -  4>(fi,m+N,n)  (|>(£,  m,n)  ,  (19) 


where 


II  (£,m,n)  =  i+^n  N  ^  [4  cos  (kR)  -4A  (kR)  ](|>(£,m,n) 

’  z 

I‘t>(«+l,m+l,Nz)  +  4>(«+l,m-l,Nz)  +  ‘t>(£-l,m+l, N^)  +  <j>(£-l, m-1, N^)  ] 


[(|>(fH,m+l,l)  +  <|)(£+l,m-l,l)  +  «j>(£-l,m+l,l)  + 


<j>(£-l,m-l,l) 


(kR)2 


sin  (kR) 


and  the  6's  are  Kronecker  delta  functions.  An  equation  of  this  general  form  was 
introduced  by  Montroll  (Ref  3)  although  he  did  not  specifically  treat  the  case  in  which 
an  abrupt  barrier  of  finite  height  V„  was  located  at  the  surface.  We  assume  that  the 
z-direction  is  the  [OOl]  direction  ana  that  the  exposed  surfaces  are  (001)  surfaces 
having  coordinates  z  =  1  and  z  =  Ng,  respectively. 

Because  of  the  finite  extent  of  the  crystal  in  the  z  direction,  we  can  no  longer 
demand  that  <t>  satisfy  the  Bloch  theorem  in  that  direction.  However,  this  requirement 
may  still  be  applied  to  the  x  and  y  directions.  A  solution  to  Eq  (18),  which  satisfies 
the  Bloch  theorem  in  these  directions,  is; 


(}>(£, m,n)  =  A  exp 


2  ffj 

N 


[Sj,t  +S2,  ml 


4>(n) 


Sf,  Sg  -  1,  2,  3,  . ..  N  . 


(20) 


If  wc  substitute  this  into  Eq  (18)  and  collect  terms,  the  result  is 
2Fj3  (j>(n)  -  (t>(nHl)  -  (t>(n-l)  =  llg  (n)  , 


(21) 
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in  which 


HbW'  -‘n.N, 


z 

2ns 


2ns, 


1  2 
Fp  =  cos  (kR)/[cos  (  )  cos  (-jj — ]  , 


2nSj 

n  =  [cos  (kR)  -A  (kR)]/[cos  (— cos  (— — )]  . 

A  difference  equation  having  the  same  form  as  Eq  (21)  has  been  solved  by  Montroll 
(Ref  4)  using  an  appropriate  Green's  function,  g(n).  hi  terms  of  g,  this  solution  may 

be  expressed  as; 

N 

(|)(n)  =  (g(j-n)Hg(j)  ,  (22) 

j=l 


where  the  Green's  function  is  defined  as: 


exp 


I  2nisn  I 

l-n^l 


<2ns, 
Fg  -cosf^ 


^23) 


For  the  |F„|  >  1,  4'(n)  decays  exponentially  as  n  refers  to  sites  further  from  the 
surface,  i.  e. ,  deeper  in  the  bulk  of  the  ciystal.  Therefore,  this  range  of  Fg  is 
the  range  for  surface  states;  and  by  substitution  of  the  solution  in  this  range  back 
into  Eq  (21),  it  may  be  shown  that  the  following  relation  holds: 


2ns. 


2ns, 


[cos^  (kR)-A^  (kR)  =  cos  (-^j-^)  cos  (— 


(24) 


This  equation  is,  for  surface  states,  the  equivalent  of  Eq  (16)  for  bulk  states  and, 
like  Eq  (16),  determines  the  density  of  states. 

The  bulk  wave  functions  may  be  normalized  as  follows.  First,  the  square  of 
the  wave  function  (square  of  the  modulus)  is  Integrated  over  each  individual  bond, 
yielding  a  quantity  proportional  to  the  fraction  of  electrons  residing  on  the  bond. 
This  fraction  is  then  summed  over  all  bonds  and  the  sum  is  set  equal  to  unity.  The 
fact  that  each  state  is  occupied  by  two  electrons  is  then  taken  into  account  at  a  later 
step  by  multiplying  the  density  of  states  function,  Eq  (29),  by  2.  The  integral  for 
the  square  of  the  wave  function  on  a  dangling  surface  bond  is,  according  to  Eq  (  9  ): 

2  r  , 

D  =|4'g(o)r  J  exp(-2px)dx=  - - .  (21 

o 


lU 


9  O  1/2 

In  the  above  equation,  p  =  (o^  -  and  |4'g  (o)j  is  the  value  of  the  wave  function 

at  the  surface  node.  The  integral  for  the  square  of  the  wave  function  on  a  bond  con¬ 
necting  two  nearest  neighbor  lattice  sites  is,  from  Eq  (10): 


D 


B 


(kx  +  6j^) 


2 


R  + 


sin  [2  (kR+  sin  [2  6g] 
2K 


(26) 


where  R  is  again  the  distance  between  nearest  neighbors.  It  should  be  noted  that 
although  both  the  above  integrals  depend  on  the  kinetic  energy  of  the  electron,  Dg  is 
also  a  function  of  the  surface  barrier  height,  while  D3  is  a  function  of  the  distance  R. 
In  a  body  centered  cubic  lattice,  each  site  has  eight  first  nearest  neighbors  and  there¬ 
fore  ei^t  radiating  bonds.  However,  for  interior  sites,  each  of  these  eight  bonds  is 
shared  by  one  other  lattice  site  and  thus  each  interior  site  has  only  four  bonds  on  the 
average. 

In  accordance  with  the  prescription  laid  out  above,  we  therefore  write  the 
following  expression: 

4D3  la^l^  N^  +  .  .  .  +0(N^  =  1.  (27) 

where  is  the  number  of  sites  in  the  lattice  and  is  proportional  to  the  number  of 
surface  sites.  All  terms  on  the  left  of  Eq  (27).  except  the  first,  are  of  order  or 
less;  and  since  the  size  of  the  network  can  be  thought  of  as  being  arbitrarily  large, 
only  the  cubic  term  needs  to  be  considered.  Consequently,  the  following  expression 
is  obtained: 


_ 

4DjN 


(28) 


It  turns  out  that  for  reasonable  values  of  the  barrier  height,  Vb'  surface  states 
on  the  (001)  plane  Lave  energies  above  the  Fermi  level  and  are  consequently 
unoccupied.  As  a  result,  surface  states  do  not  influence  the  surface  double  layer 
potential  or  the  woric  function,  at  least  not  within  the  requirements  of  the  network 
model.  Thus,  the  properties  of  the  surface  states  need  not  be  included  in  the  cal¬ 
culations  for  the  double  layer  and  work  function  which  follow.  We  can  therefore  make 
full  use  of  the  density  of  states  derivable  by  the  method  of  Jelitto  (Ref  7). 

This  density  of  states  in  reduced  momentum  sparse  is: 

dN 

p(kR)  -  Cj^sin  (kR)-^,  (29) 

where  dNg/dE  has  been  expressed  rigorously  by  Jelitto  in  terms  of  elliptic  integrals. 
It  may  also  be  approximated  using  more  elementary  functions  as  follows: 


dE 


2[1-  |in  1^16.679  +  3.6364  •  |E|  +  2.4880  •  |E|^j  ; 
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an  expression  t;ood  to  four  significant  figures  over  the  range  0.  OOfj  £  |K|  <1.0. 

K  eos  (kit)  and  Cj^  is  a  normalizing  constant  such  that  the  total  number  of  cjuantitics 
in  the  region  0  £  kit  <  tt  is  N^,  i.  c.  , 


f 


p  (kit)  d(kR)  N 


liceause  the  expression  for  p(kH)  is  rather  complicated,  the  above  integral  is  evaluated 
most  conveniently  by  numerical  means.  Unfortunately,  the  value  oi  dNjj/dU  is  infinite 
for  IK  I  0,  and  the  numerical  integration  cannot  be  performed  over  the  entire  integral 
even  though  the  integral  converges.  'I’lie  density  of  states  function,  p(kl{),  is  sj'm- 
metrie  on  the  integral  about  the  point  kit  tt/2,  and  eonsequently  the  numerical  integra¬ 
tions  need  only  be  performed  from  kR  =  0  to  kR  -  u,  where  u  <  tt/2.  There  will  be 
some  uncertainty  about  the  exact  value  of  depending  on  the  choice  of  upper  limit 
for  the  integral.  For  the  ealeulations  in  this  paper,  u  has  been  set  equal  to  0.498.1it. 
This  particular  value  was  chosen  for  two  reasons:  (1)  the  approximation  for  dNjj/dK 
is  known  to  be  reasonably  accurate  over  this  range  of  integration,  and  (2)  the  value 
of  differs  by  only  a  few  percent  from  the  value  assumed  when  the  upper  limit  on  the 
integral  is  changed  from  0.4985Tr  to  0.49999Tr. 


2.  KSTIMAllON  OF  SURFACE  CHARGE 


The  charge  on  each  dangling  surface  bond  may  be  calculated  in  a  straightforward 
manner,  being  proportional  to  Dg.  Multiplying  this  ratio  by  2  to  account  for  electron 
spin  and  finally  integrating  over  all  states  (all  possible  values  of  momentum)  gives 
the  fraction  of  the  total  number  of  electrons  residing  on  a  surface  bo^'d.  We  call  this 
quantity  G(oR): 


eos  (6jj) 


G  (oR) 


i 


-  TT, 


/2 


2p(kR)D  d(kR) 
s 


‘'o 


I  -  (kR)^ 

‘'0  sinr2(kR  +  6„ 
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/p(kH)\ 


1  + 


[2(kR  +  6p)]  -  sin  [2  6jj] 


.d(kR).  (30) 


2kH 


In  itiis  equation,  Ug  and  Dq  are  defined  by  Eq  (25)  and  (26),  respectively,  and: 


,r  V  GOS  (kR)  -  1 

tan  (Or)  - 


1'^  (o)r 


-  |a/eos^ 


These  Iv/o  relations  are  derived  from  the  requirement  that  the  square  of  the  modulus 
of  the  standing  waves  represented  by  Lq  (9-10)  must  be  the  same  at  all  lattice  sites  for 
wave  functions  associated  with  the  bulk  energy  states. 
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IBecause  of  the  diffieultics  already  discussed  in  connection  with  p(kR),  the  actual 
upper  limit  of  integration  used  in  these  calculations  was  0.  4  9  8  57t  instead  of  tt/2.  The 
bulk  wave  functions  are  assumed  to  differ  only  by  a  phase  factor  from  one  lattice  site 
to  another  even  at  the  surface.  Because  G(oR)  is  normally  a  few  tenths  of  an  electron 
or  less,  and  because  the  conductivity  of  the  metal  is  high,  all  the  positive  charge, 
balancing  the  negative  charge  located  outside  the  lattice,  is  assumed  to  reside  on 
the  crystal  surface.  The  total  positive  surface  charge  is  therefore: 


a 

S 


G(oR) 


e. 


(31) 


where  Ug  is  the  number  of  dangling  surface  bonds  and  (Ng/Ag)  is  the  number  of  sites 
per  unit  area,  while  e  is  the  charge  on  an  electron. 


4.  CALCULATION  OF  THE  DOUBLE  LAYER  POTENTIAL 


The  work  required  to  move  an  electron  from  the  surface  to  some  point  outside 
the  crystal  may  be  calculated  using  Gauss'  theorem.  Since  the  distribution  of  electric 
charge  is  uniform  parallel  to  the  cx’ystal  surface,  the  field  acting  on  any  electron  is 
perpendicular  to  the  surface  and  has  an  intensity: 


E(x)  =  (32) 

o 

where  -  1.811  x  10“^  coulombs/Newton/rnctcr^  and  o^lx)  is  the  net  amount  of 
positive  charge  lying  between  the  surface  and  the  point  x  on  the  dangling  bonds.  The 
charge  (r(x)  is  evaluated  from  the  expression: 


a(x) 


-4-  / 

2N'*  4 


cos  (6p)  exp  [-2(3x] 

J72 


_  [(oR)^  -  (kR)^] _ 

sin  [2  (kR  +  6j^^  -  sin  [26^] 


1  + 


2kR 


(33) 


in  which  we  have  combined  Eq  (30)  with  (9).  The  work  performed  against  the  double 
layer  is  then: 


E(f)  d£ 


(34) 


In  this  equation,  df  represents  distance  along  the  direction  of  the  surface  bond,  and 
Kq  is  a  factor  such  that  K^df  is  the  distance  normal  to  the  surface.  The  integral  in 
Eq  (34)  is  most  conveniently  evaluated  in  terms  of  the  reduced  distance  i/R,  and  the 
exact  value  of  the  upper  limit  3  not  very  important  bceuase  essentially  all  of  the 
contribution  to  the  double  layer  potential  occurs  within  less  than  five  times  the  nearest 
neighbor  spacing  from  the  surface.  Thus,  within  the  limitations  of  this  model, 
changes  in  the  double  layer  potential  from  one  surface  to  another  are  determined  by 
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the  variations  in  three  factors:  (1)  the  number  of  lattice  sites  jier  unit  area  of  surface, 
(2)  the  number  of  dangling  surface  bonds  at  each  site,  and  (3)  the  geometric  factor  lyj. 

5.  THE  WORK  FUNCl’lON 

As  indicated  earlier,  the  work  function  is  defined  as  work  necessary  to  move  an 
electron  from  the  point  just  inside  the  lattice  to  a  point  just  outside.  The  energy 
required  for  this  process  depends  on  the  geometry  of  the  surface  tlirough  which  tlie 
electron  passes  because  of  the  contribution  of  the  double  layer.  In  keeping  witli  the 
simplicity  of  the  network  model,  the  work  function  VVp  is  defined  as: 


W 


E 


'^11  ■ 


W 


D 


(35) 


whereas  before  Vj^  is  the  surface  barrier,  |i|.  is  the  Fermi  level,  and  now  we  denote 
the  double  layer  potential  by  VVy.  The  first  two  quantities  on  the  right  of  IJq  (35)  are 
independent  of  surface  orientation,  and  the  variation  of  work  functions  with  surface  is 
therefore  due  to  VVjj. 

The  double  layer  potential,  Eq  (34),  can  be  expressed  most  conveniently  in  the 
following  form: 

VVjj(oR)  A^  l(oH)  (3C) 

where 


and  is  only  a  function  of  the  nearest  neighbor  distance  and  surface  geometry,  while 
1  (oR)  is  an  integral  which  depends  only  on  the  praduct  qR.  The  values  of  1  (qR)  for 
various  oR  have  been  listed  in  Table  XIX.  In  Table  XX,  the  densities  of  the  alkali 
metals  and  their  first  nearest  neighbor  distances  arc  given,  as  well  as  the  values  of 
Al  for  the  (001)  and  (Oil)  planes. 

The  surface  barrier  Vjj  may  be  estimated  by  the  following  expression: 

Vj3  «  I  +  H  +  (37) 

Here  I  is  the  ionization  potential  for  an  atom  in  gas  phase,  H  is  the  heat  of  sublima¬ 
tion  per  atom,  and  pjp  the  Fermi  level  which  can  be  estimated  from  Eq  (17).  Equation 
(37)  approximates  an  equation  given  by  Bardeen  (Ref.  8)  and  has  its  origin  in  comparing 
certain  electronic  integrals  which  occur  in  the  expression  for  Vji  with  integrals  which 
also  occur  in  expressions  for  I  and  H.  I'he  barrier  height  may  also  be  estimated  by 
adjusting  Vg  so  that  the  value  predicted  by  the  network  theory  agrees  with  experiment. 
This  was  done  for  the  alkali  metals  using  the  work  function  calculated  for  the  (001) 
plane.  The  experimental  heats  of  vaporization,  ionization  potentials,  and  work  func¬ 
tions  have  been  compiled  from  various  sources  by  Wigner  and  Bardeen  (Ref.  (5). 


115 


TABLE  XIX 


all 

I  (all) 

0.6 

0. 10553 

0.  7 

0.05890 

0.8 

0.038  4  4 

0.  9 

0.02879 

1.0 

0.02215 

1.  1 

0.01765 

1.2 

0.01447 

1.  3 

0.01207 

1.4 

0.01025 

1.  5 

0.008t2 

1.6 

0.00767 

TABLE  XX 


Element 

Density 

(gm/cm^) 

Nearest 

Neighbor 

Distance 

(A) 

^L 

(001) 

1  > 

Lithium 

0.534 

3.  038 

25.82 

18.26 

Sodium 

1.007 

3.  666 

21.  39 

15.  12 

Potassium 

0.870 

4.  594 

17.08 

12.07 

Rubidium 

1.  532 

4.937 

15.88 

11.23 

Cesium 

1.873 

5.  349 

14.67 

10.  37 

In  Table  XXI,  we  compare  the  barrier  heights  estimated  from  Eq  (37)  with  those 
obtained  through  adjustment  of  Vjj  to  make  the  work  function  calculated  by  the  network 
model  agree  with  experiment.  The  best  fit  value  obtained  from  the  network  model  is 
always  about  one-half  the  value  derived  from  Eq  (37).  On  the  other  hand,  Eq  (37)  is 
itself  by  no  means  precise.  Considering  the  simplicity  of  the  network  model,  the 
results  are  not  very  disparate.  In  fact,  in  the  case  of  both  Hiickel  theory  and  the 
FEMO  method  applied  to  molecules,  tlie  same  order  of  magnitude  is  observed  for  the 
mismatch  between  theory  and  experiment. 

Also  included  in  Table  XXI  are  double  layer  potentials  for  two  different  surfaces. 
Bardeen,  (Ref.  8)  using  tlie  free  electron  model,  calculates  that  the  surface  double 
layer  for  the  (110)  plane  of  sodium  is  0.4  volt.  The  network  model  predicts  a  value  of 
0.  37  volt  (using  the  value  of  Vr  which  brings  the  theoretical  and  experimental  work 
functions  into  register)  as  compared  with  an  experimental  value  (Ref.  8)  of  0. 15  volt. 

From  these  results,  it  appears  as  though  the  simple  network  model  possesses 
about  the  same  quantitative  degree  of  validity  with  respect  to  surface  slate  problems 
that  the  FEMO  model  possesses  with  respect  to  molecules.  This  suggests  that  it  may 
be  used  to  make  ballpark  estimates  in  a  variety  of  solid  state  situations. 


TABLE  XXI 


Element 

oR 

Vd 

(volts) 

(volts) 

*^(001) 

(volts) 

*^(001) 

(volts) 

'^(001) 

(volts) 

>11) 

(volts) 

I+Il+lij, 

(volts) 

'u 

I+H+|ip, 

Lithium 

0.  979 

0.  3903 

1.018 

2.310 

0.  575 

2.479 

0.400 

8.00 

0.484 

Sodium 

1. 107 

0. 3420 

0.  099 

2.355 

0.  372 

2.404 

0.203 

0.95 

0.493 

Potassium 

1.271 

0. 2870 

0.445 

2.214 

0.217 

2.  230 

0. 153 

5.  72 

0.503 

Rubidium 

1.  334 

0.2743 

0.  385 

2.177 

0. 181 

2.  230 

0. 128 

5.42 

0.500 

Cesium 

1.  344 

0.  2373 

0.  328 

1.881 

0.104 

1.  929 

0.  no 

5.03 

0.472 
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